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The  structural  analysis  and  design  of  concrete  structures  is  often  based  on  large-scale 
finite  element  computations.  A  key  component  of  the  finite  element  method  is  the  constitu¬ 
tive  model  chat  is  usually  selected  and  calibrated  based  on  a  few  simple  tests.  The  primary 
reason  for  this  is  the  lack  of  understanding  of  the  complex  response  features  of  concrete 
and  numerical  difficulties  encountered  when  attempting  to  model  these  features.  Also,  the  >\ 
error  in  response  predictions  due  to  an  inappropriate  constitutive  model  is  difficult  to 
define  in  a  complex  large-scale  structural  analysis  problem. 
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19.  ABSTRACT  (Continued). 

The  method  of  evaluation,  as  developed  herein,  consists  of  the  following  steps: 

1.  The  design  and  execution  of  a  series  of  material  properties  tests  which  provide 
data  sufficient  for  che  calibration  of  the  constitutive  model  under  consideration. 

2 .  Calibration  of  the  model  using  the  data  developed  in  Step  1. 

3.  Design  and  execution  of  the  series  of  verification  tests  which  provide  data  suffi¬ 
cient  for  defining  key  complex  material  response  features  that  are  to  be  modeled. 

A.  Direct  comparison  of  model  predicted  response  with  experimental  measurements 
through  the  use  of  a  constitutive  driver. 

The  constitutive  models  were  evaluated,  the  Fracture  F.nergy  Based  Model  (FEBM)  and  the  Endo- 
chronic  Concrete  Plasticity  Model  (tCPM) . 

While  there  are  very  many  constitutive  models  for  concrete  currently  available  in  the 
literature,  it  was  not  possible  within  the  scope  of  this  research  project  to  evaluate  all  of 
them,  although  the  methodology  presented  should  be  equally  applicable  to  all.  The  selection 
of  the  FEBM  and  the  ECPM  is  not  intended  to  endorse  these  models  as  the  better  ones.  The 
results  show  that,  although  they  are  able  to  predict  qualitatively  some  key  response  fea¬ 
tures  observed  in  the  verification  tests,  they  fail  to  predict  accurately  other  response 
features.  These  models  were  selected  because  they  are  two  of  the  more  recent  and  comprehen¬ 
sive  ones  and  also  the  theoretical  development  of  the  two  is  significantly  different. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-Si  units  of  measurement  can  be  converted  to  SI  (metric)  units  as  follows: 
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pounds  (force)  per  square  inch 


By 
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megapascals 
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megapascals 
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Evaluation  of  Non-linear  Constitutive 
Properties  of  Concrete 
Chapter  1 :  Introduction 


I .  I  Background . 

The  complex  stress-strain  response  features  of  concrete  have 

been  studied  and  discussed  considerably  in  the  technical  literature 

with  summaries  of  these  results  presented  in  many  textbooks  (e.g., 

(n*  (2) 

Neville  ,  Chen  ,  etc.).  Some  of  the  key  response  features  which 
can  have  a  significant  impact  on  the  development  of  constitutive  mod¬ 
els  include: 

a.  The  dependence  of  yield  or  maximum  strength  surfaces  on  con¬ 
fining  stress. 

b.  The  occurrence  of  plastic  volume  change  and  in  particular 
plastic  volume  change  due  to  pure  shear  loading. 

c.  Strain  softening  under  displacement  controlled  loading 
conditions . 

The  general  effects  of  these  features  on  stress  strain  response  are 
presented  in  Figures  1.1  through  1.3.  Figure  1.1  shows  the  effects 
of  increasing  confining  stress  on  the  stiffness  and  strength  of  a 
concrete  specimen  under  increasing  axial  stress.  Figure  1.2  shows 
the  increase  in  plastic  volume  change  (compaction)  due  to  pure  shear 
loading,  while  Figure  1.3  presents  test  results  indicating  strain 
softening  at  low  confining  stress  levels.  From  a  qualitative  stand¬ 
point,  these  response  features  can  be  partially  explained  by  consid¬ 
ering  the  internal  structure  of  concrete  at  the  three-phase  level 
(aggregate,  cement  matrix,  air  voids).  Such  a  description  is  pre¬ 
sented  graphically  in  Figure  1.4a,  which  represents  a  specimen  from  a 
well  proportioned  concrete  mixture  which  has  been  properly  consoli¬ 
dated  and  cured.  Furthermore,  it  is  assumed  that  the  nominal  maximum 
size  of  the  aggregate  is  small  compared  to  the  minimum  dimension  of 

*Numbers  in  parenthesis  indicate  references  presented  at  the  end  of 
the  dissertation. 
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the  specimen  so  that  on  any  plane  cut  through  the  specimen,  the  nor¬ 
mal  and  shear  stress  distributions  can  be  effectively  represented  as 
a  constant  average  value  (made  up  of  the  contributions  from  the 
aggregates,  cement  paste,  and  voids).  In  this  unloaded  specimen, 
some  cracks  will  be  observed,  primarily  at  the  aggregate-cement 
matrix  interface.  These  cracks  develop  during  the  curing  process  and 
are  primarily  due  to  differences  in  aggregate  and  cement  paste  stiff¬ 
nesses,  shrinkage,  and  thermal  properties.  When  compressive  axial 
stresses  alone  are  applied  in  the  vertical  direction  (Figure  1.4b) 
and  monotonically  increased,  microcracks  will  begin  to  propagate 
through  the  cement  matrix  which  results  in  a  net  decrease  in  the 
stiffness  of  the  specimen.  These  cracks  will  be  oriented  primarily 
vertically  on  the  exposed  outer  surface  due  to  the  tensile  strains 
developed  in  the  circumferential  direction.  Cracks  will  coalesce 
into  longer  vertical  cracks  as  the  stress  is  increased  up  to  some 
maximum  value  where  failure  will  occur.  Actually,  if  the  load  is 
applied  in  displacement  control,  the  load  will,  from  this  maximum 
value,  decrease  as  the  axial  strain  continues  to  increase,  however, 
significant  specimen  cracking  and  damage,  will  be  very  apparent).  On 
the  other  hand,  if  hydrostatic  stresses  are  applied  to  the  test  spec¬ 
imen  (Figure  1.4c),  the  initial  cracks  (Figure  1.4a)  will  begin  to 
close  so  that  shear  stresses  can  be  transferred  more  effectively 
across  the  aggregate-cement  paste  interfaces.  As  the  hydrostatic 
stress  increases,  the  normal  component  of  stress  on  a  typical  aggre¬ 
gate  surface  element  will  increase,  and  therefore  the  frictional 
shear  strength  at  this  point  will  increase.  In  general,  this  confin¬ 
ing  stress  effect  is  to  increase  the  effective  moduli  of  the  specimen 
as  well  as  the  maximum  strength  and  yield  surfaces. 

As  the  hydrostatic  stresses  increase,  the  volume  decreases  (com¬ 
paction)  ,  and  part  of  this  volume  change  is  irrecoverable  (plastic) 
as  shown  in  Figure  1.2.  This,  of  course,  is  unlike  the  response  of 
most  metals  where  essentially  no  plastic  volume  change  occurs  as  all 
plastic  strains  are  associated  with  shear  due  to  slip  between  grain 
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boundaries.  This  difference  in  response  between  polycrystalline 
metals  and  concrete  is  primarily  due  to  the  more  heterogeneous 
internal  structure  of  concrete.  As  voids  and  microcracks  close  under 
increasing  hydrostatic  stress,  the  volume  decreases.  However,  the 
volume  at  a  particular  level  of  hydrostatic  stress  is  not  necessarily 
a  minimum  for  that  level  of  stress  due  to  the  fact  that  microcracks 
and  voids  can  still  exist  due  to  the  complicated  microstress  distri¬ 
bution  in  the  specimen  which  can  cause  bridging  around  local  discon¬ 
tinuities-  As  a  pure  shear  stress  is  applied  (at  constant 
hydrostatic  stress)  ,  a  further  decrease  in  volume  occurs  (Fig¬ 
ure  1.2).  This  feature  can  be  explained  by  approximating  concrete, 
at  a  particular  level  of  hydrostatic  stress,  as  a  loosely  packed  sys¬ 
tem  of  incompressible  spheres.  If  any  pure  shear  distortion  is 
applied  to  the  system,  the  volume  must  decrease.  Also,  one  should 
consider  the  shearing  effect  of  concrete  along  a  rough  crack  where 
each  application  of  shear  tends  to  smooth  the  crack  and  therefore 
decrease  the  volume  of  the  specimen  containing  the  crack. 

The  third  response  feature  listed  above  (strain  softening)  has 
received  considerable  discussion  in  the  past  and  continues  to  be  vig¬ 
orously  debated.  The  main  issue  is:  "Should  strain  softening  be 
considered  a  property  of  the  material  or  a  result  of  the  structural 
geometry  and  loading  conditions  or.  the  test  specimen?"  The  later 
argument  is  generally  based  on:  (a)  the  frequent  observation  of  bar¬ 
reling  of  test  specimens  which  indicates  that  unwanted  shear  stresses 
are  being  applied  at  the  specimen  boundaries  and  that  the  stress  and 
strain  state  in  the  specimen  are  not  homogeneous;  (b)  significant 
cracking  of  the  test  specimen,  which  again  implies  an  inhomogeneous 
stress  state  in  the  material;  (c)  local  inhomogeneities  are  observed 
within  the  test  specimen  and  are  usually  associated  with  shear  band¬ 
ing,  etc.  The  central  question  concerning  softening  is  localization, 
which  implies  that  the  assumption  of  a  homogeneous  continuum  is  not 
valid.  If,  on  the  other  hand,  test  results  are  interpreted  so  that 
strain  softening  is  to  be  considered  a  material  property:  (d)  the 


applied  stresses  at  the  specimen  boundaries  should  be  uniform,  and 
the  current  specimen  geometry  must  be  given  as  a  constant  times  the 
initial  geometry;  (e)  the  internal  damage  or  response  mechanisms 
(i.e.,  fracture,  void  closure,  etc.)  should  be  uniformly  distributed 
over  the  specimen  volume  and  not  localized  in  bands  of  specific 
regions.  Figure  1 .4d  presents  conceptually  what  might  occur  in  a 
test  specimen  during  a  softening  test.  The  main  point  here  is  that 
fractures  in  aggregates,  void  closures  etc.  are  uniformly  distributed 
throughout  the  specimen  volume. 

The  complex  response  features  discussed  here  present  special 
problems  for  constitutive  models.  The  pressure  sensitivity  of  the 
yield  surface  simply  means  that  yield  and  maximum  strength  surfaces 
must  be  defined  in  terms  of  the  confining  stress  level.  The  plastic 
volume  change  problem  is  usually  addressed,  wi*:h  varying  degrees  of 
success,  by  defining  yield  surfaces  which  close  on  the  hydrostatic 
axis  or  by  the  use  of  rap  models.  When  the  formei  method,  to  account 
for  volume  change,  is  used,  constitutive  models  tend  to  overpredict 
lateral  strains  on  manv  stress  paths.  To  correct  this  problem,  a 
nonassoc iated  flow  rule  can  be  used  which  in  effect  changes  the 
direction  of  the  plastic  strain  vector  and  can  therefore  reduce  the 
lateral  strain  component. 

When  strain  softening  is  modeled  as  a  material  property,  the 
approach  used  generally  is  to  degrade  or  damage  the  maximum  strength 
surface  according  to  some  rule  which  relates  decrease  in  strength 
(or,  sav,  cohesion)  to  a  softening  parameter  (e.g.  plastic  strains) 
as  total  strains  continue  to  increase  under  decreasing  load. 

Once  constitutive  models  are  developed  which  reasonably  predict 
these  complex  response  features,  other  problems  are  encountered  when 
these  models  are  implemented  and  used  in  dynamic  finite-element 
codes.  The  nonassoc iated  flow  rule  results  in  a  nonsymmetric  mate¬ 
rial  stiffness  matrix  which  can  cause  a  significant  increase  in  com¬ 
putational  time.  Also,  realistic  load  paths  can  be  defined  along 
which  the  nonassoc iated  model  will  become  unstable  according  to 
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c’rtain  material  stability  postulates.  Strain  softening  models  also 
violate  material  stability  postulates  which  can  result  in  problems  of 
numerical  instability  and  nonunique  solutions.  Generally,  softening 
and  non-associated  flow  models  tend  to  be  undesirably  sensitive  to 
small  variations  in  prescribed  initial  conditions  for  dynamics 
problems . 

If  a  particular  response  feature  is  to  be  effectively  modeled, 
then  the  basis  for  the  development  of  the  model  should  be  sound 
repeatable  material  properties  test  data  completely  def  .ling  that 
response  feature.  On  the  other  hand,  material  stability  postulates 
should  be  considered  when  appropriate  but  should  not  be  the  driving 
force  in  the  development  of  the  model.  More  specifically,  the  com¬ 
plex  material  response  features  of  concrete  to  be  used  in  a  particu¬ 
lar  structure  should  be  evaluated  through  a  carefully  planned 
laboratory  material  properties  test  program,  designed  to  subject  the 
material  to  stress  and  strain  histories  similar  to  those  whi-'h  will 
occur  at  critical  regions  in  the  structure  under  the  design  loads. 
Furthermore,  generic  tests  (e.g.,  uniaxial  strain,  unconfined  com¬ 
pression,  etc.)  should  ha  conducted  to  determine  basic  material 
response  characteristics,  validity  of  homogeneity  and  isotropy 
assumptions,  and  values  for  parameters  used  to  calibrate  constitutive 
models  of  interest.  Finally,  the  concrete  material  should  be  sub¬ 
jected  to  complex  load  path  tests,  which  can  be  used  to  evaluate  the 
consistency  and  predictive  capability  of  potential  concrete  constitu¬ 
tive  models. 

1 . 2  Ob j  ectives 

The  objectives  of  this  research  are: 

a.  The  development  of  a  methodology  for  evaluating  constitutive 
models  for  plain  concrete. 

b.  The  application  of  this  methodology  in  the  evaluation  of  two 
advanced  :onstitutive  models. 
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1 . 3  Scope 

The  method  of  evaluation,  as  developed  herein,  consists  of  the 
following  steps: 

1 .  The  design  and  execution  of  a  series  of  material  properties 
tests  which  provide  data  sufficient  for  the  calibration  of 
the  constitutive  model  under  consideration. 

2.  Calibration  of  the  model  using  the  data  developed  in  step  1. 

3.  Design  and  execution  of  a  series  of  verification  tests  which 
provide  data  sufficient  for  defining  key  complex  material 
response  features  that  are  to  be  modeled. 

4.  Direct  comparison  of  model  predicted  response  with  experi¬ 
mental  measurements  through  the  use  of  a  constitutive 
driver . 

The  two  constitutive  models  to  be  evaluated  are  the  Fracture  Energy 
Based  Model  (FEBM) (3)  and  the  Endochronic  Concrete  Plasticity  Model 
( ECPM) (4) . 

While  there  are  very  many  constitutive  models  for  concrete,  cur¬ 
rently  available  in  the  literature,  it  was  not  possible  within  the 
scope  of  this  research  project  to  evaluate  all  of  them,  although  the 
methodology  presented  should  be  equally  applicable  to  all.  The 
selection  of  the  FEBM  and  the  ECPM  is  not  intended  to  endorse  these 
models  as  the  better  ones.  The  results  show  that  although  they  are 
able  to  predict  qualitatively  some  key  response  features  observed  in 
the  verification  tests  they  fail  to  predict  accurately  other  response 
features.  These  models  were  selected  because  they  are  two  of  the 
more  recent  and  comprehensive  ones  and  also  the  theoretical  develop¬ 
ment  of  the  two  is  significantly  different. 

1 .4  Outline 

In  Chapter  2,  general  stress-strain  response  features  of  con¬ 
crete  will  be  discussed  along  with  some  of  the  implications  of  strain 
softening.  Loaders  and  test  devices  that  were  used  in  generating 
test  results  for  this  research  project  along  with  instrumentation 
used  are  then  discussed.  Finally  calibration  and  verification  test 
results  are  presented  and  discussed.  Basic  concepts  of  elasticity 
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and  plasticity  are  presented  in  Chapter  3  with  emphasis  on  generaliz¬ 
ing  the  results  observed  in  simple  material-properties  tests.  Also, 
in  Chapter  3,  a  detailed  derivation  of  the  equations  for  the  FEBM  and 
the  ECPM  are  presented  and  the  calibration  of  model  parameters  is 
discussed.  In  Chapter  4,  the  two  models  are  calibrated,  exercised 
against  the  verification  tests  results,  and  comparisons  of  tests 
results  versus  model  predictions  are  made.  Conclusions  and  recom¬ 
mendations  for  future  research  are  presented  in  Chapter  5. 


Axial  strain 

F igure 1 .  1 . Ef f ects  of  confining  stress  on  maximum  strength 
and  stiffness. 
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Figure  1.2.  Decrease  in  volume (compaction) 
under  hydrostatic  and  pure 
shear  loading. 


(a)  Unloaded  test  specimen. 


(b)Cracks  observed  in  unconfined 
compression  test. 


Figure  1 .4  .  Internal  characteristics  of  concrete (sheet  1  of  2) 
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(c)Crack  and  void  closure  due 
to  nydrostatic  stresses. 


(d)Aggregate  fracture, softening 
concept . 


Figure  1. h. Internal  characteristics  of  concrete(sheet  2  of  2) 


Chapter  2 

Stress  Strain  Response  of  Concrete 


2.1  General 

In  Chapter  1,  a  general  discussion  was  presented  concerning  prob¬ 
lems  encountered  in  large  scale  structural  analysis  of  critical  struc¬ 
tures  and  questions  many  researchers  have  regarding  various  stress 
strain  response  phenomena  of  concrete.  Considerable  experimental 
data  exist  for  concrete  subjected  to  one  dimensional  and  axisymmetric 
load.  Also  several  studies  have  been  directed  at  strain  softening, 
biaxial  stress  loading  and  fully  three  dimensional  stress  loading. 
Results  from  many  of  these  test  programs  axe  summarized  and  discussed 
by  Hegemier,  et  al.  (5).  Green  and  Swanson  (6)  conducted  an 
extensive  study  of  the  general  constitutive  properties  of  concrete  at 
intermediate  pressure  levels  (i.e.  10  -  12  ksi)  while  Van  Mier  (7) 
addressed  strain  softening  under  multiaxial  loading  conditions. 

Gerstle  et  al.  (8)  through  a  cooperative  research  effort  showed  the 
sensitivity  of  tests  results  to  test  devices  (boundary  conditions) 
and  test  procedures.  The  vast  majority  of  the  existing  test  data  is 
based  on  simple  or  proportional  load  paths  (I.e.  the  ratio  of  the 
applied  stresses  remains  constant).  In  this  Chapter,  specific 
stress-strain  response  characteristics  of  concrete  for  simple  and 
complex  load  paths  will  be  discussed.  Furthermore,  an  attempt  is 
made  to  show  that  high  quality,  repeatable,  and  consistent  test  data 
can  be  obtained  for  well  prepared  concrete  specimen  and  furthermore 
these  data  can  and  should  be  used  In  calibrating  constitutive  models. 
Nominal  unconfined  compressive  strengths  for  concretes  tested  in  this 
study  range  from  f^  =  2  ksi  to  f^  =  7  ksi.  It  is  not  intended  to 
present  herein  a  broad  discussion  of  the  many  different  complex 
response  characteristics  of  concrete,  but  rather  to  discuss  those 
features  which  are  of  primary  Importance  in  constitutive  modeling  in 
general  and  specifically  for  the  fracture  energy  based  model  and  the 
endochronic  model.  Essentially  two  types  of  tests  will  be  discussed, 
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calibration  tests  and  verification  tests.  Calibration  tests  are 
conducted  to  determine  values  for  parameters  or  coefficients  used  in 
the  mathematical  formulation  of  a  particular  model.  Verification 
tests  are  designed  to  evaluate  specific  predictive  capabilities  of  a 
constitutive  model.  Obviously  a  calibration  test  for  one  model  might 
serve  as  a  verification  test  for  another  model. 

Before  discussing  test  results  it  is  important  to  keep  in  mind 
that  in  a  test  we  are  measuring  quantities  such  as  force,  (F) ,  pres¬ 
sure  (P) ,  and  displacement  (U) ,  then  dividing  these  (i.e.  force  and 
displacement)  by  initial  areas  Aq  and  lengths  Lq  to  obtain  stress  and 
strain.  The  equations  we  use  are  of  course 

STRESS  =  F  /  A  (1) 

o 

STRESS  =  U  /  L 

o 

The  validity  of  using  these  equations  to  interpret  or  infer  material 
constitutive  properties  should  be  determined  based  on  the  following 
conditions  (Parlseau  (9)): 

1.  Test  specimen  must  be  homogeneous 

2.  A  homogeneous  state  of  stress  must  exist  in  the  specimen  at 
all  times. 

3.  No  significant  changes  in  specimen  geometry  can  occur  during 
a  test. 

When  these  conditions  are  met  one  can  reasonably  assume  constitutive 
properties  derived  from  equation  (1)  are  real  material  properties. 
Furthermore,  if  an  appropriate  number  of  tests  are  conducted  along 
load  paths  which  capture  the  key  response  features  of  the  material 
one  can  construct  rational  constitutive  models  which  effectively 
represent  the  material  response  under  a  wide  range  of  load  histories. 
If  on  the  other  hand  these  conditions  are  not  met,  one  must  be  care¬ 
ful  in  inferring  material  properties  from  equation  (1)  and  in  using 
constitutive  models  constructed  from  such  tests.  The  degree  to  which 
the  above  conditions  are  met  has  special  significance  in  interpreting 
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strain  softening  as  a  material  response  phenomenon  or  a  structural 
response  feature.  When  concrete  is  modeled  as  a  homogeneous,  isotro¬ 
pic  material  the  assumption  is  made  that  the  aggregates  and  pores  are 
randomly  distributed  throughout  the  cement  paste  matrix  and  that  the 
microcrack  surface  and  width  is  relatively  small  when  considering  a 
representative  volume  of  the  material.  This  representative  volume  in 
general  must  be  large  enough  so  that  the  effects  of  stress  concen¬ 
trations,  and  any  other  discontinuities  in  the  material  can  be 
smeared  over  the  volume  so  that  an  equivalent  homogeneous  stress  and 
strain  state  can  be  determined  which  effectively  characterizes  all 
the  important  features  of  material  response.  The  reason  for  con¬ 
ducting  material  properties  tests  is  precisely  to  measure  these 
stress  and  strain  states  along  load  paths  which  are  critical  consid¬ 
ering  the  application.  At  high  stresses  and  especially  in  the  soft¬ 
ening  region  the  validity  of  assumptions  of  homogeneity  become  more 
and  more  questionable.  Many  researchers  have  pointed  out  that  at 
ultimate  strength  and  in  the  softening  region  there  are  localizations 
of  damage  (eg.  shear  banding)  and  that  test  specimens  exhibiting 
these  effects  must  be  considered  structural  elements  and  not  as  mate¬ 
rials  subjected  to  states  of  homogeneous  stress  and  strain. 

Bazant  (10)  found  that  strain  softening  can  be  observed  only  when 
local  inhomogeneities  exist  in  a  material.  The  question  is,  how 
significant  are  the  inhomogeneities  in  terms  of  the  effects  they  have 
on  stress-strain  response  features  of  interest.  Once  this  question 
is  answered  one  can  determine  what  conditions  (e.g.  inhomogeneities) 
must  be  considered  in  constructing  a  rational  constitutive  model  for 
a  particular  problem.  Local  inhomogeneities  lead  to  strain  softening 
which  can  be  identified  as  a  material  instability  according  to 
Drucker's  postulates  (see  Chapter  3).  Drucker  (17)  points  out 
"Philosophically,  perhaps,  all  macroscopic  systems  are  stable  in  the 
sense  that  if  all  conditions  are  or  could  be  taken  into  account  the 
complete  behavior  of  the  system  could  be  followed  in  detail... 
Instability  as  normally  understood  may  arise  when  some  but  not  all  of 


the  attributes  of  a  system  are  considered."  For  strain  softening  in 
concrete  the  instability  is  due  to  the  assumptions  that  the  stress 
strain  response  is  rate  independent  and  homogeneous. 

In  this  chapter,  loaders,  test  devices,  and  instrumentation  com¬ 
monly  used  in  conducting  material  property  tests  will  be  discussed 


(Section  2.2).  The  location  of  the  loaders  and  test  devices  and  the 
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2.2  Loaders,  Test  Devices,  and  Instrumentation 

2.2.1  100-kip  Servo-Controlled  Hydraulic  Loader  (WES) 

The  100-kip  servo-hydraulic  loader  is  capable  of  applying  ten¬ 
sile  or  compressive  loads  up  to  a  maximum  level  of  100,000  lb.  The 
loader  can  be  manually  or  computer  controlled  to  produce  a  desired 
load  or  displacement  rate.  The  input  to  the  servo-control  unit  is 
produced  by  an  arbitrary  digital  function  generator  which  can  be  pro¬ 
grammed  to  produce  an  infinite  number  of  load  or  displacement  his¬ 
tories.  The  applied  force  is  measured  by  a  load  cell  while  the 
displacement  of  the  loader  head  is  measured  by  an  internal  linear 
variable  displacement  transducer  (LVDT) .  Strains  are  usually  mea¬ 
sured  by  epoxy-backed  constantan  resistance  strain  gages  bonded  to 
the  specimen  with  high-strength  strain-gage  adhesive.  Unconfined 
compression  tests,  tension  tests,  and  some  low  confinement  triaxial 
tests  can  be  conducted  in  this  loader. 

2.2.2  440-kip  Universal  Testing  Machine  (WES) 

The  440-kip  loader  can  apply  tensile  or  compressive  loads  up  to 
a  maximum  of  440,000  lb.  The  machine  is  displacement  controlled  by 
manually  adjusting  hydraulic  valves  to  achieve  the  desired  displace¬ 
ment  rate.  The  axial  load  is  measured  by  an  internal  load  cell  while 
displacements  are  measured  by  an  external  LVDT  across  the  loading 
heads.  Strains  are  usually  measured  by  strain  gages  bonded  to  the 
specimen.  This  machine  is  used  to  conduct  unconfined  compression 
tests  and  cylindrical  triaxial  compression  tests. 

2.2.3  2400-kip  Loader  (WES) 

The  2400-kip  servo-hydraulic  loader  can  apply  tensile  or  com¬ 
pressive  loads  up  to  a  maximum  of  2,400,000  lb.  The  load  or  dis¬ 
placement  rate  are  controlled  by  an  analog  function  generator  capable 
of  ramp,  sine,  or  triangular  functions.  The  load  is  measured  by  an 
internal  load  cell  while  the  displacements  are  measured  by  an 
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external  LVDT.  This  device  is  used  to  conduct  unconfined  compression 
tests,  cylindrical  triaxial  tests,  and  multiaxial  compression  tests. 

2.2.4  40-ksi  Cylindrical  Triaxial  Chamber  (WES) 

The  40-ksi  cylindrical  triaxial  test  device  is  capable  of  con¬ 
ducting  multiaxial  tests  on  nominal  2.125-inch  diameter  by  3.5-inch 
long  cylindrical  specimens  at  confining  stresses  up  to  40,000  psi  and 
axial  stresses  up  to  125,000  psi.  Figure  2.1(a)  shows  a  cross- 
section  of  the  device.  The  confining  pressure  is  developed  by  a 
manually  controlled,  air-driven  hydraulic  pump  rated  at  75,000  psi. 

The  fluid  used  in  the  device  is  a  low-viscosity  white  mineral  oil. 

Eight  7/8-inch  diameter  bolts  clamp  the  upper  and  lower  platens 
together  and  carry  the  loads  produced  by  the  fluid  pressure  against 
the  upper  and  lower  cell  caps.  The  device  has  the  capability  to  exit 
up  to  four  channels  of  instrumentation  through  Fusite-type  fittings 
in  the  base. 

In  preparing  for  a  test,  hardened  steel  end  caps  were  placed  on 
each  end  of  a  specimen.  The  specimen  was  subsequently  encased  in  an 
impermeable  60  durometer  neoprene  membrane.  Hose  clamps  were  used  to 
secure  the  membrane  at  each  end  of  the  specimen.  The  jacketed  speci¬ 
men  was  placed  in  the  device,  and  the  device  was  completely  assem¬ 
bled.  This  assembly  was  then  moved  as  a  unit  and  placed  in  the 
440-kip  universal  testing  machine.  A  swivel  was  placed  on  the  top  of 
the  ram  to  compensate  for  any  eccentricities  in  the  device  or  the 
testing  machine. 

2.2.5  30-ksi  Multiaxial  Test  Device  (WES) 

The  30-ksi  multiaxial  test  device  is  used  in  conjunction  with 
the  2400-kip  universal  testing  machine  to  test  6-inch  by  6-inch  by 
36-inch  rectangular  prismatic  specimens  under  a  wide  range  of  three- 
dimensional  compressive  stress  states.  A  vertical  cross-section  of 
the  device  is  shown  in  Figure  2.1(b).  The  test  device  incorporates 
fluid  cushion  technology  to  minimize  surface  shear  friction  and 
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(a)  40  ksi  cylindrical  triaxial  chamber 
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(b)  30  ksi  multiaxial  test  device 
Figure  2.1  WES  test  devices 
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nonuniform  pressure  distributions  normally  encountered  in  other  types 
of  test  devices.  The  axial  stress  is  developed  by  the  2400-kip  uni¬ 
versal  testing  machine  and  transferred  through  high-strength  steel 
bearing  blocks.  Strains  in  the  specimen  are  measured  in  the  central 
third  of  the  specimen.  This  location  takes  advantage  of  the  6: 1 
aspect  ratio  of  the  specimen,  thus  reducing  the  end  effects  induced 
by  the  rigid  bearing  blocks. 

The  test  fixture  is  essentially  a  pressure  vessel  which  is 
divided  into  four  internal  chambers  to  accommodate  the  individual 
bladder/seal  combinations  which  make  up  the  fluid  cushions.  The 
fluid  cushions  apply  stresses  of  up  to  30,000  psi  directly  to  the 
vertical  faces  of  the  specimen.  The  fluid  pressure  is  developed  by 
manually  controlled,  air-driven  hydraulic  pumps  rated  at  75,000  psi. 

A  low-viscosity  white  mineral  oil  is  used  as  the  pressurizing  fluid. 
Strains  are  measured  using  internal  embedded  integral  lead  strain 
gages. 

2.2.6  Eindhoven  Cubical  Cell  (UEN) 

The  Eindhoven  cubical  cell  was  developed  at  the  University  of 
Eindhoven,  Nether  1  ar.ds  [7].  The  test  device  consists  of  three 
identical  loading  frames  hanging  in  a  fourth  main  frame  structure. 

The  three  loading  frames  are  independently  suspended  by  steel  cables 
and  not  connected  to  each  other.  The  frames  are  fixed  vertically 
with  the  two  frames  in  the  horizontal  direction  free  to  move.  This 
causes  some  lack  of  symmetry  in  the  loading  system,  but  this  effect 
was  found  to  be  negligible.  Loads  are  developed  by  three  independent 
servo-controlled  hydraulic  actuators  with  nominal  maximum  capacities 
of  450  kips.  These  loads  are  transferred  to  the  test  specimen,  which 
is  normally  a  four  inch  cube,  through  brush  bearing  platens.  The 
brush  bearing  platens  are  designed  to  reduce  the  shear  restraint  at 
the  concrete  surface.  Specimen  formations  (for  three-dimensional 
tests)  are  inferred  from  LVDT  measurements  of  relative  displacement 
between  the  steel  blocks  upon  which  the  brush  rods  are  clamped.  A 
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schematic  of  one  load  frame ,  actuator  rystem  is  shown  in 
Figure  2.2(a) . 

2.2.7  Colorado  Cubical  Cell  (UCB) 

The  Colorado  cubical  cell  was  developed  at  the  University  of 
Colorado  [11]  for  testing  materials  under  multiaxial  compressive 
loads.  Pressures  are  developed  by  hand  pumps  and  applied  to  the 
specimen  through  polyurethane  membranes  filled  with  a  silicone  fluid. 
Deformations  of  the  specimen  are  measured  with  a  proximeter  probe 
system.  The  cell  is  designed  to  develop  a  maximum  stress  of  12  to 
15  ksi  on  4-inch  cuoical  test  specimens.  Key  features  of  the  cubical 
cell  are  shown  in  Figure  2.2(b).  It  is  important  to  note  that  the 
Colorado  cubical  cell  is  completely  stress  or  load  controlled. 


20 


(a)  Eindhoven  cubical  cell  (from  (18)) 


(b)  Colorado  cubical  cell 

Figure  2.2  Eindhoven  and  Colorado  test 
devices 
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2.3  Unconfined  Compression  Tests 

To  the  naked  eye  concrete  appears  mainly  as  a  two  phase  compos¬ 
ite  material  consisting  of  different  size  aggregates  embedded  in  a 
cement  paste  matrix.  One  can  usually  observe  some  small  pores  also, 
but  these  can  be  kept  to  a  minimum  if  proper  care  is  taken  in  mixture 
proportioning,  specimen  casting,  and  consolidation.  Through  a  micro¬ 
scope  one  can  observe  microcracks  primarily  located  at  the  coarse 
aggregate-mortar  interface.  As  the  magnification  of  the  microscope 
increases  one  observes  that  the  number  of  phases  in  the  concrete  com¬ 
posite  increases.  For  the  purposes  of  this  study  the  two  phase 
approximation  is  useful  in  explaining  the  stress  strain  response  fea¬ 
tures  of  concrete.  However,  the  constitutive  models  addressed  herein 
assume  homogeneous  isotropic  materials.  The  appropriateness  of  these 
assumptions  depends  on  defects  (eg.  microcrack,  pores,  etc.)  being 
spatially  distributed  in  a  random  manner  and  relatively  small  in  num¬ 
ber  compared  to  the  specific  material  volume  of  interest  as  discussed 
in  Section  2.1. 

Figure  2.3  presents  the  results  of  an  unconfined  compression 
test  of  concrete  conducted  in  displacement  or  stroke  control.  The 
region  OA  might  be  considered  the  elastic  region,  with  point  A  defin¬ 
ing  the  proportional  limit,  while  AB  is  the  hardening  region  and  BC 
and  larger  strains  could  be  considered  the  softening  region.  It  is 
generally  observed  that  the  microcracks  which  exist  at  zero  load,  at 
aggregate-mortar  interfaces,  increase  in  surface  area  and  number  as 
the  stress  is  increased  from  point  0  to  A.  At  some  point  near  A 
microcrac.ks  begin  to  propagate  through  the  mortar  and  continue  to 
propagate  and  coalesce  in  a  stable  manner  up  to  some  point  below  B  in 
the  region  AB.  At  this  point  the  cracks  begin  to  propagate  in  an 
unstable  manner  and  failure  will  be  observed  at  about  point  B  in  a 
load  control  test  device.  Considerable  discussion  is  given  to  the 
points  described  above  by  Kotsovos  and  Newman  [121.  If  at  any  point 
in  the  program  of  loading  for  Figure  2.3  unloading  occurs  followed  by 
reloading  the  load  will  return  to  essentially  the  same  point  where 
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Figure  2.3  Stress-strain  response,  for  unconfined  compressive 

test  (stroke  control! 

unloading  occurred  and  then  the  stress  strain  response  as  shown  in 
Figure  2.3  will  continue  under  monotonic  loading.  Unloading  reload¬ 
ing  moduli  will  be  somewhat  less  than  the  initial  value  and  some 
hysteresis  will  be  observed.  At  higher  and  higher  strains  the  load¬ 
ing  unloading  moduli  become  softer  and  softer  with  more  and  more 
hysteresis.  At  strains  greater  than  those  associated  with  point  B  in 
Figure  2.3  visual  cracks  will  be  detected  on  the  outside  surface  of 
the  specimen.  These  cracks  will  essentially  be  vertical  and  defi¬ 
nitely  indicate  that  the  specimen  is  no  longer  an  intact  continuum. 
This  observation  provides  the  basis  for  most  arguments  that  softening 
is  a  structural  phenomenon  as  opposed  to  a  material  property.  How¬ 
ever,  it  should  be  kept  in  mind  that  these  observations  are  for  the 
unconfined  test. 

Some  other  quantitative  observations  can  be  made  concerning  the 
material  response  features  for  the  f^  =  6.5  ksi  concrete  as  shown  in 
Figure  2.3.  For  this  concrete  specimen  the  diameter  is  2.065  in.  and 
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length  is  4.25  in.  The  coarse  aggregate  is  limestone  with  maximum 
size  of  3/8  inch,  the  water-cement  ratio  for  the  mixture  is  0.46,  and 
the  test  specimen  is  approximately  one  year  old.  Stress  is  force  per 
unit  original  area  and  strains  are  measured  by  3/4  and  1-inch  foil 
strain  gages  glued  to  the  specimen  and  by  an  LVDT  measuring  loader 
head  movement  corrected  for  the  compliance  of  the  loader  system. 

Figure  2.3  represents  a  continuous  loading  of  a  specimen  to  very 
large  strains.  In  order  to  further  explore  the  response  character¬ 
istics  of  concrete,  especially  in  the  linear  elastic  region  a  series 
of  tests  were  conducted  to  assess  the  linear  elastic  features  of  the 
material  at  stresses  below  the  "yield  poinc."  For  a  companion 
specimen  to  the  one  used  in  the  ’■est  of  Figure  2.3,  a  series  of  load¬ 
ing,  unloading,  and  reloading  tests  were  conducted  to  determine 
accumulations  of  plastic  strains,  softening  of  unloading-reloading 
moduli  and  indications  of  damage.  The  test  procedure  was  to  load  up 
to  a  specified  level  of  axial  stress  then  unload  to  zero  stress.  The 
test  specimen  was  then  removed  from  the  loader  tested  dynamically  to 
determine  its  compressional  wave  velocity,  (which  is  a  measure  of 
dynamic  modulus) ,  measured  with  a  micrometer  to  determine  permanent 
strains,  and  then  placed  back  in  the  loader  and  the  test  procedure 
repeated  to  a  higher  level  of  axial  stress.  Since  the  ultimate 
strength  of  the  concrete  of  Figure  2.3  is  f^  =  7.5  ksi  (which  is  the 
result  of  one  year  of  aging  on  the  nominal  6.5  ksi  concrete),  these 

load-unload  tests  were  conducted  at  0.3  f ' ,  0.5  f *  ,  0.75  f '  ,  0.9  f', 

c  c  c  c 

and  on  the  final  loading  program  the  specimen  was  loaded  to  an  axial 
strain  of  e  =  1.5%.  The  results  of  this  test  program  are  presented 
In  Figure  ?,4.  The  relatively  smali  effect  of  load-unload  cycles  on 
E  and  wave  velocity  Indicates  the  linear  elastic  response  character¬ 
istics  of  the  f'  =  6.5  ksi  concrete  for  stresses  below  A.  Figure  2.5 
c 

presents  a  composite  plot  of  the  f^  -  6.5  ksi  and  the  f^  =  2  ksi 
stress  strain  curves,  which  represent  upper  and  lower  bounds  of  con¬ 
crete  strengths  studied  in  this  report. 


Mean  Norma  1  Strnss 
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2.4  Hydrostatic  Compression  Test 

The  results  of  a  typical  hydrostatic  compression  test  are  pre¬ 
sented  in  Figure  2.6  with  unload-reload  cycles.  The  bulk  modulus  is 
calculated  directly  from  these  test  results  (see  Chapter  3)  as  well 
as  unloading  and  reloading  moduli.  The  hydrostatic  test  represents 
an  upper  bound  on  the  effects  confinement  can  nave  on  concrete 
response.  The  pure  hydrostatic  stress  state  tends  to  arrest  micro¬ 
crack  growth,  since  for  the  most  part  microcrack  growth  is  associated 
with  deviatoric  stress  components.  The  flat  portions  indicated  on 
Figure  2.6  at  high  stress  levels  are  due  primarily  to  creep,  which 
results  in  the  finite  amount  of  time  (here  approximately  20  sec)  in 
changing  the  load  from  increasing  to  decreasing.  This  feature  (rate 
dependence)  should  be  kept  in  mind  when  evaluating  test  results  and 
predictions  from  rate  independent  constitutive  models.  Also,  up  to  a 
mean  normal  stress  (MNS)  of  approximately  f^  the  response  is  linear 
with  very  little  plastic  volume  strains  observed  upon  unloading.  At 
MNS  =  10  ksi  the  bulk  modulus  begins  to  soften  and  approaches  a  near 
constant  value  at  approximately,  MNS  =  18  ksi. 


fi . 5-ks 1  Concrete 


Volumetric  Strain,  '/, 


Figure  2.6  Hydrostatic  test  with  unload-reload  cycles 


27 


2.5  Uniaxial  Strain  Test 

Figure  2.7  presents  results  in  terms  of  axial  stress  versus 
axial  strain  for  a  typical  uniaxial  strain  test  of  a  concrete  with 
average  nominal  strength  of  6,500  psi.  The  load  path  for  this  test 
in  the  principal  stress  space  is  also  presented  in  Figure  2.7.  The 
uniaxial  strain  test  can  be  very  important  in  calibrating  certain 
constitutive  models  because  of  its  simple  radial  displacement  bound¬ 
ary  conditions  (i.e.  e  radial  =  0  or  a  constant).  Like  the  hydro¬ 
static  test  the  uniaxial  strain  test  shows  the  significant  effects  of 
confining  stresses,  but  here  deviatoric  (shear)  stresses  exist  and 
there  is  more  of  a  tendency  for  microcrack  propagation.  In  Fig¬ 
ure  2.8  a  composite  plot  is  presented  of  unconfined  compression, 
hydrostatic  compression,  and  uniaxial  strain.  The  uniaxial  strain 
test  shows  a  stiff  response  up  to  about  an  axial  stress  of  f^_.  This 
is  seen  clearer  when  the  lower  end  of  the  test  is  expanded  (Fig¬ 
ure  2.8b)  and  the  unconfined  compression  test  results  are  compared 
with  the  uniaxial  strain,  and  hydrostatic  compression  tests.  One  can 
reasonably  assume  that  elastic  response  is  occurring  in  all  three 

tests  up  to  a  level  of  axial  stress  of  75%  to  80%  of  f'  for  this  con- 

c 

crete.  From  Figure  2.8b,  it  can  be  seen  that  all  three  tests  start 
out  with  equivalent  moduli  in  terms  of  ratios  of  axial  stress  to 
axial  strain  that  are  essentially  linear  and  continue  up  to  an  axial 
stress  of  approximately  5  ksi.  (Note  we  defined  yield  in  the  uncon¬ 
fined  compression  test  at  4.8  ksi).  At  stresses  above  this  value 
each  curve  begins  to  strain  harden,  with  the  hydrostatic  test  harden¬ 
ing  being  the  largest,  the  unconfined  test  indicating  apparent  soft¬ 
ening,  and  the  uniaxial  strain  test  in  between. 
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2.6  Triaxial  Compression  Test 

2.6.1  Conventional  Triaxial  Compression  Test 

A  Conventional  Triaxial  Compression  Test  (CTC)  is  defined  herein 
to  consist  of  two  loading  branches.  The  first  branch  consists  of  a 
pure  hydrostatic  loading  to  a  specified  level  of  hydrostatic  or  mean 
normal  stress  (MNS) .  The  second  branch  is  obtained  by  increasing  the 
axial  stress  while  holding  the  confining  stress  constant.  The  load 
(stress)  path  for  the  CTC  test  is  shown  in  Figure  2.9  in  the  Rendulic 
plane.*  It  is  important  to  note  that  path  AB 


Figure  2.9.  Conventional  triaxial  compression 
stress  path. 

for  the  CTC  test  includes  changes  in  the  hydrostatic  component  of 

stress  as  well  as  the  deviatoric  component  of  stress.  Typical  axial 

stress  versus  axial  strain  curves  for  the  fT  =6.5  ksi  concrete  are 

c 

shown  in  Figure  2.10  for  two  different  confining  stress  levels. 
Critical  information  obtained  from  these  tests  include  strength, 
ductility,  and  loading  and  unloading  moduli  as  a  function  of  con¬ 
fining  stress  level.  In  Figure  2.10  axial  stress  versus  axial  strain 


*  The  Rendulic  plane  is  the  plane  in  principal  stress  space  in  which 
two  of  the  principal  stresses  are  always  equal  (i.e.  it  is  the 
plane  which  contains  all  possible  axisymmetric  stress  states.) 
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are  presented  for  CTC  tests  on  two  similar  test  specimens  at  each 
confining  stress  level  to  demonstrate  the  repeatability  of  the  test 
results.  From  these  test  results  one  concludes  that  a  distinct  yield 
point  is  not  observed  yet  the  material  is  behaving  in  a  consistent 
strain  hardening  manner.  Furthermore,  CTC  test  results  on  similar 
specimen  at  a  higher  confining  stress  show  an  increase  in  strength 
and  loading  modulus  due  to  confining  stress,  are  also  repeatable,  and 
still  do  not  exhibit  a  distinct  yield  point.  The  effects  of  initial 
yield  surface  location  on  plastic  strain  predictions  will  be  evalu¬ 
ated  in  Chapter  4.  Test  results  presented  in  Figure  2.10  were 
carried  out  to  very  large  engineering  strains  (approximately  20%). 
Recalling  the  conditions  for  interpreting  material  properties  test 
results  (Section  2.1),  and  since  the  intent  of  this  study  is  to 
evaluate  constitutive  models  based  or  small  strain  theory,  the  axial 
stress  at  an  axial  strain  of  10%  will  be  taken  as  the  ultimate 
strength  foi  purposes  of  constructing  failure  envelopes.  It  should 
be  mentioned  at  this  point  that  failure  or  ultimate  strength  has  not 
yet  been  defined  for  concrete  subjected  to  confining  stress  condi¬ 
tions.  In  fact  for  any  of  the  confined  tests  presented  thus  far,  if 
at  any  point  (see  Figure  2.10)  in  the  loading  program  the  stress  is 
decreased,  along  the  same  path  followed  during  loading,  unloading 
will  occur  and  reloading  back  to  the  point,  approximately,  where 
unloading  started  can  be  achieved  with  little  if  any  evidence  of 
damage,  in  terms  of  measured  stresses  and  strains.  However  an  indi¬ 
cation  of  damage  which  does  occur  in  the  test  at  these  high  confining 
stresses  and  very  large  axial  strains,  is  the  decrease  in  compressive 
wave  velocity.  For  the  f^  =  6.5-ksi  concrete  the  compressive  wave 
velocity  of  an  untested  specimen  is  approximately  Cp  =  15,600  fps. 
After  testing  to  approximately  20%  axial  strain  at  15  ksi  confining 
stress  the  post  test  wave  velocity  ranged  (for  different  test  speci¬ 
men)  from  Cp  =  8,000  fps  to  Cp  =  10,800  fps.  Shear  wave  velocity, 
pretest,  was  measured  at  approximately  Cs  =  9,300  fps  while  post  test 
values  ranged  from  Cs  =  5,700  fps  to  CS  =  6,340  fps.  These 


measurements  indicate  similar  degradation  in  compressive  and  shear 
wave  velocities.  A  maximum  strength  surface  based  on  CTC  tests 
stress  levels  at  10%  axial  strain  is  presented  for  the  Rendulic  plane 
in  Figure  2.11.  Failure  stresses  at  zero  confinement  and  at  the  3 
ksi  and  4  ksi  levels  were  based  on  maximum  stress  attained  in  the 
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ligure  2.11  Maximum  strength  surface  for  6.5-ksi  concrete 
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2.6.2  Brittle  Ductile  Transition 

The  brittle  ductile  transition  (BDT)  is  defined  here  as  the  con¬ 
fining  stress  which  separates  the  apparent  softening  response  region 
and  the  strain  hardening  response  region  in  stress  space.  Lateral 
strains  in  low  confinement  tests  are  predominantly  tensile  in  the 
softening  region,  while  predominantly  compressive  lateral  strains 
occur  in  the  high  confinement  hardening  region.  This  would  seem  to 
indicate  that  the  transition  from  softening  to  hardening  response 
(i.e.  BDT)  could  be  associated  with  vanishing  lateral  strains  in  a 
CTC  test.  Based  on  this  assumption  the  uniaxial  strain  test  has  spe¬ 
cial  significance  in  determining  the  BDT,  since  this  test  is  based  on 
zero  lateral  strains.  The  uniaxial  strain  load  path  is  plotted  rela¬ 
tive  to  the  maximum  strength  surface,  as  defined  earlier,  in  Fig¬ 
ure  2.12.  It  is  interesting  to  note  that  the  intersection  of  the 
tangent  to  the  load  path  at  low  stress  and  the  higher  load  path 
tangent  is  very  close  to  the  3DT.  This  point  is  used  in  the  FEBM  for 
calibration  purposes.  Figure  2.13  presents  CTC  test  results  very 

near  the  BDT  for  the  f'  =  6.5-ksi  concrete,  which  will  be  taken  here 

c 

to  be  at  a  confining  stress  of  f^  =  4  ksi. 


Figure  2.12  Brittle  ductile  transition  relative  to 
uniax.  strain  test 


Figure  2.13  Determining  brittle  ductile  transition 

2.6.3  Deviatoric  Compression  Test 

A  deviatoric  compression  test  (DCT)  is  defined  here  to  consist 
of  two  loading  branches.  The  first  branch  consists  of  a  pure  hydro¬ 
static  loading  to  a  specified  level  of  mean  normal  stress.  The 
second  branch  is  a  pure  deviatoric  branch  along  a  stress  path,  which 
in  the  Rendulic  plane  can  be  defined  as  A =  -2Ao^  .  A  typical 
ideal  DCT  load  path  is  shown  in  Figure  2.14.  Figure  2.15  presents 
axial  stress  versus  axial  and  lateral  strains  for  DCT  conducted  on 

f'  =  6.5-ksi  concrete  at  two  different  mean  normal  stress  (KNS) 
c 

levels.  The  primary  significance  of  the  deviatoric  compression  test 
is  that  the  MNS  remains  constant  while  increments  in  deviatoric 
stress  occur.  Therefore  the  effects  of  hydrostatic  and  deviatoric 
stresses  on  hardening  can  be  studied  by  conducting  CTC  test  and  DCT 
to  the  same  point  on  the  failure  surface.  Also,  the  change  in 


plastic  volumetric  strain  during  the  pure  deviatoric  branch  is  a  mea 
sure  of  shear-volumetric  coupling.  A  composite  plot  is  presented  in 
Figure  2.16  which  shows  the  difference  in  stress  strain  response  for 
DCT  and  CTC  tests  to  essentially  the  same  level  of  maximum  stress. 
(Note,  these  tests  were  conducted  at  a  confining  stress  slightly 


below  the  BDT) . 
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Figure  2.16  Composite  plot,  to  same  maximum  stress 
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2.7  Verification  Tests 

2.7.1  General.  A  series  of  12  complex  load  path  tests  (verification 
tests)  were  conducted  on  concretes  of  different  strengths  ranging 
from  nominally  2  ksi  to  7  ksi.  Six  tests  were  conducted  at  the 
Waterways  Experiment  Station  (WES),  three  at  the  University  of  Colo¬ 
rado,  Boulder,  and  three  at  the  University  of  Eindhoven,  Netherlands. 
Reasons  for  selecting  particular  load  paths  presented  and  key  fea¬ 
tures  of  the  stress  strain  response  for  each  test  will  be  discussed 
in  the  following  sections.  Tests  conducted  at  WES  will  be  discussed 
in  section  2.7.2,  tests  conducted  at  Colorado  in  section  2.7.3,  and 
tests  conducted  at  Eindhoven  in  section  2.7.4.  The  nomenclature  for 
the  test  is  as  follows:  The  first  two  letters  indicate  verification 
tests,  the  number  or  numbers  before  the  dash  indicate  the  nominal 
unconfined  strength  of  the  concrete  and  the  number  after  the  dash 
indicates  the  number  of  the  test. 

2.7.2  WES  Tests 

Tests,  VT6.5-1  through  VT6.5-5,  were  conducted  on  f ^  =  6.5-ksi 
concrete  in  the  40  ksi  cylindrical  chamber. 

2.7.2. 1  VT6.5-1.  The  load  path  for  this  test  consists  of  pure 

hydrostatic  and  pure  deviatoric  branches  so  that  the  volumetric 
deviatoric  coupling  can  be  evaluated.  The  load  path  (Figure  2.17a) 
was  designed  to  penetrate  the  failure  surface  near  the  BDT  (i.e.  con¬ 
fining  stress  4  ksi).  The  arrows  on  the  load  path  indicate  the 
loading  and  unloading  portions  of  the  test.  The  first  branch  (0-1) 
represents  pure  hydrostatic  compression.  However,  as  shown  in  Fig¬ 
ure  2.17a  for  this  test  the  axial  and  radial  increments  are  large 
enough  so  that  the  hydrostatic  and  deviatoric  branches  are  themselves 
made  up  of  stairstep  axial  and  radial  stress  increments.  The  signi¬ 
ficance  of  this  increment  size  will  be  discussed  in  Chapter  4. 

Points  of  transition  from  one  branch  to  the  other  are  marked  on  the 
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load  path  and  the  stress  strain  plots  in  Figure  2.17.  Notice  from 
Figure  2.17b  the  material  never  softens.  For  this  test  a  clear  maxi¬ 
mum  axial  stress  of  approximately  az  =  23  ksi  was  reached  at  an  axial 
strain  of  approximately  =  0.0125  in/in.  Several  interesting  obser¬ 
vations  can  be  made  from  the  plots  presented  in  Figure  2.17b.  First, 
up  to  an  axial  stress  of  approximately  =  15  ksi  the  lateral 
strains  are  very  small  and  one  might  compare  the  response  to  a 
uniaxial  strain  test.  Also  for  the  first  four  branches  the  slope 
(i.e.  Aoz/Vez)  of  the  hydrostatic  and  deviatoric  portions  tend  to  be 
progressively  decreasing  until  branch  4-5,  where  the  material  seems 
to  be  stiffening. 

2. 7. 2. 2  VT6.5-2 

The  load  path  for  this  test  was  a  proportional  path  (the  ratio 

of  lateral  to  axial  stresses  are  kept  constant  throughout  the  test) 

designed  to  pierce  the  maximum  strength  surface  near  the  BDT  as  was 

the  case  in  VT6.5-1.  The  load  path  and  stress  strain  response  for 

this  test  are  shown  in  Figure  2.18.  For  this  test  a  clear  maximum 

axial  stress  of  approximately  a z  =  23  ksi  (the  same  as  test  6.5-1) 

was  reached  at  a  corresponding  axial  strain  of  =  0.03  in/in. 

After  reaching  maximum  stress  the  material  softens  until  the  loading 

was  stopped  and  unloading  began  at  an  axial  strain  of  approximately 

e  =0.05  in/in. 
z 

2. 7. 2. 3  VT6.5-3 

The  load  path  for  this  test  consists  of  stairstep  hydrostatic 
deviatoric  branches  just  as  was  the  case  for  VT6.5-I.  However,  the 
intent  here  was  to  penetrate  the  failure  surface  at  a  higher  mean 
normal  stress  level,  to  investigate  this  type  of  loading  in  the  hard¬ 
ening  region  (ie  above  the  BDT).  Again  each  branch  is  made  up  of 
independent  increments  of  axial  and  radial  stress.  Corresponding 
points  of  transition  from  hydrostatic  to  deviatoric  branches  are 
shown  in  Figure  2.19a  and  Figure  2.19b.  The  maximum  axial  stress 


Sigma  Z,  psi  Sigma  Z,  psi 


43 


attained  was  approximately  =  52  ksi  at  an  axial  strain  of  approxi¬ 
mately  o^  =  0.046  in/in.  Recalling  the  observations  made  for  test 
VT6.5-1,  the  lateral  strains  here  are  also  small  up  to  an  axial 
stress  of  approximately  a =  37  ksi.  Also,  the  slopes  (Ac^/Ve^)  of 
the  hydrostatic  and  deviatoric  branches  are  progressively  decreasing 
up  to  branch  4-5  where  a  significant  stiffening  occurs. 

2. 7. 2. 4  VT6.5-4 

The  load  path  for  this  test  (Figure  2.20a)  was  another  propor¬ 
tional  load  path  designed  to  pierce  the  failure  surface  near  the 
point  where  the  load  path  of  test  VT6.5-3  was  observed  to  pierce  the 
failure  surface.  In  this  test  an  unload-reload  cycle  was  executed  at 
an  axial  stress  of  approximately  a =  37  ksi.  There  is  no  clear 
maximum  or  failure  stress  for  this  test  as  the  stress  strain  response 
was  continuing  to  harden  at  an  axial  stress  of  approximately 
=  48  ksi  and  corresponding  axial  strain  of  e z  =  0.059  in/in. 

2. 7. 2. 5  VT6.5-5 

The  load  path  for  this  test  (Figure  2.21a)  consisted  of  a  pure 
hydrostatic  branch  to  a  mean  normal  stress  of  approximately,  MNS  = 

26  ksi  followed  by  a  radial  extension.  The  radial  extension  branch 
is  executed  by  holding  the  axial  stress  constant  and  decreasing  the 
radial  stresses.  This  test  was  designed  to  investigate  material 
response  and  stability  when  approaching  the  yield  and  failure  enve¬ 
lopes  along  a  load  path  significantly  different  from  the  conventional 
triaxial  compression  or  pure  deviatoric  paths.  More  specifically 
this  load  path  was  designed  to  investigate  problems  which  might  occur 
due  to  instability  as  discussed  by  Sandler  [12]  when  load  paths 
approach  the  failure  surface  in  this  manner.  Key  points  of  interests 
are  identified  (i.e.  A  through  G)  on  both  the  load  path  and  the  cor¬ 
responding  stress  strain  curve.  For  the  branch  AB  the  test  ran  very 
smoothly  as  the  axial  load  was  easily  held  constant  while  the  radial 
load  was  decreased.  However,  at  point  B  to  F  (which  is  in  the 
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(a)  stress  path 
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(b)  axial  stress  versus  axial  and  lateral  strains 

Figure  2.20  WEST  test,  VT6.5-4,  proportional  load  path  in  the 
hardening  region,  with  stress  strain  response 
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yield-failure  region  as  shown  in  Figure  2.11)  where  softening  is 
occurring  the  axial  load  could  not  be  held  constant  even  when  the 
maximum  displacement  rate  was  manually  set  on  the  440-kip  loader. 

The  load  path  was  essentially  moving  along  or  tangent  to  the  yield  or 
failure  surface  during  this  branch.  At  point  F  stability  in  the 
loading  system  was  regained  and  reloading  along  a  proportional  path 
to  point  G  was  performed  after  which  unloading  along  a  CTC  path  was 
executed  to  conclude  the  test. 

2.7.3  UCB  Tests 

Tests  VT4-1,  VT4-2  and  VT4-3  were  conducted  in  the  Colorado 
cubical  cell  on  4-inch  concrete  cubes  with  an  average  unconfined 

t 

compressive  strength  of  f  =  3.65  ksi  and  average  elastic  modulus  of 
E  -  3  x  106  psi. 

2.7.3. 1  VT4-1 

Test  VT4-1  is  a  proportional  load  path  test  to  a  final  stress 
state  with  a  hydrostatic  component  of  8  ksi  and  a  deviatoric  com¬ 
ponent  of  4  ksi.  This  test  is  similar  to  the  proportional  load  tests 
VT6.5-2  and  VT6.5-4.  The  load  path  and  resulting  stress  strain  plots 
for  test  VT4-1  are  shown  in  Figure  2.22. 

2. 7. 3. 2  VT4-2 

Test  VT4-2  consists  of  three  branches  OA  (hydrostatic) ,  AB  (con¬ 
stant  axial  stress),  and  BO  (unloading).  The  load  paths  and  corre¬ 
sponding  stress  strain  plots  are  presented  in  Figure  2.23.  This  test 
was  designed  to  reach  a  final  stress  state  (before  unloading)  with  a 
hydrostatic  component  of  8  ksi  and  a  deviatoric  component  of  4  ksi. 
This  load  path  is  similar  in  the  loading  portion  to  the  load  path  of 
test  VT6.5-5.  The  X  (lateral)  strains  measured  in  this  test  (Fig¬ 
ure  2.23b)  during  load  path  branch  AB  appear  to  be  questionable.  One 
would  expect  the  lateral  strains  to  decrease  as  the  lateral  stresses 
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Figure  2.22  UCB  test,  VT4-1,  proportional  load  path  in  the 
hardening  region,  with  stress  strain  response 
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Figure  2.23  UCB  test,  VT4-2,  hydrostatic  load  path  then  radial 
expansion,  with  stress  strain  response 
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are  decreased  from  A  to  B,  This  observation  will  be  discussed  fur¬ 
ther  in  Chapter  4. 

2. 7. 3. 3  VT4-3 

Test  VT4-3  was  designed  to  reach  the  same  final  stress  state  of 
test  VT4-1  and  VT4-2  along  a  stairstep  hydrostatic-deviatoric  path. 
The  load  path  for  this  test  and  corresponding  stress  strain  plots  are 
presented  in  Figure  2.24.  This  load  path  is  similar  to  those  of 
tests  VT6.5-1  and  VT6.5-3. 

2.7.4  UEN  Tests 

Tests  VT5-1,  VT5-2  and  VT5-3  were  conducted  in  the  Eindhoven 
cubical  cell  (section  2.2.6)  on  nominal  four  inch  cubes  of  concrete 

T 

with  an  average  unconfined  compressive  strength  of  f  =5.2  ksi  and 

6  ^ 

initial  elastic  modulus  of  E  =  5  *  10  psi.  These  tests  were 
designed  to  simulate  plane  strain  conditions  in  the  x,  z  plane, 
therefore  ey  =  0  throughout  the  test.  These  tests  were  also  designed 
to  investigate  the  effects  of  the  minor  and  intermediate  principal 
stresses  on  strength,  ductility  and  softening  in  concrete. 

2. 7. 4.1  VT5-1 

Test  VT5-1  was  conducted  along  a  proportional  load  path  in  the 
x,  z  plane  as  presented  in  Figure  2.25  with  the  ratio  of  to  a 
maintained  at  approximately  0.05.  Stress  strain  plots  are  also  pre¬ 
sented  in  Figure  2.24. 

2. 7. 4. 2  VT5-2 

Test  VT5-2  was  conducted  along  a  proportional  load  path  in  the 
x,  z  plane  with  the  ratio  of  a^/a^  held  at  0.10.  Load  paths  and 
stress  strain  plots  are  presented  in  Figure  2.26. 

2. 7. 4. 3  VT5-3 

Test  VT5-3  was  conducted  with  =  0  while  still  maintaining 
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Figure  2.24  UCB  test,  VT4-3,  hydrostatic  deviatoric  load 
path  with  stress  strain  response 
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(b)  stress  in  z  direction  versus  strains  in  z  and  x  directions 

Figure  2.26.  UEN  test,  VT5-2,  plane  strain  in  x,  z  plane 
proportional  load  path,  with  stress  strain  response. 


53 


ey  =  0.  Load  paths  and  stress  strain  plots  for  this  test  are  pre¬ 
sented  in  figure  2.21. 

2.7.5  Test  VT2-1  was  conducted  in  the  WES  30-ksi  multiaxial  test 
device  (section  2.2.5).  The  concrete  used  for  this  test  had  nominal 
strength  and  modulus  of  2  ksi  and  3.000  ksi  respectively. 

2.7.5. 1  VT2-1 

For  test  VT2-1  an  attempt  was  made  to  hold  strains  in  the  z 

direction  to  zero,  while  conducting  a  proportional  load  path  test  in 

the  x,  y  plane  with  a  =  a  .  However,  the  strain  in  the  z  direction 

x  y 

was  monitored  with  an  LVDT  which  includes  some  system  compliance. 

When  this  signal  from  the  LVDT  was  used  as  the  servo  control  signal 
for  the  2,400-kip  loader  (section  2.2.3)  some  small  strains  actually 
occurred  in  the  specimen  before  the  loader  responds.  This  is  seen  in 
Figure  2.28  which  presents  plots  of  strain  as  measured  by  the  embed¬ 
ded  strain  gages.  The  load  path  and  stress  strain  plots  for  this 
test  are  presented  in  Figure  2.28. 
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Figure  2.28  WES  test,  VT2-1,  plane  strain  test  in  the  x,  y 
plane  proportional  load  path,  with  stress  strain  response 
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Chapter  3 

Constitutive  Equations 


3 . 1  General. 

The  primary  objective  of  this  chapter  is  to  present  the  theoret¬ 
ical  concepts  and  mathematical  formulations  for  the  Fracture  Energy 
Based  Model  (FEBM)  and  the  Endochronic  Concrete  Plasticity 
Model  (ECPM) .  An  attempt  will  be  made  to  present  and  discuss  the 
motivations  for  different  concepts  used  in  each  model.  In  Section 
3.2,  the  development  of  constitutive  equations  for  linear  elasticity 
is  presented.  This  review  is  presented  primarily  to  emphasize  the 
significant  impact  of  simple  material  properties  tests  results  on  the 
development  of  sound,  rational,  and  consistent  constitutive  equa¬ 
tions.  A  general  discussion  of  conventional  plasticity  is  presented 
in  Section  3.3  with  again  emphasis  on  sound  rational  assumptions  in 
generalizing  the  results  of  simple  material  properties  tests.  Mate¬ 
rial  stability  postulates  and  uniqueness  requirements  are  presented 
in  Section  3. A  with  implications  for  strain-softening  models.  The 
Fracture  Energy  Based  formulation  is  presented  in  Section  3.5.  And 
finally,  endochronic  plasticity  theory  is  presented  in  Section  3.6 


3 . 2  Linear  Elasticity. 

Equations  of  equilibrium  and  compatibility  along  with  defini¬ 
tions  of  stress  and  strain  tensors  are  effectively  developed  in  the 
study  of  mechanics  of  continuous  media.  These  equations  and  tensor 
definitions  are  applicable  to  all  materials  which  can  be  represented 
with  sufficient  accuracy  as  a  continuous  body.  Robert  Hooke  [13]  in 
1676  presented  the  first  rough  law  of  proportionality  between  stress 
and  strain  which  is  today  commonly  referred  to  as  Hooke’s  Law. 
Hooke’s  Law  was  based  on  the  observation  that  many  engineering  mate¬ 
rials  when  stressed  in  one  direction  will  deform  in  that  direction 
and  the  relation  between  the  stress  and  strain  is  linear  and  elastic 
(i.e.  the  strain  returns  to  zero  when  the  stress  is  removed).  This 
response  can  be  presented  graphically  as  shown  below  for  the  region 
OA  of  Figure  3.1. 


Figure  3.1  One-dimensional  test  results 

This  curve  is  typical  for  many  materials  subjected  to  a  homogeneous 
stress  state  in  one  dimension.  Mathematically  the  linear 


relationship  between  stress  and  strain  in  the  region  OA  can  be 
written 
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a  =  E  e  3.1 

Where  a  is  axial  stress,  e  the  corresponding  axial  strain,  E  is 
Young's  modulus,  and  the  point  A  is  referred  to  as  the  proportional 
limit  of  the  material.  The  region  ABC  in  the  figure  will  be  dis¬ 
cussed  in  later  sections.  Since  the  state  of  stress  and  strain  is 
completely  determined  by  the  stress  tensor  o  and  strain  tensor  e 
and  since  the  one  dimensional  relationship  between  stress  and  strain 
is  linear  and  elastic  up  to  the  proportional  limit,  a  natural  gen¬ 
eralization  of  Hooke's  Law  is  obtained  by  assuming  a  one  to  one 
analytic  relationship  between  a  and  e  , 


a  =  F  (e)  3.2 

Since  F  is  analytic  it  can  be  expanded  in  a  power  series  in  terms  of 
e  as 


F=D  +  D,  e  +  D 0  e  e  + 
~  ~o  ~1  ~  ~2  ~  ~ 


D  (e) 
~n  ~ 


n 


3.3 


Further,  if  we  assume  the  initial  state  of  the  material  is  stress  and 
strain  free  and  assume  small  strain  theory  only  the  second  term  in 
(3.3)  needs  to  be  retained.  Therefore,  for  linear  elastic  small 
strain  theory  equation  (3.1)  is  generalized  to 


a  =  F  (e)  =  D  e 


3.4 


in  index  notation  this  is  expressed  as 


a 


ij 


=  Dijk£  ‘la  (i,j  ,k,  l  =  1,2,3) 


3.5 
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It  was  Cauchy  who  first  made  this  generalization  of  Hooke's  Law.  If 
the  material  is  homogeneous  the  coefficients  Dijki  will  be  indepen¬ 
dent  of  the  location  of  the  point  and  it  is  easily  shown  that 

Dijk£  =  Djiki 

|  and 

Dijki  =  Dij  £k 

due  to  the  symmetry  of  aij  and  the  symmetric  and  skew  symmetric 
decomposition  of  Dijki  (See  Sokolnikoff  [14]).  This  reduces  the 
number  of  independent  constants  in  Dijki  to  36.  The  existence  of 
the  strain  energy  density  function,  which  is  based  on  thermodynamic 
arguments,  reduces  the  number  of  independent  constants  to  21  for  the 
most  general  case  of  an  anisotropic  elastic  material.  If  the  elastic 
properties  of  the  material  are  equal  in  all  directions  it  is  said  to 
be  isotropic  and  the  number  of  independent  elastic  constants  reduces 
to  2  and  the  constitutive  equation  can  be  written  as 

crij  =  A  6  ij  +  2p  eij  3.6 

where  A , p  are  Lame7  constants.  In  this  study  only  assumed  iso¬ 
tropic  materials  will  be  considered.  The  assumptions  made  for  F  in 
equation  (3.2),  (i.e.  linear,  one  to  one,  analytic,  and  power  series 
expansion)  require  the  stress  strain  relationship  to  be  reversible. 
Contraction  of  (3.6)  leads  to 

oii  =  (3A  +  2u) 

which  when  substituted  in  (3.6)  leads  to 

>6  ij  °kk 
2u  (3  A  +  2y) 


3.7 
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The  Lame  constants  in  Equations  (3.6)  and  (3.7)  can  be  determined 
from  simple  test  results  as  follows: 


1.  Simple  tension  (or  compression):  all  stresses  except  o  are 


zero,  from  (3.7) 


X_  U11 
2y  (3X  +  2y) 


=  +  Hi  „ 

y(3X  +  2y)  11 


e22  “  e33 


2y(3X  +  2y)  °11 


recalling  from  Hooke's  law  the  ratio  of  stress  to  strain  is  defined 
as  Young's  modulus.  The  coefficient  of  <j^  in  (3.8)  can  be  set 
equal  to  1/E  such  that 


v  =  p(3X  +  2y )  . 
(X  +  y) 


Poisson's  ratio  is  defined  as  v  = 


e22  e33 


2(X  +  y) 


The  Lame  constants  are  often  expressed  in  terms  of  E  and  v  as 


(1  +  v)  (1  -  2v) 


2(1  +  v)  ’ 
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2.  Pure  hydrostatic  compression:  o,  =  o„„  =  o__  =  — —  =  a  with 

II  LL  33  3  m 

all  other  oij  =  0.  From  this  test  the  bulk  modulus  (K)  can  be  deter¬ 
mined  as  mean  stress  per  unit  volume  change 


a  =  K 
m 


e 


kk 


K  = 


3.  Pure  shear,  o  =  o?  ,  with  all  other  oij  =  0.  From  this  test 
the  shear  modulus  (G)  can  be  determined  as  shear  stress  per  unit 
change  in  shear  strain. 


a 


12 


G  yi2 


G  =»  p 


4.  Uniaxial  strain,  is  continuously  increased  while  holding  the 

lateral  strains  at  0,  i.e.  z^  =  z^  =  0.  From  this  test  the  con¬ 
strained  modulus  (M)  is  determined  as  the  ratio  of  0^  to 


M  =  (a  +  2u)  . 

Quite  often  in  the  development  and  evaluation  of  constitutive 
models  it  is  useful  to  consider  the  shear  (deviatoric)  components  of 
the  stress  and  strain  tensors  independent  from  the  hydrostatic  (volu¬ 
metric)  components.  For  an  isotropic  linear  elastic  material  this 
can  be  done  by  defining  the  deviatoric  stress  tensor  Sij  and  devia¬ 
toric  strain  tensor  eij,  as 
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Sij  =  oij  -  6ij  om 

and 

eij  =  eij  -  6ij  c 

m 


where  e  =  e,  .  /3 . 
m  kk 

As  previously  shown  contraction  of  (3.6)  leads  to 


3o  =  (3X  +  2G)  e.  , 
m  kk 


=  3  K 


e 


kk 


m 


=  K 


'kk 


3.9(a) 


which  is  the  hydrostatic  or  volumetric  component  of  (3.6).  By  sub¬ 
stituting  equation  (3.6)  into  the  deviatoric  stress  expression  along 
with  the  above  equation  for  the  deviatoric  or  shear  component  of 

the  constitutive  equation  can  be  expressed  as 


Sij  =  2G  eij  3.9(b) 

Either  equations  3.6  or  3.9  can  be  used  to  solve  constitutive 
problems  for  linear  elastic  isotropic  materials.  Given  the  strain 
tensor  the  stress  tensor  is  uniquely  defined  or  conversely  given  the 
stresses  the  strains  are  uniquely  defined.  Also,  for  finite  element 
applications,  the  same  solutions  can  be  obtained  in  rate  or  incre¬ 
mental  form  as 


=  D  c 


or 
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The  discussion  presented  in  this  section  and  derivations  of 
equations  can  be  found  in  many  texts  on  elasticity  such  as 
Sokolnikoff  [14],  Chen  and  Saleeb  [15],  etc.  The  main  reason  the 
derivation  is  presented  here,  is  to  point  out  that  in  developing  the 
governing  three-dimensional  constitutive  equations  of  elasticity,  a 
simple  one  dimensional  test  was  evaluated  (Hooke’s  Law),  the  mathe¬ 
matical  description  of  the  test  was  generalized  based  on  sound  mathe¬ 
matical  assumptions,  and  the  number  of  independent  material  constants 
required  was  determined  based  on  arguments  of  thermodynamics  and 
material  symmetry.  Also,  four  material  properties  tests  were 
described  which  provide  results  that  can  be  used  to  determine  key 
response  features  of  the  material.  It  is  not  unreasonable  to  expect 
that  more  tests  will  need  to  be  defined  and  conducted  to  determine 
key  response  features  of  materials  in  nonlinear  and  plastic  regions. 
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3.3  Conventional  Plasticity 

In  the  previous  section  a  brief  review  of  the  derivation  of 
linear  elastic  constitutive  equations  was  presented  for  materials 
which  can  be  reasonably  approximated  as  homogeneous  and  isotropic. 
This  derivation  was  based  on  the  results  of  a  simple  one  dimensional 
material  properties  test,  sound  mathematical  assumptions,  and 
rational  application  of  material  symmetry  restrictions.  In  this  sec¬ 
tion  a  derivation  of  the  equations  of  conventional  plasticity  will  be 
presented.  In  general  the  same  approach  will  be  used  here  as  was 
used  for  the  linear  elastic  case  in  terms  of  generalizing  the  results 
of  a  simple  test  based  on  rational  and,  consistent  assumptions  of 
mathematics  and  mechanics.  However,  as  is  well  known,  the  derivation 
is  not  quite  so  simple  and  direct  as  for  the  linear  elastic  case,  due 
in  part  to  the  need  to  postulate  a  specific  yield  surface,  define  a 
direction  for  plastic  strain,  and  determine  loading  and  unloading 
criteria. 

The  development  of  equations  in  this  section  essentially  follows 
that  presented  in  Martin  [16]  and  it  is  assumed  that  the  reader  is 
generally  familiar  with  plasticity  theory.  Only  material  nonlin¬ 
earity  is  considered  so  that  strains  and  rotations  are  small.  Fur¬ 
thermore,  emphasis  is  placed  on  rate  independent  hardening  plasticity 
where  there  is  no  flow  (I.e.  increasing  strain  at  constant  stress). 
Returning  to  the  simple  one  dimensional  test  shown  in  Figure  3.1, 
when  a  program  of  loading  causes  the  stress  to  reach  point  A  (yield 
point  or  proportional  limit)  initial  yielding  of  the  material  occurs. 
Up  to  this  point  (A)  if  the  loading  program  causes  the  stress  to 
decrease,  unloading  will  occur  along  the  original  linear  elastic  path 
at  constant  modulus  E  and  no  plastic  strains  are  accumulated.  If  on 
the  other  hand  the  loading  program  continues  to  cause  the  stress  (or 
strain)  to  Increase,  hardening  occurs,  and  plastic  (irrecoverable) 
strains  begin  to  accumulate.  The  definition  of  plastic  loading  is 
taken  here  to  mean  that  plastic  strains  Increase  during  plastic  load¬ 
ing.  If  in  the  hardening  region  (AB  in  Figure  3.1)  the  loading 
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Loading-unloading  response  of  a  hardening  material 


program  begins  to  cause  the  stress  (strain)  to  decrease,  unloading 
occurs  along  a  linear  elastic  path  with  modulus  E  as  shown  in  the 
figure  above.  The  material  represented  by  the  figure  above  is 
referred  to  as  a  hardening  material  because  when  unloading  occurs 
(say  at  point  I)  followed  by  reloading  a  linear  elastic  response  is 
observed  back  up  to  point  1  where  the  stress  is  greater  than 

o^,  the  initial  yield  stress.  Therefore  we  might  say  that  the 
plastic  loading  program  from  A  to  1  hardened  (or  pushed  up)  the  yield 
point  of  the  material  from  to  .  An  important  feature  of  the 

stress  strain  response  presented  is  that  once  the  material  is 
stressed  beyond  point  A,  a  knowledge  of  the  loading  history  is 
required  in  order  to  predict  the  stress  corresponding  to  a  particular 
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strain,  or  to  predict  the  strain  corresponding  to  a  particular 
stress.  The  material  is  said  to  be  history  dependent  and  a  portion 
of  its  plastic  loading  history  must  be  recorded  for  reference.  If 
the  material  is  never  stressed  beyond  point  A  there  is  no  need  to 
record  the  history  of  loading  since  elastic  response  is  not  history 
dependent.  By  observing  the  response  features  in  the  above  figure  we 
might  try  to  select  parameters  (hardening  parameters)  which  could  be 
used  to  record  the  history  of  loading  so  that  given  the  current 
strain  we  can  predict  the  corresponding  stress  or  vice  versa.  One 
hardening  parameter  might  be  a  record  of  the  highest  previously 
attained  stress  during  the  program  of  plastic  loading.  For  example 
if  the  highest  recorded  stress  is  a  and  the  current  strain  is 

then  the  current  stress  is  ^ ,  and  the  plastic  strain  is  . 

Another  choice  of  a  hardening  parameter  might  be  total  plastic 

p 

strain.  For  example  if  the  total  measured  plastic  strain  is  e ^ , 

and  the  current  strain  is  e  then  the  stress  is  again  calculated  as 

A 

ai  ' 

The  concept  of  yield  point  (i.e.  proportional  limit)  in  the  one- 
dimensional  test  can  be  generalized  directly  to  a  yield  hyper-surface 
in  stress  space.  The  region  bounded  by  this  surface  is  the  elastic 
region  and  when  the  stress  path,  corresponding  to  a  leading  program, 
lies  inside  this  region  there  is  no  accumulation  of  plastic  strain. 
The  hardening  concept  is  generalized  by  making  the  yield  surface  a 

function  of  stress  a  and  parameters  H  (a  =  l,2,....n)  which 

a 

represent  a  portion  of  the  loading  history  (recorded  history)  con¬ 
tributing  to  tite  position  of  the  yield  surface  and  plastic  strain. 

The  initial  yield  surface  is  therefore  defined  as: 

F  =  F  (a,  H  ) 
a 


Recalling  the  one-dimensional  test,  the  recorded  history  parameters 
ould  be  components  of  plastic  strain,  measures  of  work,  or  some 
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internal  variables,  because  one  could  derive  a  simple  mathematical 
model  for  this  test  and  successfully  predict  the  results  of  any  load¬ 
ing  program  as  long  as  the  total  history  of  leading  (e.g.  plastic 
strain)  is  recorded.  Figure  3.2  presents  the  different  regions  of 
stress  space  and  the  initial  yield  surface  as  discussed  above. 

I 


Figure  3.2  Yield  surface  and  elastic  region  in  stress  space 

It  is  clear  that  for  points  inside  the  yield  surface  F  <  0  ,  while 
for  points  on  the  yield  surface  F  =  0  ,  and  for  points  outside  the 
vield  surface  F  >  0  .  Due  to  the  Incremental  nature  of  plasticity 
we  are  ini  .ested  in  small  changes  in  stress  and  strain.  If  a  pro¬ 
gram  of  loading  has  resulted  in  the  stress  point  being  on  the  initial 
yield  surface  then 

F  (0,  H  >  „  -  0 

a 

If  the  loading  program  continues  and  results  in  a  small  change  in 
stress  dc  such  that  the  new  stress  point  falls  inside  the  initial 


yield  surface  then 
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dF 


_  SF  n 

H  ~  ^  d2  <  0 
a 


and  unloading  is  said  to  occur.  If  on  the  other  hand  the  loading 
program  continues  and  results  in  a  change  in  stress  do  such  that 
the  new  stress  point  is  still  on  the  initial  yield  surface  then 


f  d2‘° 


and  neutral  loading  is  said  to  occur.  There  is  no  change  in  the 
recorded  history  during  unloading  or  neutral  loading.  Finally,  if 
the  continued  loading  program  results  in  a  small  change  in  stress  do 
such  that  the  new  stress  point  lies  outside  the  initial  yield  surface 
then 


dO 


do  >  0 


and  loading  is  said  to  occur.  It  is  important  to  note  that  we  have 
defined  criteria  for  unloading,  neutral  loading,  and  loading  based  on 
gradients  of  the  yield  surface  while  Ha  remains  constant.  Also,  we 
observe  that  stress  states  outside  the  yield  surface  cannot  be 
obtained  unless  there  Is  a  change  in  the  recorded  history  Ha  . 
Therefore  when  loading  occurs,  Ha  must  change  and  we  could  write 
the  following 

F (o,  Ha)  =  0  3.10 

and 

F(o  +  do,  Ha  +  dHa)  =  0 

Equations  (3.10)  imply 

dF  =  —  do  +  dHa  =  0 

An  ~  SHa 


3.11 
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Equation  (3.11)  is  often  referred  to  as  the  consistency  condition, 
and  as  will  be  shown  later  provides  for  a  smooth  transition  between 
loading  and  unloading  conditions.  Equations  3.10  and  3.11  play  a 
very  important  role  in  the  development  of  plastic  constitutive  equa¬ 
tions.  Note  these  equations  are  differential  in  form  and  this  must 
be  kept  in  mind  when  incremental  solutions  of  the  equations  are 
sought.  Figure  3.3  presents  graphically  what  is  taking  place  in 
stress  space  under  a  loading  increment.  As  the  initial  yield  surface 
is  "pushed  out"  during  loading  the  new  yield  surface  is  referred  to 
as  a  subsequent  yield  surface  or  a  loading  surface. 


Figure  3.3  Subsequent  yield  or  loading  surfaces 


In  order  to  complete  the  development  of  a  general  plasticity 
theory  we  must  determine  the  functional  relationship  between  the 
increment  of  plastic  strain  and  stress,  hardening  parameters,  and 
increments  of  stress.  Since  we  require  that  plastic  strains  do  not 
increase  unless  the  recorded  history  increases,  it  is  reasonable  to 
assume  that  changes  in  the  recorded  history  dHa  are  the  result  of  a 
stress  increment  do  and  can  depend  on  the  current  stress  o  and 
recorded  history  Ha  .  Therefore,  we  can  write  (in  index  notation) 

dc^  =  h.  (do,  Ha)  dHa 
i  la 


3.12 
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This  general  representation  satisfies  the  requirement  that  plastic 
strains  do  not  accumulate  unless  the  recorded  history  changes.  Also 
since  we  are  concerned  here  with  hardening  behavior  only  (i.e.  no 
flow)  no  terms  can  appear  in  (3.12)  which  are  independent  of  do  and 
furthermore  since  we  consider  infinitesimal  do  ,  terms  of  second  and 
higher  order  in  do  can  be  neglected.  Therefore  we  can  assume  that 
de^  is  homogeneous  and  degree  one  in  components  of  do  .  Recalling 
the  one  dimensional  test,  it  is  reasonable  to  assume  that  do  con¬ 
trols  the  magnitude  of  dHa  ,  which  can  be  argued  by  considering 
equation  (3.11)  with  only  1  recorded  history  parameter. 


dH 


1  3F  do 
3F/9H  3  a 


Also  it  can  be  assumed  that  the  direction  of  dHa  is  a  function  of 

3  F 

o  ,  Ha  ,  and  not  do  .  Since  in  the  above  equations  —  do  is  in 

~  ~  0  <J  ~ 

effect  a  dot  product  or  scalar  quantity  and  therefore  do  contributes 
only  to  the  magnitude  of  dH  .  Therefore  we  can  express  dHa  as 


dHa  -  b  ta  (o ,  Ha) 


3.13 


where  b  can  be  evaluated  by  substituting  (3.13)  into  the  consis¬ 
tency  condition.  When  this  result  is  substituted  back  in  (3.12)  we 
have 


de 


P 

i 


hiq  tq 
3B6 


3.14 


Since  we  have  assumed  hia  ,  ta  and  3F/3Ha  each  depend  only  on  the 
state  variables  o  ,  Ha  equation  (3.14)  can  be  written 


de 


P 

i 


f i(o ,  Ha) 


dok 


71 


since  fi(o,  Ha)  is  a  vector  valued  function  of  a  ,  Ha  we  can 
express  this  vector  as 


3g(0,  Ha) 

fi(a,  Ha)  =  G(0,  Ha) - — —  3.15 

here  G(a,  Ha)  can  be  considered  a  scalar  hardening  coefficient  and 
g(0,  Ha)  can  be  considered  the  plastic  potential  function. 

In  summary  we  have  generalized  the  one-dimensional  constitutive 
equation  to  predict  the  hardening  plasticity  behavior  of  suitable 
materials  loaded  in  three  dimensions.  At  this  point  it  is  important 
to  note  that  the  development  of  the  above  constitutive  equations  was 
strongly  influenced  and  probably  dominated  by  observations  from  tests 
on  polycrystalline  metals  and  not  concrete.  The  assumption  that  the 
direction  of  the  plastic  strain  increment  should  be  independent  of 
the  stress  increment  do  is  primarily  based  on  conditions  required 
for  the  development  of  a  glide-system  of  crystal  grains.  Plastic 
strains  result  when  such  a  system  is  activated,  and  a  certain  stress 
state  is  required  to  develop  the  system.  The  required  stress  state 
is  a  macroscopic  stress  state  which  is  equivalent  to  some  statistical 
average  of  individual  grain  shears  along  critical  slip  directions. 
Once  the  stress  state  for  the  glide  system  is  developed,  along  with  a 
preferred  direction  for  plastic  strain,  the  stress  increment  do  only 
contributes  to  the  magnitude  of  the  plastic  strain.  Based  on  the 
above  arguments  and  the  observations  that  very  little  plastic  volume 
change  is  observed  during  plastic  loading  the  plastic  strain  incre¬ 
ment  direction  is  assumed  independent  of  do  and  the  plastic  poten¬ 
tial  is  assumed  independent  of  hydrostatic  stress.  The  assumption  of 
no  plastic  volume  strain  is  not  appropriate  for  concrete  and  the  role 
of  the  stress  increment  in  affecting  plastic  strain  direction  will  be 
discussed  in  Chapter  4.  Finally,  the  basic  equations  are,  in  incre¬ 
mental  form  using  index  notation: 
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strain  decomposition 


linear  elasticity 


de,  =  de®  +  de? 
iix 


de®  =  Cij  oj 

where  Cij ,  is  the  compliance  matrix  reduced  to  the  isotropic,  linear 
elastic  form  with  constant  coefficients. 

yield  surface 


F  <  0 


where  the  inequality  sign  is  used  since  we  are  considering  all  acces¬ 
sible  states  of  stress.  Finally,  the  plastic  strain  increment  de^ 
is  determined  from  the  following  conditions 


is  M 

3a,  9a, 

i  k 


0 

for 

F(a ,  Ha)  <  0 

elastic  region 

for 

F(a,  Ha)  =  0 

unloading 

0 

and 

9F 

9a  d2*° 

or  neutral  loading 

for 

F (a,  Ha)  =  0 

loading 

and 

H  d2*° 

oi  neutral  loading 

The  inequalities  are  included  for  both  unloading  and  loading  condi¬ 
tions  due  to  the  consistency  condition,  which  as  previously  pointed 
out  provides  for  a  smooth  transition  between  the  behavior  in  loading 
and  unloading.  Change  in  recorded  history 
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dHa 


ta  9F 

1L-  3°k  k 

6  », 


where  ta  =  ta  (a,  Ha) 


For  the  case  of  plastic  flow  we  must  allow  change  in  plastic 
strains  while  the  stress  remains  constant.  However,  this  is  in  con¬ 
flict  with  our  previous  assumption  that  plastic  strains  cannot  occur 
as  long  as  the  yield  surface  remains  fixed  which  implies  no  change  in 
recorded  history  parameters  Ha  .  This  problem  can  be  resolved  by 
recognizing  that  Ha  are  only  part  of  the  total  internal  parameters 
which  affect  plastic  strains.  These  additional  internal  parameters 
are  defined  here  as  JB  (B  =  1 , . . .m)  such  that 

de?  =  h.  dHa  +  h.Q  dJ 
3  ja  36  B 

where  h._  =  h.c  (03,  Ha,  J  ) 

3  p  3  p  p 

From  a  computational  standpoint  flow  can  be  allowed  by  defining  a 
maximum  strength  surface  which  bounds  all  loading  surfaces.  When  the 
stress  point  tries  to  move  outside  this  surface  it  is  forced  back  by 
some  return  procedure  to  a  new  stress  point  which  is  still  on  the 
surface  and  the  associated  plastic  strains  are  calculated. 

Quite  often  the  plastic  strain  increment  is  written  as 


where  X  is  a  scalar  hardening  parameter  defined  as 

•  3F 

X  =  G  (a,  Ha)  do 
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Recalling  the  consistency  condition  (Equation  3.11), 


3F  3F 

dF  =  r—  do  +  dH  =  0 

3o  ~  3R 


here  it  is  assumed  that  H  is  a  scalar  function  of  the  plastic 
strain  vector  £P  ,  we  can  solve  for  the  scalar  hardening  parameter. 
Substituting  the  expressions 

do  =  E  d£S 

3H  p 
dR  =  -  der 


in  the  consistency  condition  results  in 


3F  _  ,  ,  p,  ,  3F  3H  p 

5?  -  (d£  '  d£  >  +  3B  7?  dE  0 


and  since 


dEP  ’  *  If 


an  expression  for  the  scalar  hardening  parameter  can  be  obtained  as 


If  ms 


3F  3H_  3_g 

3H  3ep  3o 


21  E  & 
3o  =  3o 


Recalling 


do  =  E  dt6  -  E  (dc  -  dEP) 


the  incremental  constitutive  equation  can  be  expressed 


II  E 

3  a  ~  3o 

~  “  3F  3H  ^  W 

3H  p  3o  3o 
d  £  ~  ~ 
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3.4  Material  Stability  Postulates,  Uniqueness. 

The  stability  postulates  discussed  in  this  chapter  are  not  pre¬ 
sented  as  thermodynamic  or  physical  requirements  of  real  materials. 
These  postulates  are  often  used  as  a  basis  for  developing  constitu¬ 
tive  equations  for  certain  classes  of  materials.  Returning  to  the 
one  dimensional  test  as  shown  below,  we  observe  that  a  monotonic 
increase  in  stress  produces  a  monotonic  increase  in  strain  in  the 
region  0  A  B. 


B 


One-dimensional  test  results 


This  can  be  expressed  mathematically  as 

(a2  -  a1)  (c2  -  c1}  >0  3.16 

and  this  inequality  can  be  taken  as  the  meaning  of  work  hardening  in 
the  one  dimensional  case.  The  equal  sign  is  provided  to  allow  for 
infinite  or  0  slope.  Equation  3.16  is  simply  based  on  the  observa¬ 
tion  that  the  assumed  homogeneous,  isotropic  material  when  subjected 
to  a  one  dimensional  monotonic  increment  cf  stress  responds  with  a 
monctonic  increment  in  strain.  Inequality  (3.16)  can  be  generalized 
to  represent  multi-dimensional  states  of  stress  and  strain  in  the 
following  way.  A  monotonic  change  in  a  single  stress  component 
results  in  a  straight  line  path  in  stress  space,  and  we  require  the 
product  of  this  change  in  stress  with  the  associated  change  in  strain 
to  be  non-negative,  therefore 
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?  i  2  1 

(a:  -  at)  •  (eT  -el)  >0  3.17 

l  l  l  l 

Since  an  infinitesimal  stress  increment  is  applied  monotonically  and 
results  in  an  infinitesimal  strain  increment  inequality  (3.17)  can  be 
expressed  as 

do  de  >  0  3.18 

Another  feature  of  the  one  dimensional  test  is  that  the  unloading 
slope  in  the  region  AB  is  elastic  and  this  slope  is  greater  than  the 
slope  in  the  hardening  plastic  loading  region.  This  feature  can  be 
preserved  if  ic  is  required  that  the  complementary  work  around  a 
closed  cycle  in  stress  be  non-positive,  therefore 

e  do  S  0  3.19 

The  equal  sign  holds  for  stress  cycles  in  the  elastic  region.  If  the 
stress  cycle  produces  a  change  in  plastic  strain  Ae^  which  is 
infinitesimal  or  small,  inequality  (3.19)  is  referred  to  as  the 
second  postulate  in  the  weak  form.  If  no  restrictions  are  placed  on 
the  size  of  Ae^  inequality  (3.19)  is  referred  to  as  the  second 
postulate  in  the  strong  form.  Drucker  [17]  also  identifies  work 
hardening  for  the  one  dimensional  tests  with  monotonic  Increments  in 
stress  producing  monotonic  increments  in  strain.  However,  in  gen¬ 
eralizing  this  feature  Drucker  considers  the  work  done  by  an  external 
agency  which  slowly  applies  load  then  slowly  removes  the  load.  The 
external  agency  is  completely  Independent  of  the  loading  program 
which  produces  the  existing  state  of  stress  and  strain  in  the  mate¬ 
rial.  Drucker  then  states  that  for  all  sets  of  added  stresses  (due 
to  the  external  agency)  work  hardening  implies  that  the  material  will 
remain  in  equilibrium  and: 


1.  The  external  agency  does  positive  work  during  the  applica¬ 
tion  of  load 

do  de  >  0 

2.  The  net  work  performed  by  the  external  agency  during  a 
cycle  (application  and  removal  of  the  load)  is 
non-negative . 

do  (de  -  de6)  =  do  de^  >  0 

As  previously  mentioned,  these  postulates  are  not  intended  to  be 
physical  or  thermodynamic  requirements  for  the  behavior  of  real 
materials  but  can  be  used  to  develop  rational  and  consistent  consti¬ 
tutive  relations  for  broad  classes  of  engineering  materials.  These 
stability  postulates  prohibit  the  types  of  stress  strain  responses 
shown  in  Figure  3.4. 


Figure  3.4  Examples  of  response  features  which  violates 
stability  postulates 

Figure  3.4(a)  of  course  is  the  response  featur.  referred  to  as  strain 
softening  and  is  often  observed  as  shown  previously  in  unconfined  and 
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low  confinement  compression  tests  on  concrete  and  other  materials. 

It  was  reasoned  by  Bazant  [10]  that  strain  softening  can  only  be 
observed  on  material  test  specimen  with  micro-inhomogeneities.  Often 
we  construct  a  constitutive  model  based  on  a  softening  stress-strain 
curve  (actually  developed  from  a  force  displacement  curve)  which 
assumes  homogeneity  of  stress  and  strain  when  significant  in  homo¬ 
geneities  actually  occur  (eg.  cracking,  shear  bonding,  etc.).  Con¬ 
stitutive  models  developed  in  this  way  should  not  be  expected  to 
effectively  reproduce  material  response  over  a  wide  range  of  complex 
loading  conditions.  Therefore,  it  is  expected  that  some  models  may 
provide  good  results  for  specific  load  cases  and  be  totally  inappro¬ 
priate  and  possibly  unstable  for  other  load  cases. 

The  stability  postulates  presented  above  are  considered  as 
mathematical  restrictions  on  constitutive  relations  in  plasticity 
theory.  These  restrictions  result  in  the  requirement  that  the  yield 
surface  be  convex.  At  a  singular  point  (corner,  vertex)  on  a  yield 
surface  or  at  the  intersection  of  two  or  more  surfaces  the  plastic 
strain  increment  must  be  within  the  fan  or  hypercone  bounded  by  the 
normals  to  the  respective  surfaces  at  the  singular  point.  Further¬ 
more,  for  hardening,  satisfying  the  stability  postulates  leads  to  a 
sufficient  condition  for  uniqueness  of  solution  in  the  incremental 
sense,  where  equilibrium  equations  and  kinematic  relations  are 
linear.  Since  yielding  and  other  nonlinear  response  features  of  con¬ 
crete  are  strongly  dependent  on  the  confining  stress  level,  different 
forms  of  loading  and  failure  surfaces  are  often  used  for  different 
stress  regions.  When  these  surfaces  are  combined  in  a  constitutive 
model,  intersecting  corners  can  occur.  If  the  loading  stress  path 
intersects  these  surfaces  at  a  corner  point,  obviously  there  exist 
more  than  one  surface  normal  at  this  point  ard  the  direction  of 
plastic  strain  is  not  uniquely  defined.  If  the  plastic  strain  incre¬ 
ment  at  such  a  singular  point,  which  is  the  intersection  of  say  m 
continuous  yield  surfaces,  is  taken  as  the  linear  combination  of  the 
m  plastic  strain  increments  for  the  react!  -e  surfaces,  then  the 
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solution  is  unique  and  the  stability  postulates  are  satisfied.  Here 
reactive  implies  those  surfaces  which  are  to  be  used  for  computing 
plastic  strain  increments.  Proofs  of  uniqueness  can  be  found  in 
Martin  [16]  Pages  109-114. 
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3.5.  Examples  of  Hardening  and  Softening  Formulations. 

Consider  a  homogeneous  isotropic  material  subjected  to  loading 
under  plane  stress  conditions  such  that  =  0*  Furthermore,  assume 
the  material  behaves  the  same  under  tension  and  compression  and  that 
the  Von  Mises  yield  condition  is  appropriate  for  describing  the  ini¬ 
tial  yield  and  subsequent  (loading)  yield  surfaces  for  the  material. 
Also,  let  the  one-dimensional  response  of  the  material  to  axial 
stress  be  as  shown  in  Figure  3.5. 


Figure  3.5.  One  dimensional  response 
of  example  hardening  material 


For  the  Von  Mises  surface  initial  yielding  will  occur  when  the  state 
of  stress  satisfies  the  condition 


+  o„ 


-  0l02  = 


3.20 


and  for  any  subsequent  yield  or  loading  surface  the  stress  state  will 
satisfy  the  condition 


F  = 


Cl°2 


-2 

c 


0 


3.21 


32 


where  =  the  largest  previous  value  of 

2  2 

cr  +  a  -  oo,  assuming  o  >  o  . 


Finally,  it  is  assumed  that  an  associative  flow  rule  is  appropriate 
so  that  it  only  remains  to  determine  G  (o^_. ,  Ha) ,  the  scalar 
hardening  parameter,  before  stress  strain  response  problems  can  be 
solved.  To  this  end  we  assume  that  the  hardening  is  isotropic  which 
means  that  G  is  constant  on  any  particular  yield  or  loading  sur¬ 
face.  An  important  result  of  this  assumption  is  that  the  hardening 
characteristics  for  the  material  can  be  completely  determined  from 
one  simple  test.  With  the  assumptions  made  above  the  basic  incre¬ 
mental  equations  for  the  case  of  plane  stress  are 


de 

de 


e 

1 

P 

1 


da. 


=  G 


3F 

dO 


1 


dF 


3.22 


It  remains  to  determine  G,  which  can  be  accomplished  in  the  fol¬ 
lowing  way  for  isotropic  hardening  (i.e.  G  has  the  same  value  at 
all  points  on  any  one  subsequent  yield  or  loading  surface).  Assume 
that  a  monotonic  one  dimensional  compression  test  is  conducted  on  the 
material  discussea  above  such  that  the  stress  strain  response  is  as 
presented  in  Figure  3.5.  Since  the  basic  equations  must  apply  to 
this  test  and  since  the  loading  is  monotonic,  we  can  write 


and 


0  <  a  <  oo 


6  P 

de  =  de  +  der 


o  >  oo 
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Therefore  from  equations  (3.22)  for  d£j  , 


where  here  =  a,  and  =  0  such  that 

de  =  [|  +  4C2  g]  do  3.23 


From  Figure  3.5  it  is  seen  that  the  relationship  between  increments 
of  stress  and  strain  can  also  be  written  as 


3.24 


During  this  monotonic  loading  condition  the  stress  o  continuously 
satisfies  the  yield  condition  3.21  so  that  F  =  0  and  consequently 
o=o.  With  this  in  mind  Equations  3.23  and  3.24  imply. 


3.25 


For  this  example  we  are  treating  compression  as  positive  and  only 
concerned  with  the  positive  quadrant  of  0  2  sPace*  From  3.25 

G(£,  H«)  is  a  function,  of  the  current  stress  state  (through  o)  and 
could  also  be  expressed  as  a  function  of  plastic  strain  .  There¬ 
fore,  given  Equations  3.22  and  3.25  we  can  solve  hardening  plasticity 
problems  under  monotonic  loading  given  a  particular  stress  path  in 
plane  stress  space  (o^a^).  From  Figure  3.5  we  can  write  E^  (o)  so 
that  the  integration  of  Equations  3.22  could  be  performed  with  rela¬ 
tive  ease.  For  a  particular  continuous  proportional  loading  path 
(say  =  20o)  the  effect  on  the  loading  surfaces  is  uniform 


expansion  as  shown  in  Figure  3.6.  The  reason  for  this  is  that  as 
loading  progresses  a  increases  and  G  is  always  positive  since 
E  <  E. 


Figure  3.6  Isotropic  hardening  of  yield  surfaces 


Based  on  the  above  assumptions  the  exact  solution  to  the  problem  of 
continuous  proportional  loading  (do^  =  2dc?2)  to  a  final  stress  state 
0^  =  2 Cq  can  be  obtained.  For  this  example  we  take  the  case  where 
E  =  constant. 

2/3  /n 

At  the  end  of  the  elastic  region  (ie;  =  —j—  00,  =  —j  00)  the  elastic 

strains  are  calculated  from  (3.22)  as 


€ 


e 

1 


e 


e 

2 


S3  00 
3E 


(2-v) 


S3  00 
3E 


( I-2v) 


For  the  remainder  of  the  stress  path,  the  total  strains  can  be  deter¬ 
mined  by  integration  as 
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ao  .  ,  .63450 
£1  =  ~ E  (2“V)  +  - E  E 


e2  =  i  (1-2V) 


E  -  ET 


The  main  point  to  be  made  here  is  that  for  hardening  plasticity  the 
solution  to  this  simple  problem  can  be  obtained  in  closed  form. 

If  on  the  other  hand,  the  one-dimensional  response  to  axial  loading 
is  observed  to  be  that  presented  in  Figure  3.7  several  different 
problems  arise  in  obtaining  the  solution  of  the  plane-stress  problem. 


Figure  3.7.  One  dimensional  response  of 
example  softening  material 


Here  it  is  assumed  that  the  softening  behavior  observed  in  the 
one  dimensional  test  (Figure  3.7)  is  generalized  to  an  isotropic 
softening  in  plane-stress  space.  The  qualitative  effect  of  the 
softening  on  subsequent  yield  surfaces  is  shown  in  Figure  3.8.  Since 
the  increment  in  stress  do  for 
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Figure  3.8  Isotropic  softening  of  yield  surfaces 

plastic  loading  in  the  softening  region  will  be  negative  our  previ¬ 
ously  defined  loading  and  unloading  criteria  will  not  be  appropriate. 

Also,  since  plastic  loading  continues  as  a  and  therefore  o 

-2 

decrease,  the  largest  previously  attained  value  of  o  cannot  be 
used  as  a  recorded  history  parameter.  Furthermore,  we  cannot  use  the 
smallest  previously  obtained  value  of  o  since  a  decreases  under 
elastic  unloading  also.  Therefore  it  appears  that  plastic  loading 
should  be  defined  here  as  Increase  in  plastic  strain.  Another  prob¬ 
lem  encountered  in  softening  materials  is  demonstrated  in  the  figure 
on  the  following  page. 

For  stress  space  plasticity,  a  is  the  independent  variable  and 
strains  are  solved  for  from  the  equation 


t 
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Hardening  material.  Softening  material. 

while  for  strain  space  plasticity  e  is  the  independent  variable  and 
stresses  are  solved  for  from  the  equation 

o  *  E  e  - 

where  oc  is  a  stress  relaxation  from  the  stress  which  would  be  pre- 
dieted  under  purely  elastic  response.  For  the  softening  problem, 
conceptually,  some  inconsistencies  are  indicated  by  the  stress  space 
formulation.  First  the  end  point  (e„  ,  oD)  cannot  be  reached  by  fol- 

D  D 

lowing  a  monotonic  load  path  and  since  the  strain  must  be  specified 
in  this  test  we  do  not  know  a  priori  what  the  final  stress  o.,  will 

V 

be.  The  stress  space  formulation  could  be  used  for  the  isotropic 
softening  problem  by  developing  a  softening  function  G  in  the  fol¬ 
lowing  way.  Specializing  Equations  3.22  to  the  softening  material 
shown  above  following  the  same  procedures  used  for  the  hardening 
material  results  in 
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Therefore  considering  only  the  plastic  strain  components  for  the 

2  2  -2 

plane  stress  problem  (where  F  =  o^  ->'*  O2  -  °1°2  =  0  )  we  have 


and 


2 (2o^  doj  +  2^2  da 2) 


2 (2aj  do^  +  202  do2) 


Here  again  the  problem  is  we  do  not  know  the  stress  path  as  a  func¬ 
tion  of  strain  in  the  plastic  region  and  we  cannot  obtain  a  general 
closed  form  solution  as  we  'did  in  the  case  of  the  hardening  material 
using  a  stress  space  formulation.  For  the  special  case  of  an  iso¬ 
tropic  material  and  a  symmetric-monotonic  plastic  loading  path  (ie  e. 

P  P  * 

=  e„)  to  a  final  plastic  strain  state  of  e,  =  =  e  a  solution 

2  r  1  2  o 

for  the  softening  material  problem  can  be  obtained.  For  this  case 
aI  =  02  ’  an<^  _<*ai  =  _<*°2  ^ue  to  Therefore  from  the  above 

equations 


and 


o 


do 


1 


from  which  the  final  stress  can  be  solved  for.  The  point  to  be 
made  here  is  that  an  additional  assumption  (symmetry)  had  to  be  made 
before  this  simple  problem  could  be  solved. 
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3.6  Fracture  Energy  Based  Model 

3.6.1  General 

The  Fracture  Energy  Based  Model  (FEBM)  was  developed  at  the  Uni¬ 
versity  of  Colorado  (3)  in  an  effort  to  provide  a  rational  and 
consistent  formulation  for  predicting  the  triaxial  response  charac¬ 
teristics  of  plain  concrete  to  include  strain  softening.  The  model 
is  based  on  conventional  plasticity  in  the  sense  that  it  makes  use  of 
initial  yield  and  loading  surfaces,  decomposition  of  strain  incre¬ 
ments  into  elastic  and  plastic  components,  a  flow  rule,  and  the  con¬ 
sistency  condition.  The  loading  surfaces  expand  with  increasing 
plastic  strain  and  hydrostatic  compressive  stress  and  rotate  about  a 
point  of  equal  tension  stress.  The  model  is  developed  for  both  asso¬ 
ciated  and  nonassociated  flow  assumptions.  The  strain  softening 
formulation  is  based  on  fracture  energy  concepts  and  the  criterion 
that  strain  softening  occurs  when  the  loading  stress  state  reaches 
the  failure  surface  at  a  confining  stress  less  than  the  confining 
stress  which  defines  the  brittle-ductile  transition. 

3.6.2  Failure  Surface 

The  failure  surface  used  in  the  FEBM  was  developed  by  Leon  (18), 
and  is  based  on  forcing  continuity  between  the  Mohr-Coulomb  friction 
law  and  the  tension  cut  off  condition  of  Rankine.  Mathematically  the 
two  parameter  surface  is  represented  as 

°1 

+  Mo  -tt  -  C  =  0  3.28 

i  o 

c 

where  =  Major  principal  stress 

=  Minor  principal  stress 

Mo  =  friction  parameter 

C  =  cohesive  parameter. 

o 

f^  =  Unconfined  Compressive  Strength 


When  0^  =  0  equation  3.28  yields 


=  Co 


so  that  taking  03  =  fc  for  this  case  results  in  C^  =  1  .  This  is 

the  same  result  found  for  the  Mohr-Coulomb  surface  where  C  would 

o 

be  defined  as  the  normalized  cohesion  of  the  material  while  the  fric 
tion  parameter  Mo  is  not  exactly  the  same  as  the  coefficient  of 
internal  friction,  which  is  defined  in  the  Mohr-Coulomb  model,  it  is 
a  measure  of  the  slope  of  the  failure  surface  (and  loading  surfaces 
to  be  discussed  in  the  next  section) . 

The  failure  surface  defined  by  Equation  (3.28)  and  shown  in  Fig 
ure  3.9  represents  the  maximum  strength  of  the  material  when  the  nor 
malized  cohesive 


Maximum  strength  . /k=l  .0 


Figure  3.9  Hardening  and  maximum  strength  surfaces  for  the  FEBM 
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parameter  Cq  =  1 .  Mo  can  be  determined  by  setting  =  0  and  defin¬ 
ing  the  ^alue  of  for  this  condition  as  =  f '  which  leads  to 


Mo  = 


3.29 


where  f ’  =  unconfined  tensile  strength 

3.6.3  Loading  Surfaces  (Hardening) 

The  loading  surfaces  for  the  FEBM  during  hardening  can  be 
expressed  as 


F  (oj. 


where 


*  k) 


(1 


k2  Mo 


k 


(v 

k 


C  =  0 
o 

=  hardening  parameter,  ko  <  k  <  1 
=  value  of  k  defining  the  initial 


°1  2 

jt  ~  ko  3.30 
c 


yield  surface 


Loading  surfaces  for  different  values  of  k  are  presented  in  Fig¬ 
ure  3.9.  Three  important  observations  can  be  made  regarding  the 
loading  surfaces. 

1 .  When  k  =  I  the  loading  surface  corresponds  to  the  failure 
surface  F. 

2.  The  hardening  behavior  is  not  purely  isotropic,  since  the 
loading  surfaces  essentially  rotate  about  a  common  point  in 
the  tension  region. 

3.  The  loading  surfaces  intersect  the  hydrostat  at  angles  other 
than  90°  and  other  corners  exists  at  60°  intervals  in  devia- 
toric  sections. 


In  order  to  simulate  the  pressure  sensitive  behavior  of  concrete  the 
hardening  parameter  is  expressed  as  a  function  of  plastic  strain  and 
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confining  stress.  The  development  of  the  expression  for  k  was 
motivated  by  the  fact  that  increasing  the  confining  stress  (o^  =  o^) 
increases  the  maximum  strength  (o^)  and  the  strain  (e^)  at  which 
maximum  strength  is  reached.  Therefore,  one  could  conduct  tests  at 
low,  medium,  and  high  confining  stresses  then  normalize  the  ordinate 
(i.e.  °3^  e^^03n,ax  =  k)  and  plot  k  versus  plastic  strain  (  eP)  up  to 
strains  where  maximum  strength  is  attained.  The  general  characteris¬ 
tics  of  these  curves  are  shown  in  Figure  3.10.  To  simulate  this 
variation  of  k  with  eP,  a  monotonically  increasing  elliptic  function 
of  eP  was  selected. 


k  (cP, 


p,  i  .  1  -  ko 
X  )  =  ko  +  - - — 


-  <cV 


3.31 


where 


X*5  ~  ductility  parameter,  the  plastic  strain  when  maximum 
strength  is  attained  for  a  given  confining  stress. 


"P 

f?  = 


=  equivalent  plastic  strain  rate. 


An  empirical  expression  was  developed  for  xP  which  is  also  calibrated 
from  low,  medium,  and  high  confining  stress  test  results.  This 
expression  is 


XP  =  ^ 


3.32 


Ideally  for  a  given  concrete  three  triaxial  compression  tests  can  be 

conducted  (eg.  low,  medium,  high  confining  stress  (a^))  and  three 

values  of  y?  ,  measured.  Then  knowing  f^  and  Oj  these  equations 

,  B  ,  and  C,  can  be  determined,  and  solved. 

H  H 


in  the  unknowns 
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3.6.4  Loading  Surfaces  (Softening) 

If  the  loading  stress  path  intersects  the  failure  surface  (Equa¬ 
tion  3.28)  below  the  brittle-ductile  transition  point  (TP)  and  load¬ 
ing  continues  (i.e.  Ae^  >  0)  strain  softening  is  predicted  by  the 
FEBM.  This  softening  response  is  accomplished  by  collapsing  the 
failure  surface  through  a  continuous  decrease  of  the  cohesive  parame¬ 
ter  until  a  residual  strength  surface  is  reached.  The  initial  and 
residual  strength  failure  surfaces  are  presented  in  Figure  3.11.  The 
cohesion  is  decreased  according  to  the  expression 


where 

f^  =  tensile  strength  of  concrete 
ot  =  degrading  tensile  strength 

Uj  =  crack  opening  displacement  measured  in  a  direct  ten¬ 
sion  test 

U  =  residual  crack  opening  displacement,  crack  opening  at 
fracture. 

As  shown  in  Figure  3.11  the  decrease  in  cohesion  is  accompanied  by  an 
increase  in  M  (the  frictional  parameter)  or  say  a  frictional  harden¬ 
ing.  In  effect  the  softening  surfaces  rotate  about  the  point  TP  in 
stress  space.  The  relation  between  M  and  C  is  assumed  as 

M  =  M  -  (M  -  Mo)  C  3.34 

r  r 

Where  is  the  value  of  the  frictional  parameter  at  residual 

strength.  The  cohesion  parameter  C  can  therefore  be  thought  of  as 
the  softening  parameter  and  can  be  determined  from  a  simple  uniaxial 
tension  test  (which  is  conducted  in  stroke  control).  The  results 
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\ 

Figure  3.11 


Rotation  of  softening  surfaces  about  TP 


from  a  direct  tension  test  along  with  predictions  from  equation  3.33 
are  presented  in  Figure  3.12,  (from  Reference  3).  The  motivation  for 
using  the  direct  tension  test  for  defining  softening  in  compressive 
stress  states,  was  based  on  the  idea  that  the  strain  softening  phe¬ 
nomenon  is  dominated  by  microcrack  growth  which  is  associated  with 
dilation  (expansion)  in  compressive  tests.  Considering  the  case  of 
direct  tension  a  strain  softening  plasticity  formulation  can  be 
developed  in  the  same  manner  as  for  hardening  plasticity.  Assuming 
the  initial  yield  and  loading  surfaces  can  be  defined  by  the  maximum 
tension  strength  surface 


Figure  3.12  Softening  in  a  direct  tension  test 


F(Oj ,  ot)  =  Oj  -  ot  =  0 


where  a  =  (e^) 

=  fracture  strain 

and  that  the  plastic  (fracture)  strain  in  tension  is  defined  in  rate 
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the  consistency  condition  is  given  by 


F  = 


3F  t  E  (c  -  A  3f\  / 3ot  SefX*  •  3F  _  - 

3o  .  V~  1  3o  /  \ 3ef  3e  j  f  3o  0 


3.35 


3F  E  e 
3o 


Xf  = 


Ef  + 


/3F\C  3F 
\3o/  -  3o 


where  the  softening  modulus  Ef  is  defined  as 


3.36 


Up  to  this  point  basic  concepts  of  plasticity  (i.e.  strain  rate 
decomposition,  normality,  consistency)  have  been  applied  to  a  strain 
softening  material.  However,  when  attempting  to  compute  from 
Equation  (3.36)  a  problem  arises  due  to  the  fact  that  ot  is  defined 
as  a  function  of  not  (Equation  3.33,  Figure  3.12).  The  reason 
for  this  is  that  in  a  direct  tension  test  (in  displacement  control) 
at  peak  stress  the  microcracks  coalesce  into  a  localized  band  and  the 
response  of  the  specimen  is  better  or  more  appropriately  defined  in 
terms  of  stress  versus  displacement  (crack  opening  displacement)  as 
shown  in  Figure  3.13  (from  Reference  3).  In  the  FEBM  this  problem  is 
addressed  by  using  the  chain  rule  such  that 


3o  3ot 


Therefore,  the  softening  modulus  can  be  expressed  as 


E 


f 


3ot  3uf  9cf  3F 
3e^  3c^  3e^  3o 


3.37 
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and  the  only  term  which  cannot  be  directly  computed  is  3U^/3e^.  It 
is  here  where  an  assumption  is  made  which  will  allow  for  the  exten¬ 
sion  of  the  strain 


Equivalent  Continuum 


Figure  3.13  Localization  versus  continuum  ideas,  from  (3) 

softening  theory  to  compressive  stress  states  but  also  results  in  a 
theory  which  when  implemented  in  a  finite  element  analysis  procedure 
will  exhibit  mesh  sensitivity.  Essentially  this  assumption  is  to 
smear  the  effect  of  the  localized  crack  band  over  the  volume  of  the 
test  specimen.  This  is  accomplished  by  expressing  the  softening  mod¬ 
ulus  (for  the  direct  tension  test  as 
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where  E,  =  3a  /3U. 

a  t  r 


=  3U^/3e^  ,  height  of  tast  specimen. 


©' 


For  the  direct  tension  test  h^  can  a  so  be  expressed  as  the  elemen¬ 
tary  volume  of  the  test  specimen  div:  ded  by  the  cracked  surface  a-  ea. 


t  A 


The  theory  is  ext  nded  to  comp  essiv*  stress  states  directly  by  using 
the  failure  surf a-  2  in  the  pos  peak  region 


cl  -  a3 

p 


-  C  =  0 


where  C  and  M  ire  defined  by  Equations  (3.33)  and  (3.34),  :spec- 
tively.  Also,  in  :he  softening  region  the  equivalent  plastic  rac- 
ture)  strain  increment  is  determined  from 


<i, 


iere  the  McC;  uley  brackets  <  >  extract  only  the  positive  components 
nsion)  of  ‘  re  inelastic  fracture  r train  rate  tensor  .  The 
,y  question  remaining  for  the  compressive  formulation  is  how  to 
t  termine  3U^,  3e^  for  compressive  stiess  states.  Here,  a  similar 
r  ',ument  as  u  ed  previously  for  the  case  of  direct  tension  is  pre- 
r  ited.  The  ssuraption  is  made  that  under  compressive  states  an 
e  uivalent  ch,  racteristic  lei gth  can  be  defined  Euch  that 
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v_ 

Ac 


he 


One  might  think  of  the  volume  fraction  A c/V  as  microcrack  density  in 
compression.  While  it  was  quite  clear  how  to  determine  A  in  the 
tension  test  it  is  not  so  obvious  how  Ac  (the  fracture  surface 
area)  should  be  determined  for  low  confinement  compression  tests. 
L'nder  unconfined  and  very  low  confinement  compression  tests,  vertical 
(splitting  tensile)  cracks  form  which  are  the  result  of  mode  I  type 
fracture  (ie  fracture  due  to  tension  only)  in  the  circumferpnt ial 
direction.  If  the  area  of  these  vertical  cracks,  which  ate  the 
result  of  splitting  compression,  is  defined  as  Ac  then  the  fracture 
energy  associated  with  this  fracture  can  be  expressed  as 


il  =  C*  AAc 
cs  f 


When  results  of  compressive  splitting  tests  are  compared  to  direct 
tension  tests  (see  Figure  3.14)  it  is  seen  that  the  fracture  energy 
release  rates  in  tension  and  splitting  compression  are  quite  similar 
As  the  confining  stress  is  increased  the  tensile  cracks  tend  to 
coalesce  into  a  single  shear  band  which  is  inclined  to  the  vertical 
axis  of  the  specimen.  The  area  of  this  shear  fracture  surface  AA 

s 

is  of  the  same  order  of  magnitude  as  AA  and  this  type  fracture  is 
essentially  mode  TI  ( ie  fracture  due  to  shear).  The  energy  associ¬ 
ated  with  this  fracture  can  be  expressed  as 


..II  ..IT  a* 

;  =  G,  AA  =  G,  AA 

cv  f  s  f  t 


By  equating  the  energy  in  ompressive  splitting  tests  to  that  in  com¬ 
pressive  shear  an  approximate  expression  can  he  determined  for  AAc 
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Figure  3.1  A  Comparison  of  softening  in  tension  and  compression 

and  therefore 


h 

c 


h 


t 


where  h  here  is  the  height  of  the  test  specimen.  From  test 

C  III 

results  it  was  determined  that  the  ratio  /G^  increases  rapidly 

with  confining  stress  and  the  following  empirical  formula  was 

developed . 
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where  As  and  Bs  are  determined  (in  a  manner  similar  to  that  dis¬ 
cussed  for  the  hardening  coefficients  A^,  and  C^)  based  on  low 
and  high  confining  stress  tests. 

3.6.5  Non-Associated  Flow 


As  pointed  out  in  previous  discussions  constitutive  equations 
based  on  associated  flow  plasticity  which  do  not  incorporate  a  plas¬ 
tic  volume  change  tend  to  overpredict  lateral  strains  or  dilatancy. 
This  of  course  is  due  to  the  fact  that  in  general  these  models  make 
use  of  the  normality  principle  for  plastic  strain  direction  while  the 
effects  of  the  hydrostatic  component  of  the  stress  increment  is 
ignored.  To  reduce  this  predicted  dilatancy  a  plastic  potential 
function  not  equal  to  the  yield  surface  can  be  used.  The  main  prob¬ 
lem  here  is  that  this  produces  a  non-symmetric  material  stiffness 
matrix  which  leads  to  increased  computational  time.  It  seems  reason¬ 
able  to  assume  a  plastic  potential  function  that  is  similar  in  form 
to  experimentally  verified  yield  and  loading  surfaces  and  also  is 
consistent,  in  terms  of  predicted  lateral  strains,  with  test  results. 

For  the  FEBM  non->associated  flow  is  characterized  by  a  plastic 
potential  function  (Q)  which  is  not  equal  to  the  loading  surface  F. 
The  plastic  potential  for  non-associated  flow  hardening  response  is 
defined  as. 

2 

+  k2  p  -  k2  C  =  0  3.38 

c 


Q ( CJ" j  »  0  g  ,  h ,  M  )  — 


(1-k) 


(01  " 


°3) 


:,2 


f' 

C 


and  for  softening  response 


Q(o1 ,  o3,  k,  M  ) 


r 

C 


C  =  0 


3.39 
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From  equations  (3.38)  and  (3.39)  it  can  be  seen  that  the  only  differ¬ 
ences  in  the  plastic  potential  and  the  loading  surfaces  (i.e.  equa¬ 
tion  (3.30)  is  in  the  modification  of  the  friction  parameter  M  MQ. 
The  parameter  MQ  is  referred  to  as  the  dilatancy  function  since  it 
primarily  controls  the  amount  of  dilatancy.  For  non-associated  flow 
the  gradient  of  the  plastic  potential  defines  the  direction  of  the 
plastic  strain  increment  as 


e 

P 


X 


3Q 

3a 


where 


3Q 

3o 


where 


r  3Q_ 
30j 

3Q_ 

3a, 


3Q 

-c3' 


>  = 


1 

f  ’ 

c 
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h 


|-2(1  -  k)o1+  1]  +  k2  3Mq 

L  f’  J  3a, 


-  S. 


(1  -  k) 


f* 
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°1  -  °3 

f  ’ 
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From  the  above  it  can  be  seen  that  —  is  a  measure  of  the  lateral 
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component  of  the  plastic  strain  .vector  and  can  be  determined  from 
triaxial  compression  tests  at  different  confining  stress  levels. 
Essentially  MQ  is  determined  from  three  tests,  unconfined  tension, 
low  confined  compression  and  moderate  confined  compression.  An 
iterative  procedure  was  used  to  select  MQ  to  best  fit  the  experi¬ 
mental  data.  It  was  found  that  an  empirical  relation  for  MQ  as  a 
function  of  a^  could  be  determined  as 
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3.7  Endochronic  Plasticity  Model 
3.7 .1  Background 

The  first  endochronic  plasticity  theory  was  proposed  by 
Valanis  (19)  in  1971  and  was  based  on  the  principle  that  the  history 
of  deformation  of  a  material  could  be  defined  in  terms  of  a  "time 
scale"  which  is  not  real  time  but  is  a  property  of  the  material.  The 
time  scale  is  assumed  to  be  a  monotonically  increasing  function  of 
the  total  deformation  of  the  material  and  is  essentially  a  metric  of 
the  deformation  path  in  strain  space.  Also,  the  theory  does  not 
depend  on  the  notion  of  a  yield  surface  and  definitions  of  loading, 
unloading,  or  reloading  criteria  as  is  the  case  of  conventional  plas¬ 
ticity.  The  first  theory  attracted  considerable  review,  and  discus¬ 
sion  in  the  engineering  community.  Probably  the  most  prominent 
criticisms  of  the  theory  were  presented  by  Sandler  (20)  and 
Fivlin  (21).  Sandler's  criticism  was  based  on  the  fact  that  the 
theory  fails  to  predict  closed  hysteresis  loops  during  unloading  and 
reloading  of  a  simple  endochronic  material.  Sandler  then  showed 
where  this  problem  could  lead  to  numerical  instability  and  non¬ 
uniqueness  of  solution.  Rivlin  carefully  analyzed  the  feasibility 
and  plausibility  of  using  the  endochronic  theory  to  model  materials 
in  general.  Rivlin  was  also  critical  of  the  thermodynamic  arguments 
used  in  the  development  of  the  theory.  The  problem  of  hysteresis 
loop  closure  was  addressed  by  Valanis  (22)  and  formed  the  basis  of 
the  new  endochronic  theory  for  concrete  which  is  discussed  in  Sec¬ 
tion  3.6.3.  A  key  feature  of  the  new  endochronic  theory  is  that  the 
time  scale  (intrinsic  time)  is  assumed  to  be  a  monotonically 
increasing  function  of  the  plastic  strain  increment.  A  detailed  dis¬ 
cussion  of  the  derivation  of  the  new  endochronic  theory  is  presented 
by  Valanis  and  Read  (4). 
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3.7.2  The  First  Endochronic  Theory  (Simple  Model) 

Consider  a  one-dimensional  test  of  a  material  which  is  linear  in 
stress  versus  strain  and  whose  deformation  is  also  a  function  of 
time.  When  an  axial  stress,  which  is  a  function  of  time,  is  applied 
to  the  material,  the  strain  is  a  function  of  the  total  history  of 
loading  up  to  the  present  time.  If  the  program  of  loading,  here 
c(t)  ,  is  continuous  and  differentiable  then  in  a  small  interval  of 
time  dx  the  stress  changes  by  ^  dx  this  change  in  stress 
results  in  a  change  in  strain  de  and  the  relationship  between  these 
changes  in  stress  and  strain  during  this  short  interval  in  time  is 

de(t)  =  C(t-x)  ~(--l  dx  3.40 

at 

The  total  strain  at  time  (t)  is  then 

e(t)  =  J  C(t-x)  -1^-  dx  3.41 

o 


A  material  whose  response  can  be  reasonably  approximated  by  (3.41) 
can  be  referred  to  as  a  linear  hereditary  material.  If  in  reverse  we 
consider  a  program  of  straining  which  results  in  a  change  in  stress 
we  obtain  a  relationship 


o(t) 


fu'-oigi  dr 


3.42 


The  functions  c(t)  and  k(t)  are  referred  to  as  the  creep  function 
and  relaxation  function  respectively  and  are  properties  of  the  mate¬ 
rial.  C(t)  can  be  thought  of  as  the  strain  produced  by  a  suddenly 
applied  constant  stress  of  unit  magnitude,  while  k(t)  represents 
the  stress  required  to  produce  and  maintain  a  constant  strain  of  unit 
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magnitude.  Linear  hereditary  material  behavior  can  be  approximated 
using  a  Maxwell  model  which  is  presented  in  the  Figure  below. 


k  c 


Maxwell  Model 


The  change  in  strain  for  a  differential  change  in  time  for  the  model 
is 


de  =  dc  +  de. 

s  d 


and  since 


o  =  E  e  = 
s 


CEd 


Therefore 


e 

s 


3.43 


When  a  Maxwell  model  is  subjected  to  a  suddenly  applied  constant 
strain  eq  ,  and  the  corresponding  value  of  stress  is  Cq  equation 
(3.43)  can  be  integrated  subject  to  these  initial  conditions  to 
obtain 


o(t) 


-Et/c 
e 
o 


a 


E  c 

o 


-Et/c 

e 


3.44 
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This  defines  the  relaxation  function  for  a  Maxwell  model.  The  rate 
of  stress  change  is 

P ’)  • Et/C 

so  that  the  initial  rate  (i.e.  at  t  =  0+)  is 


The  relaxation  time  is  defined  as  the  time  for  the  stress  to  relax  to 

• 

zero  if  the  relaxation  rate  continued  at  (a(t  =  0  )  .  Under  these 
these  conditions  the  equation  for  stress  relaxation  is  easily  deter¬ 
mined  from  the  figure  below 


Stress  relaxation  in  Maxwell  Model 

Therefore  t  is  determined  as 
K 


108 


and  is  a  property  of  the  material.  Equation  (3.43)  can  be  rewritten 
as 


or 


de  _  do  o 

dt  Edt  Et_ 

K 


3.45 


de 


3.46 


If  now  we  wish  to  use  equation  (3.46)  to  model  static  loading  of  a 
rate  independent  material,  the  time  scale  is  unimportant  and  we  could 
define  another  measure  of  time,  (i.e.  intrinsic  time)  t,  such  that 
the  differential  of  intrinsic  time  could  be  written 


d£  =  |de 


Replacing  dt  in  (3.46)  with  dt, 

de  *  r  +  'de'  3-*7 

This  is  the  simplest  endochronic  model  and  is  the  form  of  the  first 
endochronic  model  proposed  by  Valanis.  From  (3.47)  it  is  seen  that 
for  loading  de  >  0  therefore 


do 

de 


de 


>  0  -  E  (l  -  o/EtR)  -  E, 


while  for  unloading  de  <  0  and 


do 

de 


de  <  0 


*  E  (l*  «rttR)  -  E„ 


so  that  for  o  >  0  and  large  enough  Et^ 

0  c  El  <  Eu 

This  formulation  leads  to  a  response  as  shown  below  and  clearly  vio¬ 
lates  Drucker's  stability  postulate  for  an  unload  reload  cycle  of 
stress.  This  problem  was  alleviated  by  Valanis  [22]  by  introducing  a 
new  intrinsic  time  measure  which  will  be  discussed  in  Section  3.6.3. 


e 


Unload-reload  for  first  endochronic  model 
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3.7.3  The  New  Endochronic  Theory 
3. 7. 3.1  Basic  Equations 

As  mentioned  previously  the  endochronic  theory  is  based  on  the 
hypothesis  that  the  current  state  of  stress  in  a  material  is  a  linear 
functional  of  the  entire  history  of  deformation  where  the  history  is 
defined  with  respect  to  an  ^trinsic  time  which  is  a  function  of  the 
plastic  strain  increment.  The  foundation  of  the  theory  is  based  on 
the  theory  of  irreversible  thermodynamics  of  internal  variables,  the 
derivation  of  some  of  the  equations  of  the  theory  will  be  presented 
in  this  section. 

For  the  development  of  the  small  deformation  isotropic  theory  it 
is  assumed  that  internal  dissipative  mechanisms  can  be  represented  by 
internal  variables  gr  (second  order  tensors,  representing  devia- 
torlc  mechanisms)  and  pr  (scalars,  representing  hydrostatic 
mechanisms).  Furthermore,  there  exists  a  free  energy  of  deformation 
function  i p  which  is  quadratic  in  terms  of  the  strain  tensor  and 

internal  variables.  Therefore  f  can  be  represented  as 

ill  =  *  + 

where 

+D  =  %  (§’  9r) 

=  pr) 

g 

The  function  can  be  further  decomposed  into  elastic  <p  and 
plastic  t|/P  parts  such  that 

=  i)>P  (eP,  gr) 

<p£  =  (eP’  pr^ 


The  thermodynamic  relations  for  the  plastic  parts  of  the  constitutive 
equations  are. 


Ill 


~  3eP 

3€P 

Due  to  the  assumptions  made  for  (ie  isotropy,  quadratic  in  strain 
and  interna1  variables)  it  can  be  reduced  to  the  canonical  forms 

=  1-2  u  !  |ep  -  gr|  |  2 

2  r 

and  3.48 

*p »  I  2.  Kr  <eP  -  pr)“ 

n  2  r  r 

A  simple  one  dimensional  model  of  an  endochronic  material  is  pre¬ 
sented  on  the  following  page. 

The  force  in  each  spring  is 

Qr  =  pr  (eP  -  qr) 

and  the  contribution  to  the  free  energy  of  endochronic  element 
r  is 


. r  1  ,  p  ,2 

^  ur  (eK  -  gr) 

and  therefore 


One-dimensional  model  with  several  endochronic  elements 


also 


^  pr  (eP  -  gr)  =  gr  =  S 

3eP  ~  r  *  - 

and 

"  Ifr  =  vr  (ep  -  gr)  -  gr 

Similar  equations  can  be  written  for  the  hydrostatic  part.  Qr  and 
Pr  are  internal  stress  type  quantities  and  gr  and  pr  internal 
strain  type  quantities  with  the  characteristics  that  they  cannot  be 
measured  as  we  normally  measure  stress  or  strain.  Since  the  entropy 
for  the  system  must  be  increasing 
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dqr  >  0  |  jdqr  |  |  >  0 


-  3^7  dpr  >  0  |dpr  |  >  0 


and  therefore  expressing  dqr  ,  and  dpr  ,  in  rate  forms  we  have. 

-  Qr  •  qr  >0 

-  Pr  •  pr  >0 

For  irreversible  thermodynamics  it  is  usually  assumed  that  consti¬ 
tutive  equations  give  the  forces  (stresses)  as  functions  of  the 
fluxes  (strains)  or  vice  versa.  This  assumption  can  be  expressed  as 


(1)  • 

Qr  =  hr  '  qr 


Pr  =  br(2)  pr 


No  sum  on  r 


In  these  expressions  br^  and  br^  play  the  same  role  that  the 
scalar  shear  and  bulk  moduli  play  in  relating  elastic  deviatoric 
stresses  to  deviatoric  strains  and  hydrostatic  stresses  to  volumetric 
strains . 

From  the  previous  expressions  for  Qr  and  Pr  the  "evolution"  equa¬ 
tions  are  developed  as 


3^ 

+  br(1)  qr  =  0 

3ar 

+  br(2)  pr  =  0 

3pr 
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Finally,  recalling  Equation  3.48. 


and 


or 


pr  <jr  +  br^  gr  =  pr  e ^ 


kr  pr  +  br^^  pr  =  kr 


br 


dQr 

^fry  2r  +  5r 


=  pr 


and 


IL-  Pr  +  4EI 

(2)  lr  dt 


Kt 


br 


del 

dT 


de* 

dt 


For  a  rate  independent  material  the  independent  variable  in  the  above 
equation  does  not  have  to  be  real  time,  but  must  only  be  a  monotoni- 
cally  increasing  parameter  which  effectively  records  the  history  of 
loading  of  the  material. 


For  the  endochronic  theory  a  parameter  is  chosen  and  defined  as 
intrinsic  time  (Z)  which  is  a  function  of  plastic  strain. 

Replacing  dt  with  dz  the  evolution  equations  become 


and 


pr 


br 


(1) 


Qr  + 


dQr 

dz 


pr 


d£ 

dz 


U  ■  Pr  +  *21 

br«> 


=  Kr 


dt* 

dz. 


These  equations  ran  be  solved  uslnv  T.aplace  Transform  techniques  with 
Z  as  the  independent  variable.  This  procedure  for  the  deviatoric 
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evolution  equation  is  presented  here.  First  the  Laplace  transform  of 
each  term  is  taken. 


{  r  1 

dQr 

V  deP) 

1  pr 

>  Qr  ,+  L  < 

.  ss  T  / 

~  1 

|br(1) 

dz 

l  SJ 

—  Lt  \ 

dz  | 

l  S  J 

Therefore  defining  ar  = 


yr 


br 


(1) 


and 


ar  Qr  (s)  +  s  Qr(s)  -  Qr(o)  =  yrs  e^(s)  -  preF(o) 


eH(o)  =  0 


Which  can  be  rearranged  to  express  Qr(s)  in  terms  of  ep(s) 


Qr  (o) 

=  ?C-r"'L T  +  Ur 


1 


(a  +  s)  Ca  +  s) 


s  ep(s) 


Taking  the  inverse  transform  of  this  equation,  and  observing  that  the 
second  term  is  the  product  of  two  transforms  and  therefore  its 
inverse  is  a  convolution  integral 


z  ,  deK 

„  _  ,  .  -ar  z  r  s  -ar  {z  -  z  )  ,  , 

Qr  =  Qr(o)  exp  s  +  ur  J  exp  s  s  -j^r  dz 

o 


and  similarly  for  the  hydrostatic  part. 


_  ,  .  -Ar  z  / H  -Ar  (z  -  z')  deP  ,  , 

Fr  =  Pr(o)  e  H  +  Kr  J  e  H  -—f  dz 

o 

Kr 


br 


(2) 


where 


Ar 
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Finally  since 


1  Or  and  o  =  1  Pr 

r  4  r 

z  deP 

“  Qr<°)  e,r  7's  +  /  p(zs  '  zs)  "dz*  dZ’ 

o 

and  3.  A  8a 

o  =  Pr(o)e  rZs  +  ®(zu  -  z')  —~t  d  z  ' 

r  I  H  d  z 

o 

where 

-<*r  z 

P  =  Ur  e  s 

r 

$  =  Xr  iXr  ZM 

r 


Kquations  (3.Afia)  are  the  general  constitutive  equations  for  the 
endochronic  theory.  It  is  important  to  note  here  that  the  develop¬ 
ment  of  the  endochronic  constitutive  equations  follows  accepted 
conventional  procedures  and  the  concept  of  intrinsic  time  or  pseudo 
time  is  not  new  (e.g.  viscoelasticity  applications).  Furthermore, 
the  exponential  forms  in  3.48a  result  from  direct  application  of 
applied  mathematics  principles  and  not  from  arbitrary  assumptions  or 
curve  fits  of  experimental  data.  The  basic  equations  of  the  new 
endochronic  plasticity  theory  can  be  expressed  as 


deP 


S  =  f'  P  (zg  -  7. '  ) 

dz' 

dz’ 

3. 49 

0  =  J~7\\  -  7  ’  ) 

dcp 
dz ' 

dz' 

3.50 

ds 

dep  =  de 

“  2U 

3.51 

dfp  =  dc 

do 

k 

3.52 
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dz 

dz 


2  =  !!deP]'2 


dz 

F 


+  k2  |deP  | 2 


dz 


H 


dz 

kF 


H 


3.53 

3.56 


where 


S 


y 

K 

F 

s 

FH 

z 


z 

s 

ZH 

k 


p(z)  ,  *(z) 


Deviatoric  stress  tensor. 

Hydrostatic  stress. 

Denotes  plastic  components. 

Elastic  shear  modulus. 

Elastic  bulk  modulus. 

Shear  hardening  parameter. 

Hydrostatic  hardening  parameter. 

Intrinsic  time  scale. 

Intrinsic  time  for  shear  response. 

Intrinsic  time  for  hydrostatic  response. 

A  constant  which  determines  the  magnitude  of  shear 
volumetric  coupling 

Weakly  singular  kernel  functions  such  that 
p(0)  =  4>(0)  =  » 

but  p(z)  and  <J>(z)  are  integrable 
in  the  domain  0  ^  Z  <  00 


Intrinsic  time  is  analogous  to  recorded  history  parameters  used  in 
conventional  plasticity.  The  form  of  dz  in  Equation  (3.53)  was 
developed  based  on  the  idea  that  the  intrinsic  time  increment  is 
given  by  the  increment  of  the  plastic  strain  path  in  plastic  strain 
space  and  therefore  could  be  written  in  general  as 


l  p  .  p 

dz  =  P  ..  ..  del 

ijki  ij  kl 


3 . 6 


where  the  metric  P  is  a  fourth  order,  isotropic  tensor  with  th 
general  representation 


Fijkl 


ij  kl 


+  K„ 


ik 


3.58 


Combining  3.57  and  3.58 


dz  kl  de?t  dejk  +  k2  dej\  de^ 
or 

dz2  =  k2  jdeP|2  +  k2  ||dep|r  3.59 

which  is  equivalent  to  Equation  (3.53)  to  within  an  immaterial  multi¬ 
plicative  constant.  The  requirement  that  p(z)  and  0(z)  in  Equa¬ 
tions  3.49  and  3.50  be  weakly  singular  kernels,  essentially  comes 
from  the  need  for  the  new  model  to  predict  closed  hysteresis  loops 
under  unloading  -  reloading  cycles  (this  was  a  major  problem  with  the 
first  model) .  This  idea  is  demonstrated  in  the  following  example. 
Consider  a  typical  plot  of  simple  shear  stress  S  versus  shear 
strain  y  as  shown  in  Figure  3.15.  If  is  assumed  in  this  figure 
that  a  yield  point  exists  but  that  the  slope  of  the  stress  strain 
curve  is  continuous  throughout.  If  we  wish  to  plot  shear  stress 
versus  plastic  shear  strain  y  ,  as  shown  in  Figure  3.15,  the  S,y 

axes  can  be  rotated  as  shown  and  the  slope  of  the  curve  at  y  =  0 

P 

will  be  infinite. 

The  endochronic  theory  for  this  case  will  yield 


■ 


Therefore 


ds 

dy1 


=  p(z  -  z’)  = 
s 

yp  -  o 


which  requires  p(0)  =  . 
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The  series  expansion  for  p 


P 


-ar 
yr  e 


z 

s 


satisfies  the  requirement  p(0)  ■  00  . 


The  parameters  which  must  be  calibrated  in  the  theory  can  be 
divided  into  three  groups:  parameters  which  describe  pure  hydro¬ 
static  response,  parameters  which  describe  the  shear  volumetric 
coupling,  and  those  parameters  which  describe  the  deviatoric  (shear) 
response.  The  calibration  of  these  parameters  will  be  discussed  in 
the  following  sections. 

3. 7.3.2  Hydrostatic  Response 

Typical  response  features  of  concrete  to  pure  hydrostatic  stress 

were  discussed  in  Chapter  2.  The  endochronic  functions  $(Z^)  and 

F  must  be  defined  and  calibrated  to  capture  these  key  hydrostatic 
H 
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stress-volumetric  strain  features  of  the  concrete  in  the  stress  range 
of  interest  for  a  particular  problem.  Assuming  that  the  reference 
state  is  the  natural  state  then  under  monotonic  pure  hydrostatic 
loading  conditions  the  following  basic  equations  apply: 


Figure  3.16  Hydrostatic  stress-strain  response  feature 


-V  •  Be 

then  expressing  e 

linear  terms  we  have 


P 


in  a  power  series  expansion  and  only  retaining 


FR  =  1  +  Bep 


From  Equations  3.53  and  3.56  it  is  easily  shown  that  (since  for 
hydrostatic  response  deP  =  dzg  =  0)  . 


d2H 


dep 

1  +  Bep 


which  upon  integration  leaves 

In  (I  +  Bep) 
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therefore 

eBzH  =  1  +  Bep  =  F 

n 

Also  it  is  noticed  that  pure  hydrostatic  response  is  independent  of 
k.  Therefore  the  constitutive  equation  for  hydrostatic  response 
becomes 


0  = 


?  / 


'H  -Ar  (z  -  z')  Bz'  , 
Kr  e  R  e  dz 


or  upon  integration 


a  = 


?  :>t  (l TTl) 


e(Ar  +»)«„.! 


Where  r  represents  the  number  of  internal  variables  K^  ,  X^  (or 
the  number  of  endochronic  elements  in  parallel)  required  to  simulate 
the  hydrostatic  response  to  the  desired  level  of  accuracy.  Therefore 
there  are  three  parameters  ,  A^  ,  anc*  &  which  must  be  deter¬ 

mined  based  on  pure  hydrostatic  test  results. 

The  calibration  of  the  parameters  ,  Af  ,  and  B  is  accom¬ 

plished  by  first  defining  a  stress  a*  such  that 


a  _  o 

FH  1  +  BeP 


and  since  for  the  hydrostatic  test 

e  Bep 

Be  =  e  we  have 


a*  = 


l 

r 


Kr 

(B  +  Ar) 


BeP  =  Bz^-  and  if  we  assume 


1  + 


1 


-(B  +  Ar) 
e 


so  the  two  curves  of  Figure  3.16  can  now  be  used.  B  is  determined 
from  the  straight  line  portion  of  the  second  branch  of  the  hydrostat 
from  the  requirement  that 


Next  a  plot  is  made  of  (a  -  a*)  versus  Z  (here  Z  =  -x  In  (1  + 
6e^))  based  on  hydrostatic  data  (note  B  is  known)  and  values  of 
and  are  selected  which  give  the  desired  fit  to  this  curve. 


1 24 


3. 7.3.3  Shear  Response. 

The  shear  response  is  defined  by  Equations  3. 49,  3.53,  and  3.56 

which  require  determination  of  the  response  functions  p(Z  )  and  F 

s  s 

along  with  the  coupling  constant  k.  As  discussed  in  Chapter  2  the 
hardening  characteristics  of  concrete  in  shear  are  very  sensitive  to 
the  confining  stress  level,  also  this  hardening  is  dependent  on  the 
direction  of  loading  in  the  deviatoric  plane  which  can  be  defined  by 
the  similarity  angle  Q  .  Therefore  it  is  assumed  that 

F  =  F  (o,0) 
s  s 

and  Fg  can  be  determined  from  a  series  of  hydrostatic  pure  devia¬ 
toric  tests  as  discussed  in  Chapter  2.  The  simplest  and  most  direct 
way  to  determine  the  shear  parameters  is  from  a  triaxial  test  in  the 
tt  plane  (ie  shear  under  zero  hydrostatic  stress).  However,  there 
were  no  tests  performed  like  this  for  the  concretes  reported  herein 
and  therefore  hydrostatic  pure  deviatoric  test  results  must  be  used. 
Changes  in  hydrostatic  plastic  strain  eP  due  to  increments  of  shear 
stress  result  due  to  change  in  Z  as  defined  in  Equation  (3.53)  and 
the  corresponding  change  in  Z  from  Equation  3.56.  If  the  kernel 
$(ZH)  in  Equation  3.50  is  approximated  with  a  Dirac  6  function 

*<V  -  K„  S(V 

Equation  3.50  results  in 


and  for  pure  hydrostatic  loading 

a  =  K  (1  +  gEP)  3.62 

o 

Therefore  Equation  (3.62)  is  essentially  a  linear  hardening  model. 
Equation  3.62  seems  reasonable  especially  for  the  initial  portion  of 
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the  hydrostat  (le  prior  to  signs  of  locking  up)  and  if  is  the 

mean  stress  during  a  pure  deviatoric  load  increment  we  can  write 

a. 

deP  =  —  dzR  3.63 

o 

Therefore  as  dZ  increases  due  to  application  of  pure  shear  incre- 
ments,  due  to  equations  3.53  and  3.56,  dcp  will  increase  as  defined 
by  Equation  3.63.  Defining  the  intrinsic  time  dZ  in  terms  of  a 
parameter  d£  such  that 


d4s  =  ! !dePl I  *  d?H  =  ldep! 

a  differential  relationship  between  d£g  and  dy  can  be  derived, 
where  dy  is  the  change  in  the  intrinsic  time  in  addition  to  the 
change  caused  by  pure  hydrostatic  loading.  This  expression  is 


and 


where 


_ 1 

1  +  2ay  | 


k  3 
o 


3.64 


y 


Z-Z^,  7}  is  the  value  of  Z  at  the  completion 
of  the  pure  hydrostatic  branch  of  loading. 


Equation  3.64  can  be  integrated  to  give 


4 

s 


1_ 

2a 


( *  +  yjx2  -  lj 


3.65 
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K  6 

OSS 

where  X  =  1  +  — —  is  the  hydrostatic  strain  caused  by 

shear  loading  at  constant  pressure. 


Finally  if  W  is  defined  as  the  increase  in  Zg  after  the 
hydrostatic  loading  such  that. 

w  =  z  -  z1 

s  s 


a  differential  relation  between  d^  and  dW  can  be  determined  as 

s 

d*s 

— -  =  p 

dW  s 

which  governs  dCg  during  loading  in  a  deviatoric  section  after 
hydrostatic  loading.  Since 


2aF  W 
s 


1  +  2aF  W 
s 


'3.66 


F  =  F  (o,0) 
s  s 

Fg  will  remain  constant  during  loading  in  a  deviatoric  section  at 
constant  0.  Finally  for  this  type  of  loading  Equation  (3.49)  can  be 
written  for  a  pure  shear  path  as 


and 


=  f  p(W  -  W1)  dW1 

o  dW 


=  FS  fW  ,  1.  2aFsW  1  .  1 

'  /o  J  (  '  I_1  +  2aFgwJ 


3.67 
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and  the  form  of  Fg  can  be  determined  by  comparing  t®  at  different 
mean  stress  levels. 

t°°  (a  ) 

- !_  =  F 

-«C  f  -  \  S 

1  v^b' 

The  functional  form  of  p  is  determined  by  solving  Equation  (3.67) 
numerically  knowing  t  and  a  (from  experiments) . 
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Chapter  4 

Comparison  of  Model  Predictions  Versus  Test  Results 

4 . 1  General 

The  Fracture  Energy  Based  Model  (FEBM)  and  the  Endochronic  Con¬ 
crete  Plasticity  Model  (ECPM)  are  calibrated  for  the  concrete  materi¬ 
als  used  in  the  verification  tests,  and  stress  strain  response 
predictions  are  discussed  and  compared  with  test  results  in  this 
chapter.  The  concept  and  operational  details  of  the  WES  Constitutive 
Driver  (WCD) ,  which  is  used  in  this  study,  is  discussed  in  Sec¬ 
tion  4.2.  Calibration  of  the  constitutive  models  and  an  analysis  of 
predicted  and  measured  results  are  presented  in  section  4.3  WES 
Tests,  4.4  Colorado  Tests  and  4.5  Eindhoven  Tests.  Most  emphasis  is 

placed  on  comparisons  for  the  WES  Tests,  because  a  more  extensive 

1 

calibration  test  program  was  conducted  at  WES  for  the  f^  =  6.5  ksi 
concrete.  In  essence  a  calibration  parameter  sensitivity  study  is 
presented  in  this  chapter.  The  objective  of  this  sensitivity  study 
is  not  to  optimize  the  values  of  calibration  parameters  but  rather  to 
discuss  the  consistency  between  parameter  values  determined  from  a 
calibration  tests  and  those  values  of  the  parameters  required  for 
better  model  predictions. 

4.2  WES-Constitutive  Driver 

The  original  concept  for  the  WES  Constitutive  Driver  (WESCD)  was 
based  on  the  need  to  develop  a  consistent  method  to  better  evaluate 
predictions  from  large  scale  dynamic  finite-element  codes.  As  dis¬ 
cussed  in  section  1.1,  the  effectiveness  and  efficiency  o'  the  con¬ 
stitutive  model  was  obscured  by  the  overall  size  and  complexity  of 
these  large  scale  problems.  Therefore,  the  concept  was  to  develop 
and  implement  a  few  rational  and  consistent  constitutive  models  in  an 
easy  to  use  modular  software  package,  with  graphics,  and  a  signifi¬ 
cant  experimental  data  base.  This  would  allow  the  objective  evalua¬ 
tion  of  different  models  and  also  provide  an  assessment  method  which 
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could  be  used  to  Improve  the  predictive  capabilities  of  specific 
models.  The  WT.SCD  consists  of  the  following  modules: 

1.  User  defined  loading  histories  for  driving  constitutive 
models  in  principal  stress,  or  strain  space  or  under  mixed 
control . 

2.  Constitutive  models. 

3.  Solution  models  for  numerical  integration  of  the  constitu¬ 
tive  equations. 

4.  Experimental  data  base. 

The  WESCD  is  written  in  Fortran  77  with  graphic  components  based  on 
Hewlett  Packard-Advanced  Graphics  Package.  The  user  interface  for 
the  WESCD  is  described  in  the  following  4  "pages". 


Page  1 :  Selection  of  Concrete  Model 

The  first  page,  shown  in  Figure  4.1  allows  the  user  to  choos. 
the  constitutive  model  of  interest. 


Figure  4.1  Constitutive  model  selection 


Selection  of  Material  Parameters 


Page  2: 

The  second  page,  shown  in  Figure  4.2,  displays  some  of  the  char¬ 
acteristic  features  of  the  constitutive  model  of  interest.  In  Fig¬ 
ure  4.2  deviatoric  and  meridian  sections  of  the  loading  and  failure 


MERIDIAN  PLANE 


DEUIATORIC  SECTIONS 


,  c-i 


DEGRADATION  OF  FA  I  LU  RESURFACE 


INPUT  FILE  NAME  > 


Material  parameter 

E  .630E+07 
POISSON  .230E+00 
F'  c  .660E+04 


WES-U I CKSBURG,  MS 
CONST.  DP.IUER 
Uers .  002-1988 


PLASTICITV 

Leon 


INPUT  FILE 


F'  t 

RESIDUAL 
DI  SP 
CRCK  SP 
INIT  HARD 
PARAMETER 
PRIN  STRS 
RATIO  AT 
TP 


.  600E+03 
.  100E-01 

.  425E+01 
.  100E+00 

.  167E+00 


UNITS 


|  US  ( Lbs. -In.) 


SI  <N.  -  MM.) 


REPEAT]  NEXT  |  EXIT 


Figure  4,2  Constitutive  model  parameters 


surfaces  for  the  FEBM  are  presented  with  indications  where  the  sur¬ 
faces  degrade  during  softening.  Also  on  this  page  the  user  is 
prompted  to  specify  certain  material  parameters  for  the  concrete.  To 
specify  or  vary  the  internal  parameters  (e.g.  hardening  parameters 
for  the  FEBM)  the  user  must  make  changes  within  the  constitutive 
model  module. 

Page  3:  Selection  of  Load  History 

This  page,  as  shown  in  Figure  4.3,  prompts  the  user  to  select 
either  stress,  strain,  or  mixed  control.  The  primary  reason  for  the 
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development  of  a  mixed  control  formulation  was  for  simulation  of 
tests  such  as  uniaxial  strain. 


Figure  4.3  Control  mode 

User  input  histories  can  be  defined  or  a  history  in  the  Experimental 
Data  Base  Module  can  be  selected.  A  variety  of  plotting  options  are 
available  for  presenting  experimental  and  predicted  results. 

Page  4:  Solution  of  Response  History  Analysis 

This  page  Figure  4.4  presents  the  results  of  the  analysis  for 
the  particular  plotting  option  selected.  For  the  example  presented 
in  Figure  4.4,  axial  stress  versus  axial  and  circumferential  strain 
are  presented  for  the  FEBM  compared  to  experimental  results.  A 
rescaling  option  is  available  to  force  identical  scales  on  different 
plots  for  comparison  purposes. 
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4.3  Comparison  of  Model  Predictions  with  WES  Tests 

4.3.1  General 

Verification  tests  conducted  at  WES  include:  VT6.5-1  through 
VT6.5-5  and  VT2-1.  These  tests  were  briefly  discussed  in  Section  2.7 
and  will  be  discussed  in  more  detail  here  along  with  comparisons  with 
predicted  results  from  the  FEBM  and  the  ECPM. 


4.3.2  The  FEBM  Test  VT6.5-1 

Test  VT6.5-1  provided  excellent  data  for  model  verification,  as 
all  four  strain  gages  were  effective  throughout  the  entire  test.  The 
as  measured  axial  stress  versus  axial  and  circumferential  strains  are 
presented  in  Figure  4.5.  From  pretest  and  posttest  micrometer  mea¬ 
surements  the  permanent  axial  strain  was  0.0133  in/in.  and  the  perma¬ 
nent  circumferential  strain  was  0.0067  in/in.  Based  on  these 
observations  the  axial  and  circumferential  strain  gages  indicated  in 
Figure  4.5  were  selected  for  use  in  model  comparisons.  The  differ¬ 
ence  in  measured  strains  in  Figure  4.5  is  due  to  normal  variation  in 
electrical  measurements,  slight  differences  in  bonding  for  the  dif¬ 
ferent  gages,  and  probably  most  due  to  the  fact  that  the  different 
1/2-in.  gages  are  bonded  to  different  percentages  of  aggregate  and 
cement  paste  at  the  different  locations.  The  effect  of  confining 
stresses  up  to  approximately  40  ksi  on  the  response  of  these  gages 
has  been  evaluated  through  calibration  tests  on  steel  cylinders  and 
found  to  be  negligible. 

Based  on  the  calibration  tests  for  the  f’  =  6.5  ksi  concrete  the 

c 

model  parameters  for  the  FEBM  were  determined  as: 


Elastic  Modulus 
Poisson's  ratio 
Compressive  Strength 
Tensile  Strength 
Rupture  Displacement 


E  =  6,300  ksi 
o 

V  =  0.23 

f'  =  6.6  ksi 
c 

f  ^  =  0.6  ksi 

U  =  0.01  in. 
r 


RXlRL  S  "RES 
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Crack  Spacing 

Initial  Hardening  Parameter  kQ  =  0.1 


S*»ECI*fJ  B.  *-2-Z 
V5SDMFS  22:  30.S'53S  .  ? 
9  '2  I  '9> 


STPRIN  (  i  n  .  /  1  n  )  . - GRGE  6.  CIRCU*. 

-  GRGE  5 .  Rxjq, 

- grge  4 .  c  i ©cum 

-  GRGE  3.  RXIR^ 


Figure  4.5  Actual  test  data  for  VT6.5-1 


Hardening  coefficients  for  calculation  of  ductility 
parameter  x^ 

Ah=  5.9 
Bh  -  13.0 
C„  *  -0.6 

H 

Coefficients  for  calculation  of  friction  parameter  (Mq)  for 
non-associated  flow. 


D  =  12.0 
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E  =  20.0  +  10  exp  (-25 ( 1 — k) ) 

F  =  -2.0  +  3.0  exp  (-25 (1-k) ) 

The  coefficients  E  and  F  which  help  define  M  were  found  to  be 

q 

dependent  on  k  ,  and  the  exponential  term  was  used  simply  to  smooth 
the  transition  from  loading  surfaces  to  the  maximum  strength  surface 
k  =  I. 

Using  the  above  calibration  an  initial  prediction  for  test 
VT6.5-I  was  made,  and  the  results  are  presented  in  Figure  4.6  for 
strain  control,  associated  and  non-associated  flow.  As  discussed 
previously  the  associated  flow  rule  forces  the  plastic  strain  incre¬ 
ment  vector  normal  to  the  loading  surface  which,  in  the  loading 
region  of  this  test,  has  a  large  component  in  the  direction  of  radial 
expansion.  This  effect,  in  a  strain  driven  format,  is  manifested  by 
overpredicting  the  confining  stress.  The  non-associated  flow  predic¬ 
tion  compares  much  better  with  test  data  as  shown  in  Figure  4.6b  and 
this  version  will  be  used  in  the  following.  Three  values  of  elastic 
modulus  were  used  to  determine  prediction  sensitivity.  These  values 

included  the  initial  calibration  value  of  E  =  6.3  x  10^  psi  plus 
6  6  ° 

E  =  5.0  x  10  psi  and  E  =  4  x  10  psi.  The  results  of  these  runs  are 
presented  in  Figure  4.7.  From  these  results  it  is  seen  that  the  5  x 
10^  psi  modulus  provides  better  overall  comparison  with  the  data, 
especially  in  the  lower  (elastic)  region  of  response.  Also,  this  is 
not  an  unrealistically  low  value  for  E  ,  for  this  concrete  when  one 
considers  the  uncertainties  or  errors  in  measuring  E  .  At  this 
point  the  best  prediction  for  VT6.5-1  is  the  non-associated  flow  with 
E  =  5  x  10^  psi  as  shown  in  Figure  4.7b.  At  about  o  =  15  ksi  the 
model  begins  to  predict  to  stiff  a  response.  Values  of  at 

transitions  from  pure  hydrostatic  to  pure  deviatoric  loading  branches 
are  indicated  in  Figure  4.7b.  Although  each  loading  branch  is  itself 
made  up  of  smaller  stair  steps  the  main  branches  are  pure  hydrostatic 
(0-1,  2-3,  4-5),  which  primarily  compact  and  stiffen  the  material, 
and  pure  deviatoric  (1-2,  3-4,  5-6)  which  shear  the  material  and  also 


(E+4  psl)  a  (E+4  psl) 


36 


cause  volume  change  (through  shear  volumetric  coupling) .At  ' 

15  ksi,  the  next  loading  increment  causes  relatively  large  increments 
of  radial  extension  strains  and  nearly  plastic  flow  in  the  axial 
direction.  The  model  captures  the  shape  of  this  stress  strain 
response  but  is  too  stiff.  Referring  back  to  the  load  path  for  this 
test  (Figure  2.16a)  the  stress  point  is  seen  to  still  be  inside  the 
maximum  strength  surface,  however  considerable  plastic  loading  has 
occurred  prior  to  this  point  as  can  be  seen  (Figure  4.8)  referring  to 
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Chapter  2.  The  final  loading  branch  (5-6)  in  general  shows  a 
decreased  average  slope  in  for  the  test  (ie  compared  to 

branches  1-2  and  3-4)  up  until  the  maximum  axial  stress  of 
* — *  •  ^  pi-oi.  xa  icucucu  at  wuxtu  pc  lli  essentially  plastic  flow  occurs* 
The  model  does  not  show  this  decrease  in  average  e slope  but 
does  predict  plastic  flow.  The  maximum  stress  state  attained  in  the 
test  is  also  compared  with  the  maximum  strength  surface  in  Fig¬ 
ure  4.8.  This  comparison  is  excellent  indicating  little  or  no  path 
dependency  for  the  maximum  strength  surface.  Also  the  plastic  flow 
indicated  by  the  model  and  test  results  is  consistent  with  the 
response  expected  near  the  brittle  ductile  transition  point,  which 
for  this  concrete  was  predicted  to  be  near  4  ksi  (see  Section  2.6.2). 
Also  plotted  in  Figure  4.8  are  the  plastic  strain  increment  vectors 
near  the  =  15  ksi  anti  =  23  ksi  stress  point.  It  is  evident 
that  the  plastic  strain  increment  is  not  directed  normal  to  the  load¬ 
ing  surfaces  in  these  regions.  For  a  test  such  as  VT6.5-1  where  no 
strain  softening  occurs  the  FEBM  predicts  a  more  ductile  response  if 
the  ductility  parameter  X^  (Equation  3.32)  is  increased.  To  see 
this  effect  the  coefficients  in  Equation  3.32  were  increased  by  10% 
and  15%  with  no  significant  effect  on  predicted  results  as  shown  in 
Figure  4.9(a).  However,  an  increase  in  of  150%  resulted  in  the 

response  shown  in  Figure  4.9b.  The  point  to  be  made  here  is  that 
small  reasonable  variations  in  X  do  not  yield  significant  improve¬ 
ment  in  the  model  predictions,  for  this  test,  and  very  large 
increases  in  ate  unacceptable  since  this  would  indicate  th< 

model  might  need  recalibrating  for  every  test. 

The  other  model  parameter  which  can  significantly  affect  ductil¬ 
ity  and  the  control  of  radial  strains  for  the  non-associated  flow 
rule  is  the  frictional  parameter  M  .  To  see  this  effect  the  value 
of  M  was  varied  from  0.5  to  5  times  the  calibrated  value.  The 

q 

best  results  from  these  predictions  are  shown  in  Figure  4.10.  The 
spike  In  the  stress  strain  results  of  Figure  4.10(a)  is  probably  due 
to  some  relative  scatter  in  experimental  data.  Also,  this  point  is 


near  the  BDT  point  and  the  model  may  be  more  sensitive  to  perturba¬ 
tion  in  input  strains  due  to  the  large  value  of  M  used  here. 

Q 

Essentially  the  same  comments  apply  for  as  were  made  for  in 

that  large  changes  of  M  are  required  to  produce  questionable 
improvements  in  the  predictions. 

Thus  far  comparisons  of  the  FEBM  predictions  and  test  results 
from  VT6.5-1  have  been  based  on  axial  stress  versus  axial  and  circum¬ 
ferential  strains.  Plots  of  octahedral  normal  stress  (c  )  versus 

o 

octahedral  normal  strain  (e  )  can  be  used  to  provide  better 

o 

detailed  insight  into  the  volumetric  response  features  of  the  mate¬ 
rial  and  those  response  features  predicted  by  the  model.  Here  a 
decrease  in  eq  (ie  more  negative)  will  be  referred  to  as  compaction 
while  an  increase  in  (i.e.  less  negative)  will  be  referred  to  as 

expansion.  In  Figure  4. II  octahedral  normal  stress  0  (mean  normal 

o 

stress)  is  plotted  against  octahedral  normal  strain  (average 

strain)  for  the  model,  initial  calibration  versus  test  results. 


Figure  4.11  FEBM-VT6 .5-1 ,  comparison  of  predicted  and  measured 
octahedral  stresses  and  strains 
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Points  of  transition  from  pure  hydrostatic  to  pure  deviatoric  loading 
are  labeled  1  through  6  while  corresponding  points  (in  stress  space) 

for  unloading  are  labeled  6^  through  1*.  The  change  in  e  between 

1  0 
6  and  6  is  due  to  plastic  flow.  During  unloading  along  pure  shear 

branches  there  is  still  significant  plastic  compaction  (i.e.  branches 
6^-5^,  4^-3^,  2^-1*).  This  behavior  is  often  observed  in  soils. 

This  behavior  can  be  explained  by  noting  that  at  the  end  of  each 
hydrostatic  branch  (loading  or  unloading)  the  volume  of  the  material 
is  not  necessarily  at  its  relative  minimum  value  because  voids  can 
still  exist  and  localized  microcracking  can  occur  during  the  hydro¬ 
static  branch  due  to  the  actual  heterogeneity  of  the  material.  When 
shear  is  aDolied  (loading  or  unloading)  the  tendency  for  the  material 
is  to  close  those  voids  and  microcracks  which  results  in  further 
reduction  in  volume.  It  is  very  important  to  note  that  the  isotropic 
hardening  model  fails  to  capture  this  response.  This  key  response 
feature  points  out  the  need  for  some  type  of  kinematic  hardening  (ie 
where  loading  surfaces  can  translate  as  a  rigid  body)  or  other  model 
response  feature  to  capture  this  unloading  behavior. 
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4.3.3  EC PM  Versus  Test  V76.5-1 

The  ECPM  was  calibrated  for  the  f  =  6.5  ksi  concrect  based  on 

c 

the  general  procedures  for  model  calibration  as  discussed  in  sec¬ 
tion  3.6.3.  The  ECPM  developed  here  could  be  considered  an  elemen¬ 
tary  or  simple  model  since  only  one  endochronic  element  is  used  (i.e. 
one  kernel  function  for  shear  and  one  for  hydrostatic  response). 

These  kernels  are  represented  here  as 


-a.  Z 

o(Z  )  =  A  e  S 
s  1 


6  j  ZH 

0(ZH)  =  Pj  e 


Furthermore  the  modc.1  includes  a  coupling  term  so  that  the  hydro¬ 
static  response  is  given  by 


and  the  kernel  F  is  approximated  as 


-Co 

r(zH)  -  r  e  °  E  0(zH) 

The  coefficients  in  these  kernel  functions  were  determined  for  the 

f'  =  6.5  ksi  concrete  as. 
c 


Aj  =  5,279  ksi 
aj  =  140 
B  =  13,010  ksi 
Pj  =  1,301 
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To  =  0.176  ksi 
Co  =  19.9 

Values  of  the  shear  modulus  (G)  and  Initial  bulk  modulus  were  deter¬ 
mined  to  be . 


G  =  2,561  ksi 
Ko  =  3,000  ksi 

The  bulk  modulus  was  assumed  to  vary  linearly  with  the  hydrostatic 
stress  as 

K  =  Ko  +  Kj  o 
where  K^  =  25 

The  shear  hardening  parameter  Fs  and  the  hydrostatic  hardening 
parameter  are  approximated  as 

Fs  =  is  +  Bs  o 


and 


F 


H 


€ 


P 


where  is 


Bs 


e 


H 


1.63  ksi 

0.689 

70 


Fs  is  continuously  evaluated  in  the  model  during  the  loading  and 
unloading  process  according  to  the  equation 
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F 

s 


T  +  B  0 
s  s 


t  +  Pr  6 
s  s 


where  a  =  current  hydrostatic  stress 

Pr  =  maximum  hydrostatic  stress  expected  for  a  particular 
problem  (here  Pr  =  30  ksi) 


And  finally  the  shear-volumetric  coupling  parameter  k  (Equation  3.53) 
is  expressed  as 


_Sk 

k  =  K  e 
o 

where  Ko  =  0.405 

6.  =  45 
k 

Using  the  above  parameters  an  initial  prediction  of  VT6.5-1  was 
made.  The  results  of  this  prediction  are  presented  in  Figure  4.12  in 
terms  of  axial  stress  versus  axial  and  lateral  strains.  The  compari¬ 
son  is  quite  good  especially  considering  that  this  is  a  simple  endo- 
chronic  model.  A  more  important  observation  is  made  when  comparing 
the  invariants  oo  and  eo  for  the  model  and  test  results,  which  is 
presented  in  Figure  4.13.  The  ECPM  simulates  (at  least  qualita¬ 
tively)  the  unloading  phase  of  the  test,  and  in  fact  predicts  compac¬ 
tion  during  pure  deviatoric  unloading  branches.  The  significance  of 
this  feature  is  seen  by  comparing  the  results  presented  in  Fig¬ 
ure  4.13  with  similar  predictions  by  the  FEBM  as  presented  in  Fig¬ 
ure  4.11.  It  is  obvious  that  failure  to  capture  the  unloading 
response  characteristics  In  a  structural  dynamics  problem  which  might 
involve  many  unload-reload  cycles  could  result  in  an  accumulation  of 
error  that  is  unacceptable  in  structural  analysis  and  design.  The 
effects  of  altering  one  of  the  hardening  parameters  (F  )  is  shown  In 
Figures  4.14  and  4.15  where  6  is  varied  from  40  to  100.  Using  the 

H 

lower  value  of  B  is  seen  to  result  in  a  softening  behavior  while 
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the  higher  value  results  in  a  stiff er  response.  In  either  case,  the 
model  still  predicts  compaction  on  the  deviatoric  unloading  branches. 
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4.3.4  FEBM  Versus  Test  VT6.5-2 

As  discussed  in  Chapter  2,  VT6.5-2  was  a  proportional  load  path 
test  designed  to  reach  a  final  stress  state  similar  to  that  of  test 
VT6.5-1.  The  stress  path  is  presented  in  Figure  2.17  and  measured 
axial  stress  versus  axial  and  lateral  strains  are  presented  in  Fig¬ 
ure  4.16.  Post  test  measurements  indicated  a  permanent  axial  strain 
of  -0.034  in. /in.  and  a  permanent  lateral  strain  of  0.036  in. /in. 
Based  on  this  observation  plus  comparisons  with  previous  calibration 
tests  near  this  stress  region  gages  3  (axial)  and  4  (lateral)  were 
selected  to  provide  test  data  for  comparison.  Gage  5  was  ruled  out 
because  of  the  unusual  early  soft  response  then  the  stiffer  response 
followed  by  what  appears  to  be  a  partial  debonding  of  the  gage  at 
maximum  stress. 


STRAIN  in. /in. 

Figure  4.16  Actual  test  data  for  test  VT6.5-2 
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Predicted  stresses  using  the  FEBM  (associated  flow)  are  pre¬ 
sented  in  Figure  4.17.  Notice  here  that  the  model  predicts  loading 
up  to  the  failure  surface  (Figure  4.17a)  then  softening  is  predicted 
in  terms  of  stress-strain  response  (Figure  4.17(b))  while  the  stress 
path  moves  down  the  failure  surface.  This  response  is  presented  in 
terms  of  ao  versus  eo  in  Figure  4.17(c).  This  is  inconsistent 
with  the  theory  of  the  model,  in  that  there  is  no  predicted  decohe¬ 
sion  of  the  failure  surface  as  softening  is  occurring.  The  signifi¬ 
cant  affect,  on  response,  of  initial  dilation  is  seen  in 
Figure  4.17(c).  Results  from  the  non-associated  flow  version  of  the 
FEBM  are  presented  in  Figure  4.18.  Here  the  predicted  stress  path 
(Figure  4.18(a))  tends  toward  the  test  stress  path  at  a  confining 
stress  of  about  -4  ksi  then  continues  to  load  up  to  an  axial  stress 
of  about  -21  ksi  (which  is  close  to  the  maximum  measured  axial 
stress).  After  reaching  the  maximum  axial  stress,  the  stress  path 
appears  to  move  downward  along  degrading  loading  surfaces  (i.e.  the 
cohesion  is  decreasing)  until  complete  unloading  occurs.  From  Fig¬ 
ure  4.18  00,  the  model  is  seen  to  predict  softening  at  an  axial 
stress  of  about  22  ksi,  which  is  slightly  less  than  the  measured 
axial  stress  of  about  23.5  ksi.  This  is  probably  due  to  the  fact 
that  the  model  does  not  solve  for  the  exact  point  where  the  stress 
path  intersects  the  failure  surface.  The  algorithm  within  this  model 
predicts  a  stress  point  outside  the  failure  surface  then  uses  a 
return  strategy  to  get  back  on  the  failure  surface.  Once  the  model 
begins  to  predict  softening,  the  softening  predicted  is  much  greater 
than  that  measured  in  the  test.  The  model  parameters  which  most 
affect  softening  predictions  are  the  residual  crack  opening  displace¬ 
ment  (Ur) ,  and  the  ratio  of  mode  II  to  mode  I  fracture  energies 
G^/Gj  .  Small  perturbations  (eg.  ±  20X  up  to  factors  of  2  changes) 
in  the  parameter  produced  no  real  effects.  Figure  4.19  presents  the 
results  of  multiplying  the  calibrated  value  of  G^/G^  by  5.  In 
Figure  4.20  the  measured  stress  path  is  compared  with  the  failure 
surface  and  the  region  of  degrading  surfaces  due  to  predicted 
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softening.  It  is  obvious  from  the  cross  hatched  region  that  the 
loading  paths  of  VT6.5-I  and  VT6.5-2  present  very  difficult  problems 
for  constitutive  models  due  to  their  close  proximity  to  the  failure 
surface  in  the  softening  region. 


157 


4.3.5  ECPM  Versus  Test  VT6.5-2 

The  ECPM  predicted  a  similar  response  to  that  of  the  FEBM  for 
the  associated  flow  version.  The  results  are  presented  in  Fig¬ 
ure  4.2.1.  '.Then  dilation  initiates,  the  model  begins  to  predict 
softening  as  can  be  seen  in  Figure  4.21b.  By  decreasing  F^  a 
stiffer  stress  strain  response  and  higher  peak  stress  can  be 

obtained.  To  obtain  this  result  8  was  decreased  by  50%  to  8  = 

s  s 

0.3445  and  the  results  are  presented  in  Figure  4.22.  For  the  pre¬ 
peak  portion  of  the  stress  path  (Figure  4.22a)  the  model  compares 
favorably  with  test  results  up  to  peak  axial  stress.  As  loading  con¬ 
tinues  into  the  softening  region  the  model  over  reduces  the  lateral 
stresses  to  account  for  the  beginning  of  expansive  volume  change  as 
shown  in  Figure  4.22c.  As  mentioned  previously,  having  only  one 
exponential  term  for  the  response  kernels  restricts  the  model  to  per¬ 
form  well  in  only  one  stress  region.  No  real  effect  was  seen  in 
reasonable  variations  of  other  parameters.  To  look  at  the  affect  of 
variation  of  the  exponential  term  in  the  shear  kernel  a  was  varied 
from  aij  =  70  to  =  250  which  represents  changes  from  about  50% 
of  the  calibrated  value  to  almost  2  times  the  calibrated  value.  The 
results  of  these  changes  are  presented  in  Figures  4.23  and  4.24.  The 
lower  value  of  a  results  in  essentially  linear  elastic  response. 

It  is  interesting  to  note  that  changes  in  a  simply  shift  the  value 
of  intrinsic  time  at  which  the  exponential  term  saturates. 
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4.3.6  FEBM  -  Test  VT6.5-3 

This  test  is  similar  to  test  VT6.5-1  except  that  the  level  of 
confining  stress  is  much  higher  here.  The  measured  axial  stress  ver¬ 
sus  axial  and  lateral  strains  are  presented  in  Figure  4.25,  (only  one 
lateral  gage  survived  the  test) .  Based  on  pre  and  post  test  microme¬ 
ter  measurements  the  plastic  axial  strain  was  recorded  as  -0.058  and 
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Figure  4.25  Actual  test  data  for  VTA. 5-3 

the  plastic  lateral  strain  as  0,05,  resulting  in  a  plastic  volume 
change  of  +0.044.  These  values  compare  well  with  strain  gage  mea¬ 
surements.  Also,  the  two  axial  gages  show  the  same  trends  in 
response,  with  gage  5  possibly  being  in  contact  with  a  higher  per¬ 
centage  of  aggregate  surface  or  possibly  not  being  bonded  quite  as 
good  as  gage  3.  Gage  3  was  selected  for  use  iii  defining  the  axial 
response  of  the  total  specimen.  Figure  4.26,  presents  the  loading 
stress  path  for  the  test  along  with  the  failure  surface  predicted  by 


0  5000  10000  15000  20000  25000  30000 

/2  #  Sigma  X  =  / 2  #  Sigma  Y.  psi 

(a)  stress  path 


(b)  axial  stress  versus  axial  and  lateral  strains 

Figure  4.26  WES  test,  VT6.5-3,  hydrostatic  deviatoric  load  path 
at  high  confining  stress  region,  with  stress  strain  response 
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the  FEBM.  The  points  of  transition  from  hydrostatic  to  deviatoric 
branches  are  defined  on  the  stress  path  as  well  as  on  the  vertical 
stress  axial  and  lateral  strain  plots.  The  response  during  loading 
from  0-1,  to  1-2,  to  2-3  is  as  expected  and  clearly  shows  the  com¬ 
pressive  increments  of  strain  for  all  strain  gages  during  hydrostatic 
branches  (0-1  ,  2-3)  and  the  expansive  increments  in  lateral  strains 
but  net  compaction  response  during  the  deviatoric  branch  (1-2) .  Near 
the  end  of  branch  3-4  the  loading  stress  path  intersects  the  failure 
surface  and  the  corresponding  stress  strain  response  (Figure  4.26b) 
shows  the  general  hardening  characteristics  as  previously  observed  in 
calibration  tests.  At  point  4  the  total  axial  strain  is  about 
-0.04  in. /in  which  indicates  a  significant  amount  of  plastic  yielding 
has  occurred.  The  hydrostatic  branch  (4-5)  is  seen  to  produce  a  very 
stiff  response  consistent  with  what  is  often  referred  to  as  locking 
up  in  soil  and  concrete  at  high  hydrostatic  stress  levels.  Finally 
the  last  branch  of  the  deviatoric  load  path  is  (5-6)  and  the  stress 
strain  response  is  as  would  be  expected.  At  point  6  and  during  the 
reversal  of  loading  process  there  is  some  near  plastic  flow  with  a 
very  large  increase  in  lateral  (expansive)  strain.  It  should  be  men¬ 
tioned  that  these  strain  gages  are  high  elongation  gages  and  rated 
linear  up  to  strains  on  the  order  of  10%.  A  plot  of  octahedral 
stress  versus  strain  based  on  test  data  is  presented  in  Figure  4.27. 
The  response  up  to  point  4  is  similar  to  that  measured  in  test 
VT6.5-1  when  compaction  is  measured  along  deviatoric  loading 
branches.  However  after  point  4  (which  intersects  the  failure  sur¬ 
face)  the  measured  response  is  dramatically  different  than  that  of 
test  VT6.5-1.  The  overall  response  of  the  specimen  can  best  be  seen 
from  the  octahedral  stresses  and  strains  presented  in  Figure  4.27. 

The  response  along  branch  4-5  indicates  essentially  infinite  bulk 
modulus,  which  if  taken  literally,  implies  that  all  the  voids  in  the 
specimen  have  been  closed.  Along  branch  5-6  the  specimen  begins  to 
dilate  unlike  the  compactive  response  seen  along  similar  branches  in 
test  VT6.5-1.  The  plastic  flow  lateral  strains  occurring  while  the 
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Figure  4.27  VT6.5-3,  octahedral  stress  vs.  strain 


axial  stress  is  near  point  6  are  very  large  and  expansive.  As 
unloading  occurs  along  6'-5'  there  is  still  more  dilation.  Branch 
5’-4'  (hydrostatic  unload)  is  accompanied  by  expansive  strains. 

There  is  some  very  small  compaction  during  the  unloading  branch  4' -3' 
while  the  unloading  along  3r-2f  is  at  a  stiffer  slope.  Again  some 
small  amount  of  compaction  occurs  during  2*-l*  and  a  reasonable 
response  is  observed  from  1'  to  0. 

Results  of  predictions  from  the  associated  flow  version  of  the 
FEBM  are  presented  in  Figure  4.28.  Due  to  reasons  previously  dis¬ 
cussed  the  associated  flow  model  predicts  a  much  stiffer  response 
than  test  measurements  indicate.  In  fact  at  about  point  3  the  model 
is  predicting  stresses  near  the  failure  surface.  Therefore ,  the 
failure  surface  is  intersected  much  too  early  in  the  loading  history 
and  since  strains  continue  to  increase  the  predicted  stress  point 
mnvpe  failure  surface  to  an  extremely  high  axial  stress  of 

about  84  ksi.  It  is  difficult  to  say  exactly  what  is  occurring  in 
the  model  at  this  point. 
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Lateral  strains  are  becoming  very  large  during  branch  3-4,  which 
has  a  significant  affect  on  the  plastic  strain  direction.  The 
response  of  the  FEBM  is  very  sensitive  to  plastic  strain  increment 
direction.  The  reason  for  this  is  that  the  plastic  stress  corrector 
is  not  generally  coaxial  with  the  plastic  strain  increment.  When  the 
strain  is  expansive  (dilatant)  the  stress  corrector  will  drift  toward 
the  compressive  side.  Conversely  when  the  strain  is  compactive  the 
stress  corrector  tends  toward  the  expansive  side.  During  the  test 
the  lateral  strains  appear  to  dominate  the  response  and  the  model 
intersects  the  failure  surface  and  the  stress  point  begins  to  slide 
down  the  failure  surface  as  shown  in  Figure  4.28a.  Figure  4.29 
presents  the  results  of  reducing  E  to  E  =  5  x  10^  psi  and  increas¬ 
ing  the  ductility  parameter  by  a  factor  of  2.  The  comparison  is 
better  but  the  problem  at  failure  still  occurs. 

Results  using  the  ron-associative  flow  version  of  the  FEBM  are 
presented  in  Figure  4.30  with  E  =  5  x  10^  psi.  As  shown  in  Fig¬ 
ure  4.30  the  comparison  with  test  results  is  similar  to  that  with  the 
associated  flow  version.  No  improvement  in  prediction  were  obtained 
by  varying  hardening  parameters  or  elastic  moduli.  Figure  4.31  shows 
the  effects  of  increasing  the  ductility  parameter.  The  main  problem 
that  this  test  poses  for  the  FEBM  is  that  when  the  model  predicts 
intersection  of  the  stress  path  with  the  failure  surface  near  point  4 
(which  actually  occurs  in  the  test)  the  hardening  parameter  has 
reached  the  maximum  value  of  1  and  cannot  decrease  during  the 
remainder  of  the  test.  Therefore,  the  model  must  either  predict  a 
purely  elastic  response  if  the  stress  point  moves  inside  the  failure 
surface  or  predict  a  response  associated  with  a  neutral  type  loading 
as  the  stress  point  slides  along  the  failure  surface.  Another  point 
to  he  made  is  that  the  model  may  have  compared  better  if  gage  3 
(Figure  4.23)  had  been  selected  as  the  axial  strain  for  comparison. 
However,  this  is  doubtful  since  the  relative  magnitude  of  the  lateral 
strains  would  be  even  larger  for  this  case. 


(b)  normal  stress  6train 


(c)  octahedral  stress  strain 
Figure  4.30  FEBM-VT6 . 5-3 ,  non-associated  flow,  E=5.0E6 
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4.3.7  ECPM  Versus  VT6.5-3 

Results  from  the  initial  calibration  of  the  ECPM  are  compared 
with  test  results  for  VT6.5-3  in  Figure  4.32.  Model  predictions  com¬ 
pare  reasonably  well  up  to  branch  3-4  where  the  model  tends  to  over¬ 
predict  the  response.  At  an  axial  strain  of  about  4%  the  model 
begins  to  predict  softening  until  the  end  of  the  test.  The  ECPM  is 
also  reaching  some  type  of  limiting  surface  at  an  axial  strain  of 
about  41%  and  sliding  down  the  surface  throughout  the  rest  of  the 
test.  A  remarkable  similarity  in  the  predicted  results  of  the  two 
models  for  this  test  can  be  seen  by  comparing  Figures  4.30  (FEBM)  and 
4.32  (ECPM).  The  results  of  reducing  the  hardening  parameter  F 

H 

by  varying  over  a  broad  range  of  values  is  shown  in  Figure  4.33. 

From  Figure  4.33  it  is  seen  that  F^  essentially  does  not  impact 
response  predictions  until  the  stress  path  reaches  branch  3-4.  The 
effects  of  the  shear  hardening  parameter  Fs  was  studied  by  varying 
8s  from  half  the  calibrated  value  to  1.5  times  the  calibrated  \alue. 
These  results  are  presented  in  Figure  4.34.  The  sensitivity  of 
response  to  variation  in  C  and  F  are  shown  in  Figure  4.35. 


(b)  normal  stress  strain 


(o')  octahedral  stress  strain 
Figure  4.32  ECPM-VT6 . 5-3 ,  initial  calibration 
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4.3.8  FEBM  -  VT6.5-4 

As  discussed  in  Chapter  2,  VT6.5-4  was  a  proportional  load  path 
test  designed  to  reach  a  final  stress  state  similar  to  that  of 
VT6.5-3.  The  stress  path  for  this  test  is  presented  in  Figure  2.19. 
Measured  axial  stress  versus  axial  and  lateral  strains  are  presented 
in  Figure  4.36.  From  micrometer  measurements  the  permanent  axial 
strain  was  measured  at  -0.O66  in/in. ,  and  lateral  strain  at  +0.057. 
Predicted  results  using  the  associated  flow  version  of  the  FEBM  are 
pre5,-ented  in  figure  4.37  for  the  initial  calibration.  As  shown  in 
the  figure  the  model  intersects  the  failure  surface  at  relatively 
small  axial  strair  and  then  slides  down  the  failure  surface  as  soft¬ 
ening  is  predicted.  The  model  does  reproduce  the  unload-reload  cycle 
which  occurs  (for  the  predicted  response)  during  softening  (see  Fig¬ 
ure  4.37b).  For  the  associated  flow  version  the  stress  point  slides 
down  the  failure  surface  during  softening  so  that  no  damage  (or  deco¬ 
hesion)  is  predicted  as  was  predicted  in  the  true  softening  test 
VT6.5-2.  Predictions  for  the  non  associated  flow  version,  initial 
calibration,  are  presented  in  Figure  4.38.  Here  the  model  again  pre¬ 
dicts  to  stiff  a  response  until  the  failure  surface  is  reached  and 
then  softening  response  is  predicted  as  the  stress  point  slides  down 
the  failure  surface.  For  this  test  no  affects  were  observed  by  vary¬ 
ing  any  of  the  ductility  parameters.  Figure  4.39  shows  the  small 
effect  of  changing  y  from  0.5  times  the  calibrated  value  to  3 
times  the  calibrated  value.  This  is  the  only  test  where  the  response 
of  the  model  was  so  little  sensitive  to  parameter  variations.  Once 
again  the  breakdown  in  the  model  (here  predicting  softening)  occurs 
when  the  octahedral  strains  begin  to  dilate  (see  Figure  4.38(c)). 


(c)  octahedral  stress  strain 

Figure  4.38  FEBM-VT6 . 5-4 ,  non-associated  flow,  initial  calibration 
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4.3.9  ECPM  -  VT6.5-4 


Predictions  from  the  ECPM  are  presented  in  Figure  4.40  for  the 
initial  calibration,  where  it  is  seen  that  the  model  significantly 
under  predicts  the  response.  The  effects  of  reducing  8g  and 
increasing  8  are  shown  in  Figure  4.41.  As  was  the  case  for  the 
FEBM  reasonable  variations  in  model  parameters  had  little  affect  on 
predicted  response  and  in  no  case  would  the  model  predict  response 
beyond  the  unload-reload  cycle  (which  occurs  just  after  dilation 
initiates) . 
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4.3.10  FF.BM  and  ECPM  -  VT6.5-5,  VT2-1 

These  tests  were  not  conducted  under  the  same  test  program  as 
VT6.5-1  through  VT6.5-4,  and  only  digitized  records  of  stresses  and 
strains  were  obtained  for  the  tests.  Only  the  initial  calibration 
will  be  used  for  the  FEBK  and  ECPM  in  predicting  these  tests  results. 
Calibration  parameters  for  the  lower  strength  concrete  discussed  here 
and  in  sections  4.4  and  4.5  are  presented  in  Wiliam  (3)  for  the  FEBM 
and  Valanis  (4)  for  the  ECPM. 

Predicted  responses  for  test  VT6.5-5  from  the  FEBM  are  presented 
iri  Figure  4.42  and  4.43  while  Figure  4.44  presents  results  from  the 
ECPM.  The  associated  flow  version  over  predicts  the  axial  stress 
then  fails  abruptly.  The  non-associated  version  also  over  predicts 
toe  axial  stress  but  then  softens.  Also  Figure  4.43a  shows  fairly 
good  agreement  between  model  and  test  results  in  terms  of  oo  vs  to 
(at  least  qualitatively) .  The  ECPM  over  predicts  the  axial  stiess 
hue  here,  it  seems  that  the  ECPM  is  simply  too  stiff  in  hydrostatic 
response  by  a  constant.  Once  the  peak  stress  is  reached  the  model 
predicts  abrupt  softening  then  essentially  plastic  flow  at  an  axial 
stress  of  about  o z  =  18.5  ksi.  Qualitative  comparisons  with  the 
octahedral  stress  strain  plot  in  Figure  4.44a  is  quite  good. 

Predictions  from  the  FEBM  and  the  ECPM  are  compared  with  test 
results  from  test  VT2-1  in  Figures  4.45  and  4.46.  There  was  no  dif¬ 
ference  in  the  associated  flow  and  non-associated  flow  version  for 
this  test  so  only  the  results  from  the  associated  flow  version  are 
presented . 
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4.4  FEBM  and  ECPM  -  VT4-1 ,  VT4-2 ,  VT4-3 

These  tests  were  conducted  at  the  University  of  Colorado  in  the 
CU  tour-inch  cubical  cell.  Strains  are  measured  in  this  device 
through  the  use  of  displacement  gages  which  are  in  contact  with  a 
point  at  the  center  of  the  specimen  face.  The  associated  and  non- 
associated  flow  versions  of  the  FEBM  produced  essentially  the  same 
results  for  test  VT4-1  .  Also  results  from  the  ECFM  were  very  similar 
to  the  FEBM  and  therefore  only  the  results  from  the  associated  flow 
FEBM  are  presented  in  Figure  4.47.  Unloading  data  was  not  provided 
and  therefore  little  insight  into  material  response  features  can  be 
gained  from  this  test. 

Furthermore,  both  versions  of  the  FEBM  and  the  ECPM  predicted 
similar  results  for  VT4-2  and  therefore  only,  the  FEBM  results  are 
presented  in  Figure  4.48  Studying  the  test  data  (VT4-2)  and  model 
predictions  in  more  detail  reveals  some  inconsistencies.  First  just 
considering  the  measured  strains,  after  point  A  is  reached  in  Fig¬ 
ure  4.48(a)  ex  =  ey  =  constant  while  c z  increases  to  a  maximum 
value  of  fz  =  -0.0135  in/in.  This  branch  (ie  AB  in  Figure  4.48(a)) 
is  then  equivalent  to  a  uniaxial  strain  test  (see  Section  2.5). 

Based  on  these  observations,  the  model  is  predicting  the  stresses 
consistent  with  this  strain  path.  In  other  words  when  the  model 
receives  constant  strain  input  it  predicts  a  oz  ez  path  which  is 
similar  to  the  uniaxial  strain  response  as  shown  in  Figure  2.8. 

Based  on  these  observations  it  is  concluded  that  there  is  some  error 
in  the  test  data  beyond  point  A.  One  would  expect  a  response  like 
that  measured  In  test  VT6.5-5,  where  strain  gages  were  used  to  record 
strains . 

There  was  little  difference  in  the  FEBM  and  ECPM  in  predicting 
the  response  for  VT4-3  (Figure  4.49)  which  is  a  stair  step 
hydrostatic-deviatoric  test  similar  to  VT6.5-1,  and  VT6.5-3,  except 
that  the  unloading  branches  were  not  conducted.  The  improvement  in 
r'-'Tt-l  test  comparisons  using  the  mixed  control  option  can  be  seen  in 
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4.5  FFRM  and  EC  PM  -  VT5-1  ,  VT5-2  ,  VT5-3 

These  tests  were  conducted  in  the  Eindhoven  cubical  cell  as  dis¬ 
cussed  in  Chanter  ?.  The  associated  flow  version  FEBK  failed  to  con¬ 
verge  on  a  solution  after  the  early  portion  of  the  strain  path  for 
VTS-1 .  Tiiis  is  probablv  due  to  the  severe  strain  softening  measured 
in  the  test.  Results  for  the  FEBM  and  the  F.CPM  are  presented  in  Fig¬ 
ure  4.12.  Except  for  the  maximum  strength  prediction  and  stiffer 
predicted  response  the  FEBM  performed  quite  well  for  this  test  espe¬ 
cial  lv  in  the  softening  region.  Model  test  comparisons  for  VT5-2  and 
VT5-3  are  presented  in  Figures  4.53  and  4.54  respectively.  Except 
for  the  consistentlv  high  predicted  maximum  stress  and  stiffer  pre¬ 
dicted  prepeak  response  all  che  models  performed  quite  well  in  these 
tests.  The  models  certainly  appear  to  capture  this  overall  softening 
response  (ie  Tests  VT5-?  and  /T5-3)  much  better  than  in  test  VT6.5-2. 
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Chapter  5 

Conclusions /Recommends  t ions 

5  . 1  Summary 

An  effective  methodology  has  been  developed  for  evaluating  the 
predictive  capability  of  potential  constitutive  models.  This  method¬ 
ology  consists  of  designing  and  conducting  a  series  of  tests  which 
are  used  to  calibrate  the  model,  designing  and  conducting  a  second 
series  of  tests  which  are  used  to  verify  the  mode],  and  finally  exer¬ 
cising  the  calibrated  model  along  the  strain  or  stress  paths  of  the 
verification  tests  so  that  predictions  and  test  results  can  be  com¬ 
pared.  The  stress  and  strain  paths  of  the  verification  tests  should 
be  consistent  with  all  design  load  paths  for  critical  regions  in  the 
structure  of  interest.  Also,  stress  and  strain  paths,  which  demon¬ 
strate  key  complex  response  features  of  the  material,  should  be 
included  in  the  series  of  verification  tests.  Model  predictions 
should  be  carefully  compared  with  test  measurements  for  stress  path, 
strain  path,  stresses  versus  strains,  octahedral  stresses  versus 
octahedral  strains,  etc.  It  is  critical  that  invariant  quantities 
such  as  the  octahedral  components  be  studied  since  subtle  features 
like  shear  compaction  on  pure  deviatoric  load  paths  are  not  so 
noticeable  when  looking  at  normal  stress  and  strain  plots. 

Two  of  the  most  recent  constitutive  models  for  concrete  (the 
Fracture  Energy  Based  Model  FEBM,  and  the  Endochronic  Concrete  Plas¬ 
ticity  Model  ECPM)  were  then  tested  using  the  proposed  methodology. 

A  number  of  tests  were  specifically  designed  and  conducted  for  the 
evaluation.  In  addition  data  from  tests  conducted  at  the  University 
of  Eindhoven,  The  Netherlands  and  the  University  of  Colorado,  Boulder 
were  also  used.  The  results  indicated  that  the  models  would  indeed 
reproduce  some  of  the  features  of  the  response  but  failed  in  others. 


5.2  Conclusions 

A  number  of  conclusions  can  be  drawn  from  the  tests  conducted 
and  the  evaluation  of  the  models  in  relation  to  compaction  under  pure 
deviatoric  stress  paths,  isotropic  versus  kinematic  hardening,  strain 
softening  and  associated  versus  non-associated  flow  rules. 

5.2.1  Shear  Volume  Response  (Compaction) 

A  critical  feature  of  concrete  response  was  demonstrated  in  test 
VT6.5-1  where  compaction  was  observed  during  both  loading  and  unload¬ 
ing  along  pure  deviatoric  stress  paths.  An  explanation  for  this 
phenomenon  is  presented  in  Figure  5.1.  Here  at  the  end  of  a  hydro¬ 
static  branch  (either  compactive  or  expansive)  the  material  is  at 
some  volume  Vi.  As  shear  is  applied  (with  hydrostatic  stress  held 
constant)  fracture  or  crushing  will  occur  at  points  where  stress  con¬ 
centrations  exist.  These  points  often  will  be  near  voids,  small 
cracks  or  small  gaps  between  aggregate  surfaces  and  cement  and  the 
tendency  will  be  to  close  or  fill  these  openings  resulting  in  a  final 
volume  (V^)  which  is  less  than  V^.  As  the  number  of  deviatoric 
cycles  and  or  the  level  of  hydrostatic  stress  increases  the  volume 


Hydrostatic 


V. 

i 


V 


f 


Figure  5.1  Shear  compaction  concept 
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change  (ie  -  V^)  should  approach  zero.  This  was  observed  in  test 
VT6.5-3  which  was  conducted  at  a  much  higher  hydrostatic  stress. 

5.2.2  Conventional  Isotropic  Hardening  Plasticity 

It  was  pointed  out  in  Chapter  4  that  the  "isotropic"  hardening 
model  used  in  the  FEBM,  results  in  failure  of  the  model  to  predict 
important  material  response  features  (i.e.  compaction  along  certain 
deviatoric  paths).  Also,  since  the  leading  surface  does  not  contract 
the  model  only  predicts  elastic  response,  after  the  maximum  strength 
surface  is  reached,  and  the  stress  point  moves  back  into  the  assumed 
"elastic"  region.  In  fact  this  is  also  the  case  for  cap  models  which 
make  use  of  isotropic  hardening  only.  While  this  assumption  is  rea¬ 
sonable  in  many  cases,  for  metal  plasticity,  it  is  inappropriate  for 
concrete  since  significant  plastic  response  can  occur  after  the  maxi¬ 
mum  strength  surface  is  intersected  and  loading  continues,  in  the 
direction  toward  the  interior  of  the  loading  surface,  which  was  pre¬ 
cisely  the  case  in  test  VT6.5-1  and  VT6.5-3. 

5.2.3  Strain  Softening 

The  argument  that  strain  softening  is  a  structural  as  opposed  to 
a  material  property  is  usually  based  on  observations  made  during  post 
test  examinations  of  test  specimen  geometry.  Typical  examples  of 
these  observations  for  unconfined  compression  tests  and  confined 
tests  of  rocks  are  that  large  cracks,  barreling  or  shear  handing 
usually  occur  during  the  test  which  violate  the  assumptions  of  homo¬ 
geneity  for  the  test  specimen  and  the  internal  stresses.  Typical 
post  test  views  of  these  test,  specimens  are  shown  in  Figure  5.2. 
However  when  low  confinement  compression  tests  are  conducted  (e.g. 
test  VT6.5-2)  softening  can  be  observed  in  test  measurements.  When 
typical  test  specimens  are  cut  open  and  viewed  with  the  naked  eye  no 
discernible  difference  is  observed  between  pre  test  and  post  test 
conditions.  However  when  a  small  (say  1  inch  square)  window  of  the 
cut  specimen  is  examined  under  magnification  (Figure  5.3)  the  differ¬ 
ences  in  pre  and  post  test  internal  structure  of  the  specimen  can  be 
clearly  seen.  In  Figure  5.3(b)  the  fracture  of  aggregates,  void 


Figure  5.2.  Post  test  views,  barreling,  shear  banding. 

closure,  and  separation  of  aggregate  from  cement  paste  at  some  points 
can  be  seen.  However,  the  post  test  specimen  is  still  intact  and 
similar  specimens  have  been  retested  in  unconfined  compression  and 
seen  to  still  possess  from  75%  to  100%  of  their  initial  strength  and 
elastic  moduli.  The  conclusion  drawn  here  is  that  the  material  of 
Figure  5.3b  probably  satisfies  the  assumptions  for  a  macroscopic  con¬ 
tinuum  just  as  well  as  the  material  in  Figure  5.3a.  Furthermore,  if 
a  well  planned  carefully  executed  series  of  calibration  tests  are 
conducted  on  the  material  of  Figure  5.3b,  hardening  and  softening 
parameters  can  be  determined  to  effectively  calibrate  a  rational  con¬ 
stitutive  model  which  can  predict  the  softening  response  features. 
This  is  not  to  say  that  problems  will  not  arise  in  terms  of  unique¬ 
ness  and  stability  of  solutions.  The  main  point  is  that  material  at 
critical  regions  in  certain  concrete  structures  can  exhibit  softening 
and  there  will  be  a  local  redistribution  of  stresses  and  external 
loads.  For  the  region  where  softening  is  occurring,  continued  load¬ 
ing  will  appear  to  be  in  displacement  control.  Constitutive  models 
which  are  rationally  developed  to  simulate  softening  can  be  effec¬ 
tively  used  to  predict  the  stress  redistribution  in  these  structures. 
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Figure  5.3.  6X  magnification,  pre  and  post  test, 


5.2.4  hon-Associated  Flow 

From  an  experimental  standpoint  results  presented  in  this  study 
fe.g.  Figure  4,8)  clearly  indicate  that  concrete  response  does  not 
necessarily  satisfy  conditions  required  for  associated  flow  plastic¬ 
ity.  The  non-associated  flow  version  of  the  FEBM  generally  improved 
the  predictive  capability  of  the  model  compared  with  the  associated 
flow  version. 

5.2.5  Endochronic  Model 

The  main  conclusion  drawn  for  the  ECPM  is  that  the  model  seems 
to  be  based  on  just  as  sound  and  rational  assumptions  as  plasticity 
models  and  is  not  based  on  a  variety  of  curve  fits.  The  capability 
of  the  model  to  predict  compaction  on  all  deviatoric  branches  in  test 
VT6.5-1  is  very  impressive  for  such  a  simple  model.  The  simple  model 
presented  here  does  not  take  advantage  of  the  concept  of  adding 
together  several  endochronic  elements  in  parallel  so  that  a  wide 
range  of  material  responses  can  be  modeled.  This  is  one  of  the 
strong  features  of  the  endochronic  theory.  However,  the  ECPM  did  not 
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perform  well  in  the  softening  region  as  shown  in  Figure  4.22.  This 
is  partly  due  to  the  use  of  onl}  one  endochronic  element  and  partly 
due  to  the  fact  that  the  model  is  not  calibrated  to  predict  softening 
response . 

5.3  Recommendations 

A  considerable  amount  of  research  work  remains  to  be  done  in  the 
area  of  constitutive  modeling  of  concrete,  rock,  or  soils.  The 
methodology  presented  in  this  work  should  serve  as  a  basis  for  cali¬ 
bration  and  evaluation  of  present  and  future  models.  Since  the  pres¬ 
ent  study  was  limited  tc  the  evaluation  of  only  two  models  it  appears 
that  some  of  the  other  models  available  or  present  should  be  sub¬ 
jected  to  the  same  process.  Some  recommendations  on  how  to  apply  the 
proposed  methodology  and  specific  areas  of  behavior  that  need  further 
examination  are  discussed  next, 

5.3.1  Application  of  the  Methodology 

When  the  methodology  developed  in  this  study  is  used  to  evaluate 
constitutive  models,  calibration  tests  and  verification  tests  should 
be  repeated  to  demonstrate  uniformity  and  consistency  of  measurements 
as  well  as  variation  in  test  results  due  to  experimental  error.  Mod¬ 
els  should  be  driven  under  strain  control,  stress  control,  and  mixed 
control  for  a  full  evaluation. 

5.3.2  Shear  Volume  Response 

A  series  of  hydrostatic-aeviatoric  tests  should  be  conducted  at 
different  deviatorlc  sections  to  determine  the  effects  of  hydrostatic 
stress  level  on  shear  compaction.  Also,  similar  tests  where  the 
deviatoric  path  is  cycled  several  times  should  be  conducted  to  deter¬ 
mine  this  effect  on  limiting  values  for  volume  compaction. 

5.3.3  Hardening 

It  has  been  shown  in  this  study  that  conventional  isotropic 
hardening  concepts  alone  cannot  capture  the  nonlinear  material 
response  of  concrete  when  the  stress  point  reaches  the  maximum 
strength  surface  then  moves  back  into  the  "elastic  region."  This  is 
the  case  for  both  associated  and  non-associated  flow.  The  simple 
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endochronic  model  did  reasonably  well  in  predicting  this  response  and 
should  be  further  studied  and  developed  in  this  area.  Also,  kine¬ 
matic  hardening  concepts  should  be  further  studied  and  developed  in 
this  context  to  determine  their  capabilities.  Further  testing  should 
be  conducted  where  the  stress  point  moves  out  to  the  maximum  strength 
surface  (along  different  paths)  then  moves  back  into  the  elastic 
region  (along  different  paths).  Measured  plastic  strain  increments 
should  be  plotted,  along  these  paths,  which  will  indicate  the  orien¬ 
tation  and  position  of  a  kinematic  surface  (if  one  exists).  These 
tests  should  be  performed  along  load  paths  in  the  rendulic  plane 
(Figure  5.4a)  as  well  as  along  fully  three  dimensional  load  paths 
(Figure  5 . 4b)  , 

5.3.4  Strain  Softening 

A  carefully  planned  test  program  should  be  conducted  to  study 
the  issue  of  strain  softening.  This  test  program  should  include  the 
following  considerations: 

(a)  Careful  preparation  of  test  specimen  to  insure  uniform 
consolidat ion . 

(b)  Minimize  friction  on  loaded  surfaces  of  the  specimen. 

(c)  Repeated  softening  tests  at  different  rates  of  loading. 

(d)  Careful  dissecting  of  specimen  for  internal  microscopic 
examination. 
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